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INTRODUCTION 


The  energy  in  the  exhaust  of  a  jet  engine  is  immense.  To  illus¬ 
trate,  with  an  assumption  of  only  a  25%  conversion  efficiency,  the 
thermal  energy  in  the  exhaust  of  a  J79  engine  running  at  intermediate 
power  is  tantamount  to  the  energy  required  to  generate  about  280,000  lb/hr 
of  150-psi  steam.  To  provide  a  perspective  for  this  magnitude,  the 
average  peak  steam  consumption  of  the  entire  Naval  Air  Station  (NAS.), 
Miramar,  Calif.,  on  a  winter  day,  is  about  35,000  lb/hr. 

New  jet  engines  and  engines  undergoing  maintenance  or  repair  are 
examined  in  test  cells*  prior  to  operational  use.  In  these  test  facil¬ 
ities,  the  jet  energy  is  "wasted";  the  jet  exhaust  is  diluted  and  dis¬ 
charged  to  the  atmosphere.  Yet,  test  cells,  by  their  confining  geom¬ 
etry,  provide  an  opportunity  and,  if  waterwalls  are  installed  or  if  heat 
exchangers  are  located  somewhere  in  the  path  of  the  jet,  possibly  a 
means  of  obtaining  part  of  a  vast  energy  supply. 

Operation  of  Test  Cells 


Jet  engine  test  facilities  have  multiple  functions:  to  supply  a 
clean,  distortion-free  airflow  to  the  engine  inlet;  to  attenuate  noise 
and  chemical  pollution  by  the  jet;  to  decrease  the  kinetic  energy  of  the 
jet  exhaust;  and  to  cool  the  exhaust  sufficiently  before  dumping  either 
to  the  atmosphere  or  to  acoustical  or  environmental  treatments. 

Although  the  physics  of  test  cells  is  extremely  complex,  their 
concept  is  simple;  the  jet  exhaust  entrains  the  cooler  surrounding  air, 
losing  its  energy  to  this  secondary  flow,  which  it  carries  along  through 
the  augmenter.  This  concept  is  shown  schematically  on  Figure  1. 

The  actual  mass  flow  rate  of  the  augmentation  airflow  is  important 
to  the  functions  of  the  test  cell.  Too  little  augmentation  airflow: 

(1)  allows  the  exhaust  gas  to  recirculate  back  around  to  the  engine 
inlet,  resulting  in  hot  spots  and  performance  degradation;  (2)  allows 
the  test  cell  to  heat  up;  and  (3)  also  results  in  an  excess  density  of 
visible  emittants,  possibly  violating  Ringelmann  number  restrictions. 

Too  much  air  entrainment  induces  excessive  pressure  gradients  along  and 
across  the  test  cell  augmenter,  creating  structural  problems.  Too  much 
augmentation  airflow  can  also  cause  errors  in  engine  thrust  measurements 
(Ref  1). 


*No  distinction  will  be  made  between  hush  houses,  where  the  engine  remains 
mounted  on  the  aircraft,  and  test  cells,  where  the  engine  is  tested  off 
the  aircraft. 
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Scope 


For  several  reasons,  very  little  thought  has  been  given  to  test 
cell  energy  recovery.*  First,  jet  engines  are  rarely  tested  for  long 
durations.  A  few  hours  per  day  of  test  cell  usage  is  maximum.  Second, 
the  dilution  of  the  jet  exhaust  with  the  augmentation  air  decreases 
temperature  gradients  and,  therefore,  the  potential  for  heat  transfer 
from  the  exhaust.  Even  a  25%  energy  conversion  efficiency  is  unlikely 
with  most  test  cells.  Finally,  there  are  side  effects  to  energy  recovery 
that  affect  the  ability  of  a  test  cell  to  perforin  its  functions. 

Waterwalls  will  decrease  jet  exhaust  temperatures,  beginning  exactly 
at  the  nozzle.  The  environmental  ramifications  are,  generally,  bad.  If 
combustion  was  complete,  lower  gas  temperatures  would  decrease  equilib¬ 
rium  concentrations  of  both  NO  and  CO.  At  the  higher  power  settings, 
however,  combustion  continues  out  into  the  augmenter.  (Afterburning 
jets  exhaust  as  much  as  10%  of  their  fuel  (Kef  1).)  Premature  cooling 
of  the  exhaust  will  slow  any  combustion  reactions  still  taking  place. 
Unburned  hydrocarbons,  carbon  monoxide,  even  raw  fuel,  could  be  carried 
into  the  atmosphere.  The  effect  on  particulates  is  uncertain;  it  can  be 
presumed  that  changing  the  combustion  chemistry  will  change  the  parti¬ 
culate  concentration. 

Lowering  gas  temperatures  also  changes  the  flow  characteristics  of 
the  jet,  changing  kinetic  energy,  noise  generation,  and  augmenter  pres¬ 
sure  gradients. 

Some  side  effects  of  augmenter  energy  recovery  are  beneficial. 
Waterwalls  will  eliminate  burnout  problems  of  the  liners  and  acoustic 
pads.  Waterwalls  will  also  decrease  structural  thermal  shock  by  damping 
out  temperature  fluctuations  at  startup,  shutdown,  and  while  changing 
power  settings. 

Convection,  water  tube  heat  exchangers  located  close  to  the  nozzle 
will  introduce  the  same  problems  as  waterwalls.  In  addition,  the  pres¬ 
sure  gradient  across  the  device  could  induce  engine  back  pressure  prob¬ 
lems.  Convection  heat  exchangers  located  further  down  the  augmenter 
tube  will  have  little  effect  on  augmenter  gas  temperatures,  but  high 
engine  back  pressures  still  have  to  be  faced. 

Regardless,  such  a  large  energy  source  should  be  investigated. 

This  report  summarizes  the  first  step:  an  analytical  study  of  the 
economics  of  test  cell  energy  recovery.  The  study  is  limited  to  thermal 
energy  recovery;  other  forms  of  energy,  such  as  kinetic  energy,  are  not 
included.  Side  effects  are  also  neglected.  The  analyses  of  the  side 
effects  will  be  lengthy;  if  test  cell  energy  recovery  is  not  economi¬ 
cally  feasible,  an  examination  of  side  effects  is  moot. 

Both  steam  and  electricity  are  considered  as  products  of  energy 
recovery-steam  directly  and  electricity  off  a  Rankine  cycle.  Four  pos¬ 
sible  heat  exchanger**  configurations  are  analyzed  (see  Figure  1): 


*FFV,  the  Swedish  aerospace  company,  has  done  some  work  in  this  area, 
including  the  construction  of  a  small  energy  recovery  system. 

**The  terms  "heat  exchanger"  and  "boiler"  will  be  used  interchangeably. 
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1.  Convection,  water  tube  heat  exchanger  located  in  the  exhaust 
path  at  the  augmenter  exit 

2.  Waterwalls  located  near  the  augmenter  inlet  where  the  primary 
mode  of  heat  transfer  would  be  radiation 

3.  Waterwalls  located  near  the  augmenter  exit  (for  this  configu¬ 
ration,  the  primary  mode  of  heat  transfer  would  probably  be 
convect ion) 

A.  Convection,  water  tube  heat  exchanger  located  near  the  nozzle 
of  the  engine 

The  plan-of-attack  is  to  select  several  test  cells  in  the  San  Diego 
area,  evaluate  their  energy  recovery  potential,  and  then  generalize  the 
results,  reinforced  with  parametric  studies. 


TEST  CELL  SIMULATION 

The  major  burden  is  determining  how  much  heat  can  be  transferred 
from  the  jet  exhaust  to  waterwalls  or  to  a  water  tube  heat  exchanger. 

Once  determined,  steam  generation  and,  subsequently,  the  economics  of 
energy  recovery  can  easily  be  computed.  Sizes,  shapes,  wall  construc¬ 
tions,  and  flow  rates  of  test  cells  are  almost  limitless.  The  obvious 
approach  to  the  analyses  is  to  formulate  and  utilize  a  flexible  mathe¬ 
matical  model. 

The  complexities  in  modelling  test  cells  are  manifest.  Re-examine 
Figure  1  and  consider  the  following  augmenter  characteristics:  multiple, 
nonconcentric  flows,  each  at  a  different  velocity;  turbulence;  flow 
separation;  all  three  modes  of  heat  transfer,  simultaneously;  extreme 
changes  in  temperatures  and,  thus,  variable  properties;  and  nonsymmetric , 
three  dimensional  geometries.  In  addition,  some  tests  are  of  short 
duration ,  resulting  in  a  highly  transient  situation.  These  complexities 
lead  to  a  transient  nonlinear  problem  with  coupled  energy,  mass,  and 
momentum  transfer.  Numerical  techniques  are  necessary. 

Even  so,  to  make  a  test  cell  simulation  attainable,  assumptions 
must  be  made.  The  first  assumption  is  to  restrict  the  analyses  to  the 
augmenter  duct.  Jet  and  augmentation  flows  are  treated  as  boundary 
conditions.  For  this  present  work,  the  emphasis  is  on  heat  transfer. 

The  modelling  tactics  will  be  to  describe  fluid  flow  empirically  but  use 
first  principles  to  determine  heat  transfer  rates.  With  this  in  mind, 
the  following  additional  assumptions  are  made: 

1.  Flows  are  one-dimensional,  passing  through  a  cylindrical 
augmenter 

2.  Steady  state,  steady  flow  is  stipulated 

3.  Both  the  jet  and  augmentation  air  are  ideal  gases 


3 


a. 


Flow  within  the  jet  and  within  the  augmentation  air  can  be 
considered  turbulent  and  perfectly  mixed 


5.  Pressure  gradients  along  the  augmenter  are  negligible  compared 
with  temperature  gradients  and  can  be  neglected 

6.  The  jet  nozzle  is  located  exactly  at  the  inlet  to  the  augmenter 
tube 

7.  The  augmenter  walls  are  black,  the  jet  is  gray,  and  the  aug¬ 
mentation  air  is  transparent  to  thermal  radiation 

8.  Axial  radiation  along  the  jet  has  a  negligible  effect  on 
temperatures;  all  radiation  is  from  jet  to  wall  or  from  wall 
to  wall 

9.  Heat  transfer  between  the  jet  and  the  augmentation  air  is 
negligible  compared  with  enthalpy  changes  due  to  mixing 

10.  Radiation  from  the  jet  to  any  heat  exchanger  located  past  the 
augmenter  exit  is  negligble 

Jet  flow  through  the  augmenter  is  empirically  modelled  using  the 
curve  fits  of  Becker  et  al.  (Ref  2).  Their  data  are  from  cold  jets,  and 
their  test  rig  is  smaller  than  test  cells,  necessitating  some  modification 
to  the  empirical  constants.  Specifically,  Becker’s  relationships  have 
been  adjusted  to  be  able  to  predict  the  arrival  of  the  jet  to  the  wall 
at  locations  observed  experimentally  at  the  hush  houses,  NAS  Miramar 
(Ref  3  and  A). 

To  model  the  heat  transfer,  the  augmenter  tube  is  divided  into 
axial  segments,  with  the  jet,  augmentation  flow,  and  inner  and  outer 
liners  considered  as  separate  sections.  Temperature  is  assumed  constant 
throughout  each  section-segment.  Flows  are  perfectly  mixed.  Conserva¬ 
tion  of  energy  is  then  applied  to  each  segment  of  each  section,  and  the 
resulting  equations,  along  with  the  equations  of  continuity,  are  solved 
simultaneously  to  determine  temperatures  and  heat  fluxes  along  the 
augmenter  assembly.  Due  to  the  nonlinearities,  an  iterative  technique 
is  necessary;  the  Gauss-Seidel  method,  with  relaxation,  is  employed  for 
this  purpose. 

A  detailed  derivation  of  the  model,  including  a  discussion  of  the 
assumptions  and  the  accuracy  of  the  simulation,  is  enclosed  as  Appendix  A. 


ENERGY  RECOVERY  ECONOMICS 

By  supposition,  a  waterwall  or  a  convection,  water  tube  heat  exchanger 
is  added  to  the  augmenter.  Once  potential  test  cell  steam  (or  electricity) 
generation  has  been  determined,  its  accumulated  value  is  computed  and 
compared  with  the  initial  cost  of  the  waterwalls  or  water  tube  boiler, 
i.e.,  a  savings-to-investment  ratio  (SIR)  is  calculated. 

Appendix  B  includes  a  description  of  the  relationships  used  to 
calculate  SIR. 


A 


SAN  DIEGO  TEST  CELL  ENEKGY  RECOVERY 


Five  test  cells  located  in  the  San  Diego  area  are  examined  to 
determine  their  potential  for  energy  recovery.  These  cells  are  selected 
to  provide  a  variety  of  sizes  and  applications,  ranging  from  a  hush 
house  down  to  a  2-foot-diam  cell  used  for  testing  helicopter  engines. 
Specifications  and  operating  characteristics  are  summarized  in  Appendix  C. 

For  these  sites,  energy  is  recovered  in  the  form  of  350°F  saturated 
steam.*  The  energy  recovery  potential  of  each  cell  is  illustrated  hy 
means  of  two  figures.  First,  the  accumulated  present  value  of  the 
generated  steam  is  plotted  against  the  economic  life  of  the  steam  gen¬ 
eration  equipment.  A  Navy  average  is  used  for  the  value  of  the  steam. 
Predictions  of  discount  and  escalation  rates  are  employed  (Ref  5).  The 
second  approach  is  to  estimate  the  cost  of  the  energy  recovery  equipment 
and  plot  the  energy  recovery  potential  as  an  SIR  (again  against  the 
economic  life  of  the  equipment). 

Hush  House  No.  1,  NAS  Miramar 

Figures  2  and  3  show  the  potential  for  energy  recovery  in  the 
NAS  Miramar  Hush  House  No.  1.  The  augmenter  is  large  and  approximately 
elliptical,  with  an  overall  length  of  90  feet  and  a  cross-section  minor 
axis  of  about  11  feet.  A  water  tube  boiler  at  the  augm?nter  exit 
(Configuration  No.  1  of  Figure  1)  and  both  waterwall  configurations  are 
examined.  These  figures  reflect  testing  of  both  the  TF30  (F-14A)  and 
J79  (F-4)  at  intermediate  or  afterburning.  Although  other  engines  are 
tested,  the  TF30  and  J79,  together,  account  for  about  75%  of  the  total 
hush  house  usage. 

To  summarize:  the  generation  of  steam  in  this  facility  is  not 
economically  feasible.  Waterwalls  produce  very  little  steam  and,  in 
addition,  generate  this  small  amount  very  uneconomi ca 1 1 y .  A  large  water 
tube  heat  exchanger,  such  as  the  220,000  ft2  device  studied  here,  can 
generate  plenty  of  steam  during  afterburning  (A/B) ,  as  much  as  95,000  lb/hr, 
but  A/B  testing  is  limited,  and  these  large  boilers  are  very  expensive. 

This  particular  heat  exchanger's  estimated  cost  is  about  $6  million; 
thus,  the  SIR  is  low. 

"A”  Test  Cell,  NAS  Miramar 

The  NAS  Miramar  "A"  test  cell,  with  an  augmenter  cross-sectional 
diameter  of  about  6  feet,  can  be  classified  as  a  "mid-sized"  facility. 

The  economics  of  thermal  energy  recovery  from  this  test  cell  is 
summarized  on  Figures  4  and  5.  Waterwalls  located  near  the  augmenter 
inlet  are  not  economically  feasible.  Waterwalls  placed  along  the  last 
20  feet  of  the  augmenter,  however,  are  feasible  although  the  steam 
generation  is  low  (see  Figure  4).  For  this  configuration,  the  SIR 
reaches  one  after  about  22  years.  Water  tubes  placed  near  the  augmenter 


*Steam  generation  at  350°F,  138  psia  is  assumed  throughout  this  work. 
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exit  are  economically  feasible,  and  the  steam  generation  rate  is  high, 
averaging  480,000  lb/hr  at  A/B  and  about  40,000  lb/hr  at  intermediate 
power  settings.  The  SIR  reaches  one  after  7  years. 

Steam  loads  at  NAS  Miramar  average  about  30,000  lb/hr.  Therefore, 
480,000  Ib/hr  is  considerably  more  steam  than  the  base  can  consume. 

Thermal  energy  is  not  easily  stored;  a  customer  for  the  excess  steam  is 
necessary.  If  not  found,  much  of  the  steam  generated  during  A/B  testing 
would  not  be  used,  and  the  above  water  tube  SIR  is  inflated. 

The  TF30  is  the  only  engine  tested,  yet  this  cell  has  had  a  very 
high  frequency  of  A/B  testing  -  62.6  total  hours  in  1982  -  and  this 
accounts  for  much  of  the  positive  energy  recovery  economics.  The 
practicality  of  installing  some  type  of  energy  recovery  system  depends, 
to  a  very  large  extent,  on  continuing  this  schedule  into  the  fu*  ° . 

Test  Cell  11,  NARF  North  Island 

Test  Cell  11  is  a  small  cell  used  to  test  the  T58  helicop  engine. 
The  diameter  of  the  augraenter  is  about  2-1/2  feet;  its  length  ' "  feet. 

Figures  6  and  7  summarize  the  energy  recovery  capability  i  »s 
cell.  Again,  the  installation  of  waterwalls  is  not  economically  feasible. 
The  use  of  water  tubes  placed  at  the  exit  to  the  augmenter  is  economic¬ 
ally  feasible,  with  the  SIR  reaching  one  after  10  years.  Approximately 
8,000  lb/hr  of  steam  is  generated  in  a  6,600  ft2  water  tube  heat  exchanger, 
only  a  dribble  compared  with  the  steam  load  at  NAS  North  Island. 

Test  Cell  14,  NARF  North  Island 

Test  Cell  14  is  a  mid-sized  facility  used,  primarily,  to  test  the 
J79  gas  turbine  engine.  The  specifications  of  this  cell  are  similar  to 
those  of  the  NAS  Miramar  "A"  test  cell. 

The  economics  of  energy  recovery  from  cell  14  are  illustrated  by 
Figures  8  and  9.  This  cell  is  used  i nfrequent ly,  and  any  form  of  energy 
recovery  is  economically  unfeasible. 

Test  Cell  19,  NARF  North  Island 

Test  Cell  19  is  also  a  mid-sized  cell  used  to  test  the  J79.  It  has 
an  augmenter  cross-sectional  diameter  of  about  8  feet  with  a  length  of 
almost  80  feet. 

Figures  10  and  11  show  the  potential  for  energy  recovery  in  Test 
Cell  19.  Both  waterwalls  and  water  tubes  are  economically  marginal 
(i.e.,  within  possible  error  in  estimating  the  cost  of  the  equipment). 

The  steam  generation  rate  is  moderate.  A  116,000  ft2  water  tube  boiler 
will  generate  an  average  of  20,000  Ib/hr  of  steam  during  intermediate 
power  tests  and  about  310,000  lb/hr  during  A/B. 

A  comparison  with  the  NAS  Miramar  "A"  test  cell  is  informative. 
Although  Test  Cell  19  is  used  more  frequently,  210*  hr/yr  compared  with 
about  135  hr/yr  for  the  Miramar  test  cell,  much  less  energy  recovery 


intermediate  and  afterburning  only.  Idle  and  part  power  settings  do 
not  produce  temperature  gradients  high  enough  for  thermal  energy 
recovery  and  are  not  considered. 
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potential  exists.  The  J79  exhausts  somewhat  less  than  the  TF30,  hut  the 
major  difference  is  the  low  A/B  testing,  only  18  hr/yr  in  Test  Cell  19 
compared  with  more  than  62  hr/yr  in  the  Miramar  "A"  test  cell. 


PARAMETRIC  STUDY  OF  TEST  CELL  ENERGY  RECOVERY 

A  perusal  of  the  five  cells  just  analyzed  plus  the  others  described 
in  Appendix  C  suggests  that  the  test  cells  can  be  segregated  into  three 
size  categories  for  parametric  examination:  large,  mid-sized,  and 
small.  The  distinction  is  based  upon  augmenter  diameter.  Large  test 
cells  will  be  represented  by  NAS  Miramar  Hush  House  No.  1,  mid-sized 
cells  by  the  NAS  Miramar  "A"  test  cell,  and  small  cells  by  NARF  North 
1 s land  Test  Cell  11. 

Again,  both  the  accumu 1 ated  present  value  of  the  steam  and  SIR  will 
be  used  as  measures  of  feasibility.  Perhaps  somewhat  redundant,  the 
value  of  the  steam  is  an  indication  of  steam  generation  magnitudes  while 
SIR  is  an  indication  of  economics.  Both  figures  are  included  also  for 
documental  ion  purposes.  They  can  be  used  to  evaluate  test  cells  not 
specifically  covered. 

Large  Test  Cells 

It  was  not  emphasized  in  the  discussion  of  NAS  Miramar  Hush  House 
No.  1,  but  energy  can  be  recovered  in  the  large  test  cells,  by  water- 
walls  or  by  downstream  water  tubes,  only  durit  ■>  afterburning  tests.  All 
other  engine  power  settings  are  nonproductive.  This  is  illustrated  by 
Figures  12  and  13,  which  isolate  the  energy  recovery  capabilities  of  the 
different  power  settings.  Augmenter  exit  water  tubes  are  employed  in 
these  figures.  This  is  the  most  efficient  energy  recovery  configuration 
for  Hush  House  No.  1.  Note  that  some  of  the  hush  houses  are  capable  of 
simul taneously  testing  two  engines. 

Although  the  A/B  steam  generation  rate  is  potentially  high,  perhaps 
100,000  lb/hr  while  testing  the  TF30,  the  total  hours  of  A/B  testing  is 
very  low.  A  testing  frequency  of  180  hr/yr,  or  an  average  of  0.5  hr/day, 
was  assumed  for  Figures  12  and  13.  This  schedule  is  characteristic  of 
intermediate  power  tests  in  the  San  Diego  cells,  but  afterburner  tests 
are  conducted  less  frequently  and  are  of  a  much  shorter  duration.  An 
A/B  test  frequency  of  30-40  hr/yr  is  a  better  estimate,  thus,  a  25-year 
SIR  of  0.1,  rather  than  0.5  (see  Figure  13). 

The  absence  of  any  potential  for  energy  recovery  during  gas  turbine 
intermediate  power  tests  can  be  attributed  to  the  high  augmentation  air 
flow  rates  inherent  in  the  large  test  cells.  Augmentation  air  flows  are 
typically  1,200  to  1,500  lb/sec  compared  with  400  to  500  lb/sec  through 
mid-sized  cells.  The  engine  exhaust  is  cooled  down  to  about  200°F  by 
the  time  it  reaches  the  wall;  it  is  even  cooler  when  it  reaches  the 
augmenter  exit.  At  this  temperature,  the  exhaust  cannot,  of  course,  be 
used  to  generate  350°F  steam. 

Energy  recovery  will  never  be  economically  feasible  in  the  large 
test  cells  unless  a  system  capable  of  recovering  energy  during  inter¬ 
mediate  power  testing  is  employed.  Either  of  two  methods  can  accomplish 
this : 
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1.  Recover  the  energy  from  the  gas  turbine  exhaust  before  it  is 
mixed  with  the  augmentation  air  and  cooled 

2.  Transfer  the  energy  from  the  200°F  exhaust/air  mixture  to  some 
highly  volatile  working  fluid  and  use  this  fluid  to  generate 
elect  r i c i t y 

The  first  method  requires  placing  the  heat  exchanger  up  near  the 
nozzle  of  the  engine.  Configuration  No.  4  of  Figure  1.  Figures  14 
and  15  parametrically  examine  this  configuration  when  employed  on  large 
test  cells.  Figure  14  shows  the  potential  of  a  500  ft2  heat  exchanger. 
Approximately  85  2-inch-diam  water  tubes  stretch  across  the  augmenter 
just  behind  the  test  engine,  probably  near  the  maximum  number  of  tubes 
that  could  reasonably  be  placed  there.  Note  that  it  takes  an  average  of 
4  hr/day  of  testing  to  make  this  heat  exchanger  feasible. 

Figure  15  sets  the  daily  test  cell  operation  at  0.5  hr/day  and 
illustrates  the  effect  of  increasing  heat  transfer  area  (i.e.,  increas¬ 
ing  the  number  of  water  tubes).  It  takes  about  4,000  ft2  of  surface 
area  to  be  economically  feasible.*  This  is  equivalent  to  700  tubes 
stretched  across  the  augmenter.  In  addition  to  the  design  problems 
innate  to  this  configuration,  the  pressure  drop  across  700  tubes  would 
probably  induce  unacceptable  engine  back  pressures.  Configuration  No.  4 
is  not  workable.** 

The  second  method  of  recovering  energy,  employing  a  working  fluid 
with  a  low  boiling  point  to  generate  electricity,  is  the  topic  of  a 
subsequent  section. 

Mid-Sized  Test  Cells 


The  examination  of  the  San  Diego  cells,  particularly  the  NAS  Miramar 
"A"  test  cell,  showed  that  energy  recovery  in  mid-sized  test  cells  could 
be  economically  feasible  with  a  favorable  test  schedule.  Steair  can  be 
generated  from  the  exhaust  of  gas  turbines  tested  at  intermediate  power 
in  mid-sized  cells.  Augmentation  air  flow  rates  are  moderately  low,  and 
the  exhaust  temperature  is  still  about  400°F  to  450°F  when  it  reaches 
the  augmenter  wall. 

Both  water  tubes  and  downstream  waterwalls  can  generate  steam;  but 
the  upstream  waterwalls  (Configuration  No.  2)  are  of  no  use.  Radiation 
from  jet  to  the  upstream  walls  is  maximum;  however,  convection  heat 
transfer  from  these  walls  to  the  cold  augmentation  air  is  also  maximum 
(see  Figures  A-4,  A-5,  and  A-6).  There  is  little  net  heat  transferred. 

Figures  16  through  19  parametrically  show  the  potential  for  energy 
recovery  from  engines  being  tested  at  intermediate  power  in  the  mid-sized 
cells.  A  comparison  of  Figures  16  and  17  with  Figures  18  and  19  discloses 
the  superiority  of  the  water  tubes  over  the  waterwalls.  Figures  16 


*Probably  more  since  the  cost  of  the  heat  exchanger  was  not  increased 
with  heat  transfer  area. 

**The  mid-sized  and  small  test  cells  already  have  the  potential  to 
generate  steam  at  intermediate.  For  them,  there  is  no  need  to 
consider  such  a  radical  configuration. 


8 


and  18  show  the  water  tubes  generate  about  60  times  more  steam  (actually, 
40,000  lb/hr  compared  with  only  700  lb/hr  by  the  waterwalls).  Figures  17 
and  19  show  the  water  tube  cr  riguration  has  an  SIR  about  four  times 
higher.  The  superiority  of  water  tubes  exists  with  all  sizes  of  test 
cells  and  can  be  attributed  to  the  smaller  heat  transfer  surface  area  of 
the  waterwa 11s. 

Yet,  for  even  the  downstream  water  tube  heat  exchanger  to  be  econom¬ 
ically  feasible  for  mid-sized  test  cell  energy  recovery,  about  1  hr/day 
of  intermediate  testing  is  required.  This  frequency,  365  hr/yr,  is 
about  twice  the  usage  averaged  by  the  San  Diego  facilities.  More  than  a 
trace  of  A/B  testing  is  still  required  in  order  to  make  energy  recovery 
feasible . 

A  huge  amount  of  steam  can  be  generated  when  gas  turbines  are 
tested  at  A/B  power  settings  in  a  mid-sized  cell.  Figure  20  is  analogous 
to  Figure  17  except  that  the  engines  are  being  tested  while  afterburning 
rather  than  at  intermediate.  Examining  these  two  figures  together  is 
enlightening.  They  show  the  combinations  of  intermediate  and  A/B  testing 
that,  roughly,  make  energy  recovery  feasible  in  mid-sized  cells,  for 
example,  0.5  hr/day  (180  hr/yr)  of  intermediate  power  testing  with  an 
SIR  of  about  0.3  after  25  years  combined  with  0.05  hr/day  (18  hr/yr)  of 
A/B  testing  with  an  SIR  of  about  0.7. 

Small  Test  Cells 


Energy  recovery  in  small  test  cells  involves  different  limitations 
than  energy  recovery  in  the  larger  cells.  There  is  little  or  no  augmen¬ 
tation  air;  therefore,  temperature  gradients  are  higher.  However, 
engines  tested  in  these  cells  are  much  smaller;  there  is  much  less 
thermal  energy  available  for  steam  generation.  Also  there  is  no  after¬ 
burning.  It  follows  that  steam  can  usually  be  generated  in  these  cells 
but  not  very  much  of  it. 

Figures  21  and  22  summarize  the  energy  recovery  potential  of  small 
test  cells.  Steam  is  being  generated  at  a  rate  of  about  8,000  lb/hr. 

The  25-year  break-even  SIR  occurs  when  test  cell  usage  reaches  about 
36  hr/yr.  A  steam  generation  rate  of  only  8,000  lb/hr  for  perhaps  1/2 
hr/day  would  not  contribute  much  toward  satisfying  the  steam  load  on  the 
base  and,  despite  any  favorable  SIR,  would  not  warrant  the  tribulations 
involved  in  maintaining  the  boiler. 

Thermal  energy  recovery  in  small  test  cells  might  become  more 
attractive  if  several  of  them  could  be  combined.  The  steam  generation 
rate  would  not  improve,  unless  several  cells  could  be  scheduled  for 
testing  at  the  same  time,  but  the  steam  output  would  be  steadier,  and 
the  total  steam  contribution  much  higher  (for  example,  cells  9,  10,  11, 
and  12  at  NARF  North  Island  are  all  within  a  space  of  60  feet  and  could, 
conceivably,  all  be  ducted  to  a  single  heat  exchanger). 


ENERGY  RECOVERED  AS  ELECTRICITY 

The  energy  recovery  analyses  presented  to  this  point  were  conducted 
with  steam  as  the  product.  Electricity  could  also  have  been  produced. 
Local  conditions  might  make  one  or  the  other  more  valuable  at  a  particu¬ 
lar  site,  but,  by  using  Navy  averages  for  the  values  of  energy,  supply 
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and  demand  can  be  expected  to  prevail.  Similar  results  would  be  antici¬ 
pated  from  similar  analyses  using  either  steam  or  electricity  as  the 
product  of  energy  recovery. 

Two  situations  have  been  encountered,  however,  in  which  the  genera¬ 
tion  of  electricity  would  seem  to  be  preferable.  The  first  occurs  when 
trying  to  recover  energy  from  a  large  cell  testing  a  gas  turbine  at 
intermediate  power.  Engine  exhaust  temperatures,  diluted  by  the  high 
augmentation  air  flow,  are  too  low  to  generate  350°F  steam.  Since 
energy  recovery  during  intermediate  power  testing  is  essential  if  the 
concept  is  to  be  economically  feasible,  some  other  product  must  be 
found.  It  was  suggested  that  the  energy  in  the  jet  exhaust  could  be 
transferred  to  some  highly  volatile  working  fluid  (at  high  pressures  but 
low  temperatures)  and  this  fluid  expanded  through  a  turbine  to  generate 
electrici ty . 

The  second  situation  in  which  the  generation  of  electricity  would 
appear  to  be  preferable  occurs  during  A/B  testing  in  mid-sized  cells 
located  on  the  smaller  Navy  bases.  The  problem  is  excess  steam.  If  a 
large  water  tube  boiler  is  employed,  as  much  as  500,000  lb/hr  of  350°F 
steam  can  be  generated  from  the  exhaust  of  a  TF30  during  afterburning 
tests.  This  is  more  steam  than  many  facilities  can  consume.  Yet,  for 
energy  recovery  to  be  feasible  in  mid-sized  cells,  the  energy  of  after¬ 
burning  engines  must  tie  utilized.  Either  the  steam  must  be  stored  or  a 
customer  found  for  the  excess.  Neither  is  practical.  Thermal  energy  is 
not  easily  stored,  and  potential  customers  for  the  steam  would  probably 
not  be  interested  in  such  an  erratic,  unpredictable  source.  The  genera¬ 
tion  of  electricity  rather  than  steam  may  be  a  solution  to  the  problem. 
Many  utility  companies  will  buy  (or  compensate  for)  electricity  gener¬ 
ated  locally. 

Thus,  the  feasibility  of  energy  recovery  in  the  form  of  electricity 
is  examined  for  the  large  and  mid-sized  test  cells.  These  analyses  are 
conducted  by  again  employing  the  numerical  model.  It  is  assumed  that 
electricity  will  be  generated  off  a  Rankine  cycle,  with  energy  supplied 
from  the  jet  exhaust  through  waterwalls  or  through  a  water  tube  heat 
exchanger.  A  hypothetical  organic  fluid  is  substituted  for  the  steam-; 
otherwise,  the  analyses  are  identical  to  previous  efforts.  The  heat 
transfer  efficiency  is  inherent  in  the  model.  Efficiencies  of  the  other 
Rankine  cycle  components  are  stipulated. 


n  = 

pump 

0.75 

^turbine 

0.85 

^generator 

0.90 

and  multiplied  by  the  energy  recovered  from  the  exhaust  to  provide  a 
magnitude  for  electricity  generation. 

*The  organic  fluid  is  assumed  to  vaporize  at  120°F  with  a  heat  of 
vaporization  of  145  Btu/lb. 


Figures  23  and  24  show  the  potential  For  energy  recovery,  as  elec¬ 
tricity,  from  jet  engines  being  tested  at  intermediate  power  in  the 
large  cells.  An  improvement  over  steam  is  shown  as  would  be  expected 
(energy  recovery  in  the  form  of  350°F  steam  was  negligible).  The 
improvement  is  small,  however.  A  25-year  break-even  requires  a  4  hr/day 
intermediate  power  test  frequency,  about  eight  times  the  average  of  the 
San  Diego  cells.  Too  much  A/B  testing  is  still  necessary  in  order  to 
make  energy  recovery  economically  feasible  in  the  large  cells. 

The  generation  of  electricity  in  mid-sized  cells  shows  somewhat 
more  promise.  When  testing  at  intermediate,  there  is  little  difference 
in  the  economics  of  energy  recovery  with  electricity  or  with  steam  as 
the  product.  Figures  16  and  17  can  be  used  to  evaluate  energy  recovery 
within  reasonable  accuracy.  With  the  working  fluid  assumed*  for  this 
study,  the  generation  of  electricity  from  the  energy  in  the  exhaust  of 
afterburning  engines  is  less  favorable.  Figures  25  and  26  show  the 
potential  for  generating  electricity  in  mid-sized  cells  testing  engines 
at  A/B.  The  SIR  for  electricity  generation  is  roughly  one-third  the  SIR 
calculated  for  an  analogous  system  generating  steam.  Still,  the  energy 
recovered  as  electricity  from  afterburning  engines  is  equivalent  to  the 
energy  required  to  generate  160,000  lb/hr  of  steam.  This  is  quite 
significant  if  all  this  energy  can  now  be  used. 

Economic  feasibility  is  determined  by  energy  recovery  during  both 
intermediate  and  afterburning  tests.  An  example  was  provided  earlier  to 
illustrate  the  feasibility  of  energy  recovery  as  steam  in  mid-sized 
cells.  It  is  informative  to  repeat  this  same  example  with  electricity 
as  the  product.  Generating  steam,  0.5  hr/day  of  intermediate  testing 
combined  with  0.05  hr/day  of  afterburning  tests  yields  an  energy  recovery 
SIR  of  one  after  25  years.  Generating  electricity  (from  Figures  17 
and  26),  0.5  hr/day  of  intermediate  testing  must  be  combined  with 
0.15  hr/day  (54  hr/yr)  of  A/B  testing  to  yield  the  same  SIR. 

The  generation  of  electricity  in  the  NAS  Miramar  "A"  test  cell,  a 
cell  with  a  high  frequency  of  afterburning  tests,  would  be  expected  to 
be  feasible,  and  it  is.  This  cell  is  used  an  average  0.2  hr/day 
(72.9  hr/yr)  for  intermediate  power  tests,  yielding  an  SIR  of  0.25  after 
25  years;  and  used  approximately  0.17  hr/day  (62.6  hr/yr)  for  A/B  tests, 
yielding  an  SIR  of  0.85.  Thus,  the  total  savings-to-investment  ratio 
reaches  about  1.1  after  25  years. 


CONCLUSIONS 

Thermal  energy  recovery  is  not  economically  feasible  in  hush  houses 
and  large  test  cells.  The  augmentation  air  flow  rates  are  very  high, 
rapidly  cooling  the  engine  exhaust  and  eliminating  temperature  gradients 
necessary  for  steam  generation.  Electricity  generation  off  an  organic 
Rankine  cycle,  although  economically  better,  still  requires  too  much 
test  cell  usage  to  be  feasible. 


*The  selection  of  an  optimum  organic  fluid  is  itself  a  major  task, 
the  choice  varying  with  cell  geometry,  test  schedule,  and  other 
conditions . 
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Energy  recovery  in  most  small  test  cells  is  unfeasible.  These 
facilities  are  normally  used  for  testing  a  single  class  of  gas  turbine; 
the  hours  of  operation  are  limited.  Even  if  the  SIR  is  favorable,  the 
engines  are  small  and  steam  or  electricity  generation  is  so  low  that  it 
would  not  be  worth  the  effort.  An  exception  might  occur  if  energy 
recovery  systems  for  several  small  cells  could  be  combined  (e.g.,  ducting 
the  augmenter  exit  gases  of  several  cells  to  a  single  heat  exchanger). 

Energy  recovery  in  the  mid-sized  cells  is  marginal,  with  the  feasi¬ 
bility  depending  upon  cell  usage  and  energy  loads.  Steam  generation 
rates  can  be  excellent  because  augmentation  flows  are  lower  than  in  the 
large  cells,  yet  mid-sized  cells  are  often  used  for  testing  the  larger 
engines.  If  the  Naval  facility  is  small,  steam  generation  rates  may 
even  exceed  the  steam  demand,  particularly  when  testing  afterburning. 

When  this  occurs,  the  generation  of  electricity  is  an  alternative  to  be 
considered. 

A  water  tube  boiler  located  outside  the  augmenter  exit,  Configuration 
No.  1  of  Figure  1,  is  the  only  practical  choice  for  an  energy  recovery 
heat  exchanger.  Augmenter  waterwalls  and  water  tubes  inside  the  augmenter 
have  too  small  a  heat  transfer  surface  area.  The  steam  generation 
potential  is  much  too  low. 


COMMENTS 

Boiler  design  problems  have  not  been  examined,  but  they  will  tend 
to  discourage  thermal  energy  recovery.  Boilers  do  not  function  particu¬ 
larly  well  when  subjected  to  cyclic  operation.  Not  only  is  the  incessant 
mechanical  expansion  and  contraction  of  the  pressure  components  a  diffi¬ 
cult  design  problem,  but  the  water  swell  occurring  at  each  start  up 
presents  a  control  problem  and  certainly  results  in  a  significant  quantity 
of  boiler  water  wastage.  Tube  or  waterwall  chemical  corrosion  can  be 
expected.  Surface  erosion  is  also  probable.  A  gas  velocity  of  150  ft/sec 
is  about  the  practical  limit  that  the  heat  exchanger  components  could 
withstand.  If  the  heat  exchanger  is  outside  the  augmenter,  the  engine 
exhaust  can  be  slowed  to  acceptable  velocities  before  entering  the 
device.  For  waterwalls  and  for  water  Lubes  inside  the  augmenter,  heat 
transfer  surface  erosion  would  be  a  very  serious  problem.  In  addition, 
water  tubes  inside  the  augmenter  would  have  to  withstand  blasts  of 
3,000°F  exhaust  whenever  afterburning  tests  were  conducted. 

Although  not  intentionally  included  in  these  analyses,  some  knowledge 
of  potential  side  effects  was  acquired  as  a  byproduct  of  the  heat  transfer 
calculations.  The  only  energy  recovery  configuration  that  has  a  signifi¬ 
cant  effect  on  flow  characteristics  inside  the  augmenter  is  upstream 
water  tubes  (Configuration  No.  U  of  Figure  1).  It  is  unlikely  that  this 
configuration  will  ever  be  employed;  the  design  problems  are  presumably 
insurmountable.  Very  little  energy  is  removed  from  the  jet  by  water- 
walls,  and  side  effects  will  probably  be  negligible  (compare  stack  gas 
temperatures  on  Table  1  with  and  without  waterwalls).  The  downstream 
water  tubes,  Configuration  No.  1,  have  an  appreciable  effect  on  flow 
characteristics,  but  only  outside  the  augmenter  tube.*  Thus,  side 
effects  to  test  cell  energy  recovery  may  not  be  a  major  factor. 


*The  effect  on  engine  back  pressures  could  still  be  a  problem. 
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Stark  gas  temperatures  are  another  byproduct  of  the  heat  transfer 
calculations.  Table  1  summarizes  the  effect  of  energy  recovery  boilers 
on  the  temperature  of  the  jet  exhaust.  As  suggested  above,  the  effert 
of  the  waterwalls  is  very  small.  Convection,  water  tube  configurations, 
however,  signi f icant 1 y  lower  stack  gas  temperatures,  a  beneficial  side 
effect  when  further  environmental  or  acoustical  treatments  are  necessary. 
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NOMENCLATURE 


qCOND 

qCONV 

qRAD 

qSTM 

r* 


R 

R(Z) 


Surface  area 
Specific  heat 

Craya-Curtet  Number,  Equation  A-l 
Assumed  value  of  steam.  Equation  B- 1 
Diameter 

Equivalent  diameter 
Discount  rate 
Boiler  effectiveness 

Escalation  rate  of  the  cost  of  steam  or  electricity,  Equation  B-2 

Radiation  shape  factors 

Enthalpy 

Convective  heat  transfer  film  coefficient 
Coefficient  of  thermal  conductivity 
Mass  flow  rate 

Augmentation  flow  entrained  by  the  jet 

Initial  cost  of  waterwalls  or  heat  exchanger.  Equation  B-3 

Present  value 

Conduction  heat  transfer 

Convection  heat  transfer 

Radiation  heat  transfer 

Heat  transferred  from  jet  and  augmentation  airflow  to  the  steam 

Ratio  of  jet  nozzle  area  to  cross-sectional  area  of  augmenter, 
Equation  A-l 

Radius 

Radius  at  distance  Z  from  the  augmenter  inlet 
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r 


SIR  Savings-to-investment  ratio,  Equation  B-3 

T  Temperature 

TIME  Average  utilization  of  the  test  cell,  Equation  B- 1 

UA  Product  of  boiler  overall  heat  transfer  coefficient  and  heat 

transfer  area 

VALUE  Value  of  the  steam  or  electrictiy  generated 
V  Velocity 

Z  Distance  from  the  augmenter  inlet 

a  Radiation  absorptivity 

Ahf^  Heat  of  vaporization  of  steam 

AT^  Logarithmic  mean  temperature  difference.  Equation  A-10 

e  Radiation  emissivity 

M  Viscosity 

p  Density 

o  Stefan-Boltzmann  constant 

r|  Efficiency 


Subscri  £tS 

AUG  Refers  to  augmenter  tube 

GAS  Refers  to  the  mixed  jet  and  augmentation  flow 

I,  J  Segment  node  designations,  see  Figure  A-i 

IL  Refers  to  inner  liner  of  augmenter 

I-»J  From  segment  "I"  to  segment  "J",  etc. 

JET  Refers  to  jet  exhaust 

JET-»IL  From  the  jet  to  the  inner  liner  of  augmenter,  etc. 

N  Number  of  axial  segments  of  the  augmenter  tube,  see  Figure  A-l 

NN  Economic  life  of  waterwalls  or  energy  recovery  heat  exchanger, 

Equation  B-2 


OL 


Refers  to  outer  liner  of  augmenter 


REF  Reference,  used  for  defined  properties 

SEC  Refers  to  augmentation  (secondary)  air  flow 

STACK  Refers  to  jet  and  augmentation  flow  leaving  the  convection  boiler 

STM  Refers  to  steam  generated 

TUBE  Refers  to  water  tubes  of  convection  heat  recovery  heat  exchanger 

»  Refers  to  ambient  conditions 

Abbreviations 

A/B  Afterburning 

INTER  Intermediate 

ORC  Organic  Rankine  cycle 
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Table  ].  Effect  of  Gas  Turbinp  Test  Cell  Energy  Recovery  on  Stack  Gas  Temperatures 


Engine 

Type  of  Energy 
Recovery  Equipment 

Total  Gas  Flow 
Out  the  Stack, 

Jet  Exhaust  and 
Cooling  Flow 
(lb/sec) 

Avg  Stack  Gas 
Temp  Without 
Energy  Recovery 
(°F) 

Avg  Stack  Gas 
Temp  With 
Energy  Recovery 
(°F) 

|  a.  Large  Test  Cel] 

J79  afterburning 

(1)  With  water  tubes  at 
augmenter  exit, 
250,000  ft2 

1 ,700 

420 

370 

(2)  With  waterwalls 

along  last  50  feet 
of  augmenter 

1 ,700 

420 

415 

TF30  afterburning 

(1)  With  water  tubes  at 
augmenter  exit, 
220,000  ft2 

1  ,450 

650 

425 

(2)  With  waterwalls 
along  last  60  ft 
of  augmenter 

]  ,450 

650 

640 

b.  Mid-5 

>ized  Test  Cell 

TF30  at  intermediate 

CD  Water  tubes  at  aug¬ 
menter  exit, 

80,000  ft2 

660 

425 

370 

(2)  Waterwalls  along 
last  20  ft  of  aug¬ 
menter 

660 

425 

420 

TF30  afterburning 

(1)  Water  tubes  at  aug¬ 
menter  exit, 

80,000  ft2 

660 

1  ,220 

560 

(2)  Waterwalls  along 
last  20  ft  of  aug¬ 
menter 

660 

1  ,220 

1  ,205 

c .  Sma 1 1  Test  Cel  1 

T58  at  intermediate 

(l)  Water  tubes  at  aug- 
mente r  exit, 

6,600  ft2 

20 

820 

465 

(2)  Waterwalls  along 
las*.  38  ft  of  aug¬ 
menter 

20 

820 

810 
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Figure  2.  Generation  of  steam  in  gas  turbine  engine  test  facilities:  NAS  Miramar  Hush  House  No. 


1)  TF30  AND  J79  AT  INTER  OR  A/B 

2)  TEST  CELL  UTILIZATION: 

TF30SINTER  150  HRS/YR  TF30(?A/B  15  HRS/YR 
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Figure  3.  Economics  of  gas  turbine  engine  test  facility  energy  recovery:  NAS  Miramar  Hush  House  No. 
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Figure  4.  Generation  of  steam  in  gas  turbine  engine  test  facilities:  NAS  Miramar  Test  Cell 
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Economics  of  gas  turbine  engine  test  facility  energy  recovery:  NAS  Miramar  Test  Cell  " 


Figure  6.  Generation  of  steam  in  gas  turbine  engine  test  facilities:  NARF/NAS  North  Island  Test  Cell  11. 
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'  1 )  J 79  AT  INTERMEDIATE 

.2)  TEST  CELL  UTILIZATION-  15.7  HRS/YR 

!3)  COST  OF  HEAT  RECOVERY  EQUIPMENT: 

WATER  TUBE  BOILER  @  $1,700,000. 
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Figure  9.  Economics  of  gas  turbine  engine  test  facility  energy  recovery:  NARF/NAS  North  Island  Test  Cell  14. 


NOTES: 

l?ooooo -|  (i)  J79  flT  INTER  OR  fl/B 

(2)  TEST  CELL  UTILIZATION: 196.3  HRS/YR  INTER 

17.7  HRS/YR  A/B 


0X0 


Economic  Life,  Years 

Figure  10.  Generation  of  steam  in  gas  turbine  engine  test  facilities.  NARF/NAS  North  Island  Test  Cell  19. 
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Accumulated  Present  Value  of  Steam  Generated,  $ 

Figure  11.  Economics  of  gas  turbme  engme  test  facility  energy  recovery.  NARF/NAS  North  island  Test  Cell  19. 


f  1 )  WATER  TUBE  HEAT  EXCHANGER  AT  AUGMENTER  EXIT 
12)  HEAT  TRANSFER  AREA  -  250000  SQ  FT 


Economic  Life,  Years 

Figure  12.  Effect  of  engine  power  setting  on  the  generation  of  steam  in  a  typical 
large  gas  turbine  test  cell. 
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Figure  13.  Effect  of  engine  power  setting  on  the  economics  of  energy  recovery  in  a 
typical  large  gas  turbine  test  cell. 
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heat  transfer  area  on  energy  recovery  in  a  typical  large  test  cell. 
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Economic  Life,  Years 

Figure  21.  Steam  generation  in  a  typical  small  gas  turbine  test  cell  using  water  tubes 
at  the  augmcnter  exit. 
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Economic  Life,  Years 

Figure  25.  Generation  of  electricity  from  thermal  energy  recovered  in  a  typical  mid-sized 
gas  turbine  test  cell  using  water  tubes  at  the  augmenter  exit. 
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Economic  Life,  Years 

Figure  26.  Economics  of  electricity  generation  from  thermal  energy  recovered  in  a  typical 
mid-sized  gas  turbine  test  cell  using  water  tubes  at  the  augmenter  exit. 


Appendix  A 


SIMULATION  OF  A  JET  ENGINE  TEST  CELL 
AUGMENTER  TUBE 


The  problems  inherent  in  simulating  the  operation  of  test  cell 
augmenter  tubes  are  extremely  complex.  Re-examine  Figure  1  and  consider 
the  following  augmenter  characteristics:  multiple,  nonconcentric  flows 
at  different  velocities;  turbulence;  flow  separation;  all  three  modes  of 
heat  transfer,  simultaneously;  extreme  changes  in  temperature  and,  thus, 
variable  properties;  and  nonsymmetric  three-dimensional  geometries.  In 
addition,  some  tests  are  of  short  duration,  resulting  in  a  highly  tran¬ 
sient  situation;  e.g.,  testing  with  afterburners  might  involve  30-second 
cycles  with  3,000°F  exhaust.  These  complexities  lead  to  a  transient, 
nonlinear  problem  with  coupled  energy,  mass,  and  momentum  transfer. 

To  make  an  augmenter  simulation  attainable,  assumptions  are  necessary. 
For  this  present  work,  the  emphasis  is  on  heat  transfer.  Can  enough 
steam  be  generated,  via  energy  transferred  from  the  jet  exhaust,  to 
warrant  the  purchase  of  a  boiler?  The  modelling  tactics  will  be  to 
simulate  the  fluid  flow  (mass  and  momentum)  empirically  but  use  first 
principles  to  determine  heat  transfer.  For  this  plan-of-attack ,  the 
following  assumptions  are  made: 

1.  Flows  are  one-dimensional,  passing  through  a  cylindrical 
augmenter 

2.  Steady  state,  steady  flow  is  stipulated 

3.  Both  the  jet  and  augmentation  air  are  ideal  gases 


A- 1 


4.  Flow  within  the  jet  and  within  the  augmentation  air  can  be 
considered  turbulent  and  perfectly  mixed 

5.  Pressure  gradients  along  the  augmenter  are  negligible  compared 
with  temperature  gradients  and  can  be  neglected 

6.  The  jet  nozzle  is  located  exactly  at  the  inlet  to  the  augmenter 
tube 

7.  The  augmenter  walls  are  black,  the  jet  is  gray,  and  the  aug¬ 
mentation  air  is  transparent  to  thermal  radiation 

8.  Axial  radiation  along  the  jet  has  a  negligible  effect  on 
temperatures;  all  radiation  is  from  jet  to  wall  or  from  wall 
to  wall 

9.  Heat  transfer  between  the  jet  and  the  augmentation  air  is 
negligible  compared  with  enthalpy  changes  due  to  mixing 

10.  Radiation  from  the  jet  to  any  heat  exchanger  located  past  the 
augmenter  exit  is  negligble 

The  validity  of  several  of  the  assumptions  and  the  accuracy  of  the 
resulting  simulation  are  discussed  in  a  subsequent  section. 


SIMULATING  THE  JET  EXHAUST 

The  jet  is  empirically  modelled  using  the  curve  fits  of  Becker 
et  al.  (Ref  2).  Their  data  are  from  cold  jets,  and  their  test  rig  is 
smaller  than  test  cells,  necessitating  some  modification  to  the  empirical 
constants.  Specifically,  these  relationships  have  been  adjusted  to  be 
able  to  predict  the  arrival  of  the  jet  to  the  wall  at  locations  observed 
experimentally  at  the  hush  houses,  NAS  Miramar  (Ref  3  and  4). 
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Following  Becker's  technique,  jet  flow  geometry  is  extrapolated 
from  the  NAS  Miramar  data  by  means  of  a  nondimensional  similarity  param¬ 
eter  usually  referred  to  as  the  Craya-Curtet  number,  C  ,  a  ratio  of  the 
kinematic  to  dynamic  mean  inlet  velocities, 


r*2vJET«»  *  ('"-‘Hec 


r*2vJET(0)  •  (0.5-r*2)V2EC  -  0.5 


r"  VJET(0^  +  (l'r"  ,VSEC 
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where:  r*  =  ratio  of  jet  nozzle  area  to  cross-sectional  area  of 

the  augmenter  tube 

Vjet^O)  =  velocity  of  jet  at  the  augmenter  inlet 

VgEc  =  velocity  of  augmentation  (secondary)  airflow 
Then,  for  0  <  5  0.7, 


Rjet(Z)  =  Rjft(°)  +  0.131  Z 


where:  Z  =  4.07  exp  (3.54  C  ) 


1  - 


Z/(ZR  W 
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JET 

R 


(Z)  =  radius  of  jet  at  a  distance  Z  from  the  augmenter  inlet 


AUG 


=  radius  of  augmenter  tube 


and  for  0.7  <  , 


RJET(2)  *  Rjet(0)  *  0.126  Z 


(A-2 ') 


The  geometry  of  the  jet  is  fully  defined  by  Equations  A-l  and  A-2. 
Although  it  appears  that  fluid  properties  are  being  neglected,  it  should 
be  noted  that  the  augmentation  airflow  rate  is  input  directly,  rather 
than  calculated.  This  partially  introduces  properties  of  both  the  jet 
and  augmentation  air  through  the  variable  V<jE£- 

A  final  flow  relationship  is  acquired  by  employing  the  assumptions 
of  an  isopiestic  ideal  gas 
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constant 


(A-3) 


p  T  = 
where:  p  =  density 

T  =  gas  temperature 


SIMULATING  THE  HEAT  TRANSFER 


To  model  the  heat  transfer,  the  augmenter  tube  is  first  divided 
into  axial  segments,  with  the  jet,  augmentation  flow,  and  inner  and 
outer  liner*  sections  lumped  separately  as  illustrated  on  Figure  A- 1 . 
Temperature  is  assumed  constant  throughout  each  section-segment;  flows 
are  perfectly  mixed. 

Mass  flow  rates  through  each  segment  are  determined  from  the  geometry 
of  the  jet;  i.e.,  from  Equation  A-2. 
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=  M  -  AM 

SEC, I  1+1 

augmentation  flow  entrained  by  the  jet  in  segment  "I" 

mass  flow  rate  of  jet  through  segment  "I" 

mass  flow  rate  of  augmentation  air  through  segment  "I* 

density  of  augmentation  air 

radius  of  segment  "I"  of  the  jet 


Conservation  of  energy  is  then  applied  to  each  segment  of  each 
section.  For  the  jet  exhaust  (see  Figure  A-2), 


*Some  augmenters  do  not  have  liners  but  are  of  solid  construction.  For 
this  configuration,  the  inner  and  outer  "liners"  in  the  model  represent 
the  inside  and  outside  of  the  wall. 


•  f  JET.I+l  • 

n JET ,1+1  \  JET.I+l  '  REF  2  /  qRAD, I , JET^IL 


m  L  h  JET,  I 

MJ£T,I  V*  JET , I  hR£F  2 


-  \  SEC ,  I  -  hR£F  +-¥iI 


or,  using  the  ideal  gas  approximation  after  applying  conservation  of 
mass  to  eliminate  the  reference  enthalpies. 
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where:  T  =  temperature  of  jet  in  segment  "I* 
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JP, JET, I 


specific  heat  of  the  jet  exhaust  in  segment 


T  T  =  temperature  of  augmentation  air  in  segment  "I' 
ofct ,  1 


"P,SEC,I 


specific  heat  of  augmentation  air  in  segment 


The  radiation  from  jet  segment  "I"  to  the  wall,  qDAn  T  TI ,  is  deter- 

KAU y 1 ,  J  L  l  I  L» 

mined  by  summing  the  net  radiation  between  this  segment  and  all  sections 
of  the  wall, 
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where:  T 


1L,  J 


JET,1 
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l.J 


£ 
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N 


temperature  of  augmenter  inner  liner  segment  "J" 
surface  area  of  jet  segment  "I" 

Stefan-Boltzmann  constant 

shape  factor  for  radiation  from  jet  segment  "I"  to 
inner  liner  segment  "J" 

£(1,Tjet  j)  =  emissivity  of  jet  segment  "I" 
absorptivity  of  jet  segment  "1" 

£<tjet,i/tii.,j>°’5  <Ref  6> 

number  of  axial  segments  of  the  augmenter  tube 


For  the  augmentation  airflow  (see  Figure  A-3), 
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where:  <lC0Nv ,  I ,  I L^SEt 
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AAUG 


convection  heat  transfer  from  the  inner 
liner  to  augmentation  air  segment  "I" 


hAAUG  (TIL,I 


-  T  ) 

sEc,r 


h  (I,  V  )  =  convective  film  coefficient 

t  i 

of  segment  "I" 


AAUG  I  =  surface  area  of  liner  in  contact 
with  augmentation  air  segment  "I" 


Conservation  of  energy  across  segments  of  the  inner  and  outer  liners  are 
handled  in  a  similar  manner.  For  the  inner  liner, 
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where:  jl  =  conduction  along  the  inner  liner 

1  kAIL(TlL,l  '  TlL,I-tV4Z 

qCOND , I , I L-K)L  =  conduction  between  inner  and  outer  liners 


qRAD  I+J  IL  =  radiation  between  different  segments  o 
'  ’  inner  liner 

qRAD,I-J,IL  =  °  AAUG  FlI,J(1  '  e)  (T f L , I  '  TIL,jj 


f  the 


where:  FI 


lj  j  =  shape  factor  for  radiation  from  inner  liner  segment  "] 
’  to  inner  liner  segment  "J" 

k  =  effective  coefficient  of  thermal  conductivity, 

includes  liner,  acoustic  pillows,  concrete,  etc. 

(where  applicable) 

A  =  cross  sectional  area  of  inner  liner,  acoustic  pads, 


AZ  =  spacing  between  nodes 


Finally,  for  the  outer  liner, 
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where:  q 


CONV ,  I  ,OL-x»  = 


convection  heat  transfer  off  outer  skin  of 
augmenter 


h  A  (T  -  T  ) 
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Equations  A-A  through  A-7  for  all  sections,  along  with  Equation  A-3 
and  the  Equation  of  Continuity,  are  solved  simultaneously  to  determine 
temperatures  throughout  the  augmenter  assembly.  When  the  jet  reaches 
the  wall,  conservation  of  augmentation  flow  energy  is  deleted,  and  the 


other  equations  are  modified  accordingly.  Due  to  the  nonlinearities,  an 
iterative  technique  is  necessary,  and  the  Gauss-Seidel  method,  with 
relaxation,  is  employed  for  this  purpose.  Figure  A-4  illustrates  typical 
predicted  augmenter  temperature  profiles. 

The  heat  transfer  charact eri s t i cs  perhaps  need  further  discussion. 
Radiation  was  modelled  by  considering  each  jet  segment  as  a  separate  gas 
mass  and  applying  the  mean  beam  length  approximation  suggested  by  Hottel 
(Ref  6).  With  this  approximation,  radiation  from  irregular  geometries 
is  related  to  radiation  from  a  gas  hemisphere  by  extrapolation  using  a 
fictitious  dimension  called  the  mean  beam  length. 

Emissivity  of  the  jet  is  a  function  of  its  geometry  (mean  beam 
length),  composition,  and  temperature.  Geometry  is  easily  acquired  from 
the  node  spacing  and  from  Equation  A-2.  Composition  of  each  segment  is 
determined  by  inputting  the  weight  fraction  of  the  radiators,  CO^ ,  1^0, 
and  CO,  at  the  engine  nozzle  and  then  progressively  diluting  the  jet 
with  augmentation  air.  Temperature  is  the  problem;  the  interdependency 
with  temperature  necessitates  including  the  emissivity  calculations  in 
the  temperature  iteration.  The  gas  emissivity  data  of  Hottel  were 
curve-fitted  to  provide  a  continuous  relationship  as  a  function  of  mean 
beam  length,  composition,  and  temperature.  Figure  A-5  shows  typical 
emissivities  calculated  in  this  manner. 

Very  little  work  has  been  done  in  the  area  of  convection  heat 
transfer  to  confined  jets.  Kang  et  al.  (Ref  7)  have  experimentally 
determined  some  Nusselt  numbers,  but  their  tests  are  not  comprehensive 
enough  to  provide  the  foundation  for  any  correlation.  Thus,  conventional 
relationships  (Ref  8)  for  turbulent  flow  through  a  cylinder  had  to  be 
employed,  substituting  the  equivalent  diameter  of  the  augmentation  flow 
for  the  diameter  of  the  cylinder, 
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SEC 
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where:  h 


kSEC 


film  coefficient  of  convective  heat  transfer  between 
augmentation  airflow  and  the  inner  liner 

thermal  conductivity  of  augmentation  air 
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This  convection  coefficient  varies  with  equivalent  diameter  of  the 
augmentation  flow  and,  therefore,  varies  from  segment  to  segment. 
Properties  of  the  augmentation  flow  vary  with  temperature,  requiring  the 
film  coefficient  calculations  to  also  be  included  in  the  iteration. 
Figure  A-6  shows  typical  calculated  coefficients. 


BOILER  SIMULATION 

To  simulate  waterwalls,  the  temperature  of  the  outer  liner  is  set 
equal  to  the  saturation  temperature  of  the  steam;  and  all  heat  resistance, 
except  the  metal,  is  removed  from  between  the  two  liners.  This  is 
equivalent  to  neglecting  the  resistance  to  heat  transfer  on  the  boiling 
water  side  of  the  walls,  a  good  approximation  since  this  resistance  is 
several  orders  of  magnitude  less  than  the  resistance  to  convection  and 
radiation  from  the  jet  exhaust.  The  steam  generation  is  then  determined 
directly  from  a  summation  of  the  energy  transferred  to  the  outer  liner. 

Simulation  of  convection,  water  tube  boilers  is  more  difficult. 
Conservation  of  energy  is  applied  to  the  jet  and  augmentation  flow 
through  the  heat  exchanger,  to  the  water/steam  flow  through  the  heat 
exchanger,  and  to  the  overall  heat  exchanger.  These  relationships  are, 
in  the  same  order, 

%TM  =  MGAS  CP,GAS(TGAS  *  ?STACK^ 
qSTM  =  ”STM  ^fg 
qSTM  =  UA  ATM 


(A-8) 
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where:  M 


GAS 


STM 


jSTM 


sum  of  jet  and  augmentation  air 
steam  generation 

heat  transferred  from  jet  and  augmentation  airflow  to 
the  steam 


Tgas  =  temperature  of  gas  entering  the  boiler 


tstack 

UA 


temperature  of  gas  exiting  the  boiler 

product  of  boiler  overall  heat  transfer  coefficient 
and  heat  transfer  area 
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heat  of  vaporization  of  steam 
logarithmic  mean  temperature  difference 


where : 
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STM 


saturation  temperature  of  steam 


resulting  in  Equations  A-8,  A-9,  A-10  and  the  four  unknowns:  q  _  , 

a  2)lfl 

MSTM’  ^STACK’  an<1  UA"  ^he  relationship  is  determined  by  stipulating 

a  boiler  effectiveness*,  e, 


TGAS  ~  tstack 

T  -  T 
GAS  STM 


If  boiler  effectiveness  is  known,  steam  generation  and  the  other  unknowns 
can  be  calculated.  Furthermore,  if  the  boiler  tube  configuration  is 
also  stipulated,  the  heat  transfer  characteristics  can  be  estimated. 

For  example,  if  the  boiler  is  to  be  of  a  staggered  tube  configuration 
(Ref  8), 


TUBE 


=  0.287 


JET 

3TUBE 


PJET  VJET  °TUBE 
MJET 


*A  value  e  =  0.76  is  assumed  throughout  these  analyses. 
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where : 


TUBE 


=  convection  film  coefficient  between  jet  exhaust  and 
water  tubes 


TUBE 


=  outer  diameter  of  water  tubes 


It  should  be  noted  that  the  resistance  to  heat  transfer  through  the  tube 
walls  and  to  the  boiling  water  are  being  neglected.  With  this  equation, 
an  estimate  of  boiler  heat  transfer  area  is  also  acquired  because  UA  is 
known.  Properties  are  determined  by  assuming  the  jet  exhaust  behaves 
identically  to  air  and  applying  the  Eucken  equations  (Ref  9), 
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ACCURACY  OF  THE  AUGMENTER  MODEL 

To  evaluate  the  simulation,  predicted  augmen ter  inner  liner  temper¬ 
atures  are  compared  with  measured  values.  The  accuracy  of  these  temper¬ 
atures  makes  a  valid  benchmark  since  they  indicate  the  accuracy  of  the 
simulation  of  heat  transfer  to  the  walls  and,  by  supposition,  to  a 
waterwall  heat  exchanger.  The  same  reasoning  applies  to  the  augmenter 
exit  gas  temperature,  but  here  the  modelling  involves  only  simple  mixing 
of  gases,  and  little  error  is  likely.  Regardless,  these  temperatures 
can  be  examined  concurrently  by  noting  that,  due  to  high  wall  resistance, 
there  is  little  difference  between  gas  and  wall  temperatures  near  the 
augmenter  exit.  This  trend  is  illustrated  in  Figure  A-4. 

An  abundance  of  data  are  not  available  with  which  to  make  this 
comparison.  The  hush  houses  at  the  NAS  Miramar  have  been  studied  exten¬ 
sively.  These  facilities  were  used  to  develop  the  empirical  model  of 
the  jet  exhaust,  however.  It  is  suspect  to  use  the  same  sites  to  evaluate 
the  model.  Nevertheless,  they  will  be  used.  The  hush  house  at  the 
Marine  Corps  Air  Station  El  Toro,  although  studied  less  extensively, 
provides  a  better  accreditation. 


The  comparisons  are  made  on  Figures  A-7  through  A-9.  Predicted  and 
measured  temperatures  are  within  ±50°F  of  each  other.  The  exception  is 
the  TF10  engine  undergoing  afterburning  tests  where  the  simulation  of 
inner  liner  temperatures  is  in  error  by  more  than  200°F.  Not  enough 
temperatures  have  been  measured  near  the  augmenter  inlet  to  verify  the 
existence  of  the  radiation  peaks  shown  on  the  figures. 

Several  sources  of  error  are  obvious,  each  related  to  the  assump¬ 
tions.  First,  the  hush  houses  are  not  cylindrical  but  elliptical. 
Although  the  distance  from  the  nozzle  to  the  nearest  augmenter  wall  was 
used  as  an  effective  radius,  the  effects  of  this  cross  section  on  flow 
characteristics  could  not  be  accounted  for.  The  steady  state  assumption 
is  possibly  not  valid  for  the  afterburning  tests  which  tend  to  be  of 
short  duration.  Actual  inner  liner  temperatures  could  be  much  higher  or 
much  lower  than  steady  state  temperatures,  depending  upon  wall  thickness, 
specific  heat,  and  thermal  conductivity  as  well  as  time.  Jet  exhaust 
and  secondary  entrained  air  are  certainly  not  perfectly  mixed  at  the 
nozzle.  A  hot  core  of  exhaust  will  survive  for  some  distance,  tending 
to  damp  out  both  the  radiation  peak  and  convection  valley  predicted  by 
the  model  and  illustrated  on  Figures  A-7  through  A-9.  Finally,  the  jet 
nozzle  is  not  located  exactly  at  the  augmenter  inlet  but  some  10  feet 
back  inside  the  hush  houses. 
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Figure  A-2.  Conservation  of  energy  across 
segment  of  jet  exhaust. 
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Figure  A-3.  Conservation  of  energy  across 
segment  of  augmentation  air. 
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Figure  A-4.  Typical  predicted  temperature  profiles  along  the  augmcnter. 


Notes:  (1)  NAS  Miramar  hush  houses 

(2)  J79  at  MAX  A/B 
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Figure  A  5.  Typical  emissivities  of  jet  exhaust. 
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(2)  J79  at  MAX  A/B 
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Figure  A-6.  Typical  coefficients  of  convective  heat  transfer  between 
inner  liner  and  augmentation  airflow. 
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Figure  A-7.  Evaluation  of  augmenter  model  by  comparison  with  temperatures 
measured  during  J  79  testing  in  the  hush  houses  at  NAS  Miramar. 
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Appendix  B 


METHOD  OF  CALCULATING  THE  ECONOMICS  OF 
TEST  CELL  ENERGY  RECOVERY 


Once  test  cell  steam  generation  has  been  determined,  its  value  is 
estimated , 


VALUE  =  Motm  x  COST  x  TIME  (B-l) 

a  FM 

where:  M„„u  =  mass  flux  of  steam  generated 

bln 

COST  =  assumed  value  of  one  unit  of  steam 

TIME  =  average  utilization  of  the  test  cell  (see  Table  C-l) 

and  accumulated  over  the  economic  life  of  the  heat  exchanger  (or  life  of 
the  test  cell,  if  shorter), 
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NN 


=  present  value  of  steam  generated  over  a  period  of  NN 
ESC  =  assumed  escalation  rate  of  the  cost  of  steam  (Ref  5) 
DISC  =  assumed  discount  rate  (Ref  5) 
allowing  a  savings/investment  ratio  to  be  predicted, 
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where  PRICE  is  the  cost  of  the  waterwalls  or  heat  exchanger. 

The  economics  of  the  Rankine  cycle  generation  of  electricity  are 
analogous,  substituting  the  value  of  electricity  for  the  value  of  steam. 


COST  OF  ENERGY  RECOVERY  HARDWARE 

A  major  obstacle  is  the  determination  of  PRICE.  Very  few  boilers/ 
heat  exchangers  large  enough  for  this  application  are  manufactured,  and 
these  few  are  all  custom  jobs. 


To  acquire  an  estimate  of  the  cost  of  a  test  cell  heat  exchanger, 
twenty  boiler  manufacturers  were  asked  to  provide  a  preliminary  design. 
Table  B-l  describes  a  water  tube  heat  exchanger  composited  from  their 
responses.  It  is  presented  without  documentation.  For  purposes  of 
applying  this  design  to  a  particular  test  cel  1 /engine/power  setting,  a 
scaling  criterion  was  developed. 

The  downstream  convection  heat  exchanger  (Configuration  No.  1  of 
Figure  1)  cost  is  ratioed  linearly  from  Table  B-l  using  the  total  gas 
flow  as  the  independent  variable. 


PRICE 


PRICE 


B1 


TOTAL 


TOTAL, B1 
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where  the  subscript  B1  refers  to  the  values  in  Table  B-l  and  TOTAL 
refers  to  the  sum  of  the  jet  and  secondary  flow.  The  overall  heat 
transfer  coefficient  is  kept  at  approximately  9  Btu/hr-ft2oF,  equivalent 
to  a  flow  velocity  of  about  80  ft/sec  over  2-inch  tubes. 

When  convection  heat  exchangers  are  inside  the  augmenter  near  the 
engine  nozzle  (Configuration  No.  A  of  Figure  1),  the  heat  transfer  area 
is  input  directly  while  the  overall  heat  transfer  coefficient  is  allowed 
to  increase,  following  the  high  velocities  encountered.  The  price  of 
the  heat  exchangers  is  established  in  a  manner  analogous  to  Equation  B-A, 
using  heat  transfer  area  in  place  of  total  gas  flow, 
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This  relationship  is  applied  for  lack  of  a  better  approach.  The  high 
gas  velocities  would  be  expected  to  create  some  serious  tube  erosion 
problems,  requiring  thicker  tubes  and  probably  special  alloys.  In 
addition,  the  tubes  would  have  to  withstand  occasional  bursts  of  A/B. 

It  follows  that  these  heat  exchangers  would  be  more  expensive  than 
equivalent  downstream  configurations. 

No  waterwall  designs  were  submitted  in  response  to  the  survey  of 
boiler  manufacturers.  Waterwall  areas  are  input  and  heat  transfer 
characteristics  built  into  the  model,  howeve. ,  neither  presents  a  problem. 
The  use  of  waterwall  prices  will  be  avoided,  to  the  extent  possible,  in 
the  comparisons.  When  necessary,  Equation  B-5  will  be  employed.* 


*Again,  for  lack  of  a  better  approach.  This  is  tantamount  to  the 
assumption  that  the  fabrication  of  the  tubes  is  the  major  expense 
in  the  manufacture  of  each  device. 
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Table  B-l.  Composite  Heat  Exchanger  for  Energy  Recovery 
From  Gas  Turbine  Engine  Test  Cell 


Type 

Convection,  Shell  and  Tube 

Gas  flow 

2.3  (10)6  lb/hr 

Gas  temperature  in 

420°F 

Gas  temperature  out 

367°F 

Steam  flow 

30,300  lb/hr 

Steam  temperature 

350°F 

Steam  pressure 

135  psia 

Heat  transferred 

30  (10)6  Btu/hr 

Heat  exchanger  width 

50  ft 

Heat  exchanger  height 

15  ft 

Heat  exchanger  length 

72  ft 

Tube  diameter 

2  inches 

Tube  thickness 

0 . 12  inch 

Fin  segment 

0.156  inch 

Fin  length 

0.5  inch 

Fin  thickness 

0.035  inch 

Number  of  fins 

6/in . 

Number  of  circuits 

45 

Number  of  passes 

12 

Heat  transfer  surface 

118,000  ft2 

Gas  AP 

5.85  inches  H2O 

Cost 

$2,000,000. 
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Appendix  C 

JET  ENGINE  TEST  CELL  CHARACTERISTICS 


Table  C-l  summarizes  the  geometries  and  flow  characteristics  of 
test  cells  used  in  the  analyses.  With  the  exception  of  the  hush  houses 
at  the  MCAS  El  Toro,  all  of  these  facilities  are  located  in  the  San 
Diego  area.  Geometries  and  cell  utilization  were  acquired  from 
Reference  10.  This  reference  also  describes  the  cells  in  greater  detail. 
Flow  rates  and  temperatures  were  acquired  from  References  3  and  U.  The 
compositions  of  the  engine  exhaust  gases  at  the  nozzle  were  determined 
by  back-calculating  from  stack  gas  compositions  measured  on  the  Miramar 
"A"  test  cell  and  NARF  North  Island  Test  Cell  19  (Ref  11,  12  and  13). 

It  was  assumed  that  the  chemical  composition  of  the  jet  exhaust  at  the 
nozzle  was  dependent  only  on  the  type  of  engine  and  power  setting;  the 
same  values  were  used  for  the  other  test  cells. 


C-l 
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Appendix  D 

COMPUTER  MODELLING  OF  AUGMENTER  TUBE 


The  relationships  in  Appendixes  A  and  B  were  programmed  in  FORTRAN 
IV.  Figure  D-l  is  a  flow  chart  of  the  program.  A  listing  and  some 
typical  results  follow. 


D-l 


j 


Figure  D-l.  Flowchart  of  computer  modeling  of  augmenter  tube 
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PROGRAM  LISTING 


PROGRAM  JETEST  ( I NPUT * OU TPUI , T APE 5 = I NPUT , T APES =OU TPUT > 

C  PROGRAM  TO  CALCULATE  THE  HEAT  TRANSFER  FROM  A  JET  EXHAUST  TO  THE 
C  TEST  CELL  AU3MENTER  WALLS  OR  TO  BOILERS  WRAPPED  AROJND  THE  WALLS 
C  OR  TO  BOILERS  PLACE!  IN  THE  JET  STREAM  AT  THE  AUGMENTER  EXIT  OR  BOTH 
C 

C  STEADY  STATE  AXISYMMETRIC  FLOW  THRU  A  CYLINDRICAL  DUCT  IS  ASSUMED 
C  IN  ADDITION,  IT  IS  ASSUMEO  THAT  BOTH  THE  JET  AND  SECONDARY  FLOW  ARE 
C  TURBULENT,  PERFECTLY  MIXED  IDEAL  GASES 
C 

REAL  MJETC51) , MSEC (51 ) ,K DUCT , KAI NT , K AEXT , KP AD, N2C 50 ) , L C 50 ) , K AI R  , 

1  MU A  IR ,MSTM1C50 >  ,K 
INTEGER  TYPE 

COMMON  02,N2,C02,H20,CO,H2,LUM 

DIMENSION  TC  50  *  4) , TEMP C 5 0 , 4 ) , HCO NV C 50  >,ARADC50),rC50,50),DELMC50> 
1,02C50) »AJETC50),RJET  <50 ) , RAD S T < 5 0 ) , AB 0 1 LC 2 ) , TSTM C 2 ) , PSTM C 2 > 
2,HSTM(2),RADIN(50) ,RADOUTC50> ,Z(50) ,C02 ( 5 0 ) , H20 C 5 0 1  ,COC50) ,H2C50) 
3,3CONVC5D> , EMI  SC  50) ,F1C50,50> ,UJETC51  ) 

DIMENSION  PVC25),SIRC25)  ,P18C50> ,P43C50>  »P33C50> 

OAT A  PI ,PRESS»RAIR, SIGMA, TQATUM/3.1416,2116. ,53.36,0. 1 714E-3  ,520. / 
pATA  I  TER , T  OL /I 00 *1./ 

OATA  R/-0.5/ 

DATA  EMPIR1,EMPIR2/14.,14./ 

C 

C  I N3  J  T  THE  NUMBER  OF  AXIAL  SEGMENTS  THAT  THE  AUGMENT  ER  TUBE  WILL  BE 
C  01 V l DED  INTO,  N.LE.50 
READ! 5,1)  N 

C  INPJT  THE  AUGMENTER  GE OME TR Y .  . . . THE  AREA  OF  THE  JET  NOZZLE  AND  THE 
C  AREA  OF  THE  AUGMENTER  INLET,  BOTH  IN  SOFT. ...THE  LENGTH  OF  THE 
C  AUGMENTER  IN  FEET. ...AND  THE  THICKNESSES  OF  THE  DUCT  AND  THE  ACOUSTIC 
C  PILLOWS  IN  INCHES 

RE ADC  5, 2 1  A JET  Cl) , A AUG, AUGL, 0 ELR , DELP AD 
C  I N3 J T  SECONDARY  FLOW  RATE  IN  LBS/SEC  AND  ITS  TEMPERATURE  IN  DESF 
REA DC  5, 3)  MSEC  Cl) ,  T  f  1  ,  2  ) 

C  INPJT  JET  FLOW  RATE  IN  LBS/SEC,  TOTAL  TEMPER  AT  JRE  AND  TOTAL  PRESSURE 
C  OF  JET  IN  OEGF  AND  PSIA,  AND  RATIO  OF  JET  SPECIFIC  HEATS 
R E ADC  5,4)  MJETC1) , TT JE T , PT JET , GA MMA 
C  INPJT  THE  CHEMICAL  COMPOSITION  OF  THE  JET ...  .PERCE  NT  WEIGHT  OF 
C  OXYSEN,  NITROGEN,  CARBON  DIOXIDE,  CARBON  MONOXIDE,  WATER  VAPOR,  AND 
C  HYDROGEN 

READC  5,6)  02C1) ,N2C1) , C  0  2  C 1 ) , C  0  C 1 ) , H  2  0  C 1 )  ,  H2  C 1 ) 

C  INPJT  THE  LENGTH  OF  THE  LUMINOUS  PDRTION  OF  THE  JET  EXHAUST  IN  FEET 
REA0C5,3>  FLAME 

C  INPJT  THE  COEFFICIENT  0=  THERMAL  CONDUCTIVITY  OF  THE  DUCT  WALLS  AND 
C  THE  ACOUSTIC  PILLOWS  IN  B TU /HR-F T -DE GF . . .  .  I F  THE  PILLOWS  ARE  NOT 
C  USED,  INPUT  OELP AD  =  0 •  AND  KPAO^IOOOO. 

READC 5, 7)  KDUCT ,KPAD 

C  INPJT  BOUNDARY  CONDITIONS  TO  SIMULATE  THE  HOUSE  STRUCTURE  ACTING  AS 

C  A  HEAT  SINK . THE  EFFECTIVE  TEMPS  AND  THERMAL  CONDUCTANCES,  IN  DEGF 

C  AND  8 TU/ HR-QEGF ,  AT  THE  AUGMENTER  INLET  AND  EXIT,  IN  THAT  ORDER 
READC5»4>  TINT, KAINT, TEXT, KAEXT 

C  INPJT  AMBIENT  COND I T I ONS . . . .T HE  TEMP  IN  DEGF  AND  THE  NATURAL 
C  CONVECTION  HEAT  TRANSFER  COEFFICIENT  IN  B TU/ HR- SOFT- DEGF 
REA0C5, 3)  T  AMB»HAMB 

C  INPJT  THE  HEAT  EXCHANGER  C  B  OILER )  CONFIGURATION 
C  "1"  SIGNIFIES  NO  HEAT  EXCHANGERS 

C  "2"  SIGNIFIES  WATERWALLS  WRAPPED  AROUND  THE  AUGMENTER  TUBES 

C  "3"  SIGNIFIES  CONVECTION  BOILER  IN  THE  AUGMENTER  GASES  AT  EXIT 
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c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 

c 

c 


c 

c 


c 

c 


"4"  SIGNIFIES  BOTH  WATERUALLS  AND  CONVECTION  BOILERS 
REA0C5,1)  TYPE 
IF  (TYPE.EQ.l)  GO  TO  10 
INPUT  THE  APPLICABLE  80ILE*  CHARACTERISTICS 

INPUT  THE  TOTAL  NUMBER  OF  HRS  PER  DAY  THAT  THE  TEST  CELL  IS  UTILIZED 
WITH  THIS  ENGINE  AND  AT  THIS  POWER  SETTING  AND  ALSO  INPUT  IND2.GE.7, 
FOR  A  PARAMETRIC  STUDY  INPUT  TIME=0.25  AND  I ND2  =  0 .  .. . I N3UT  PURCHASE 
PRICE  OF  THE  WATER-WALL  AND  CONVECTION  HEAT  RECOVERY  HEAT  EXCHANGERS 
IN  THAT  ORDER.. ..IF  THERE  IS  ONLY  ONE*  INPUT  ZERO  IN  THE  SPACE 
FOR  THE  OTHER 


RE ADC  5. 73)  T I  ME *1 ND2* PR  I CE 1  * PR  ICE 2 

INPUT  THE  CURRENT  COST  OF  THE  STEAM  IN  S/MBTU,  THE  DISCOUNT  RATE*  AND 
THE  LONG  TERM  ESCALATION  RATE  OF  THE  COST  OF  STEAM  IN  PERCENT 
PER  YEAR 

REA DC  5*4)  COST  *DI SC  *ESC 
IF  CTYPE.EQ.3)  GO  TO  5 

IN3 JT  THE  WATERWALL  CH AR ACT ER I S T I C3 . . . . S TE A M  TEMP  IN  DE3F*  STEAM 
PRESSURE  IN  PS  I A .  HEAT  OF  VAPORIZATION  IN  BTU/LB ♦  AND  THE  TOTAL  HEAT 
TRANSFER  AREA  IN  SOFT 

READC5*4)  TSTM(l) ,PSTM(1 )»HSTM(1)»AB0IL<1) 

TSTMC1)=TSTMC1)*460. 

INPUT  THE  LOCATION  OF  THE  U A TER- W ALLS.  .  . . FR ONT  AND  BACK  IN  TERMS  OF 
FEET  FROM  THE  AUGMENTER  INLET 
READC5,3)  FRONT, BACK 
IF  CTYPE.EQ.2)  GO  TO  10 

IF  BOTH  BOILERS  ARE  UTILIZED,  INPUT  WATERWALL  CHARACTERISTICS  FIRST 
I N3  J  T  THE  CONVECTION  BOILER  CH AR ACT ER I S T I CS . . . . S TE AM  TEMP,  PRESSURE* 
AND  HEAT  OF  VAPORIZATION,  TOTAL  HEAT  TRANSFER  AREA,  AND  THE  OUTER 
DIAMETER  OF  THE  TUBES  IN  INCHES 
5  READ(5,8)  TSTM t 2 ) ,PST M C 2 ) »HST MC 2 ) » AB 01 L < 2  )  * DTUBE 
TSTMC2)=TSTMC2)«-460. 


PRINT 

ALL 

DIMENSIONS 

10 

UR 

ITE 

CS 

,11) 

UR 

ITE 

CS 

,12) 

AJETC 

UR 

ITE 

CS 

,13) 

UR 

ITE 

CS 

,14) 

MJETC 

UR 

ITE 

CS 

,1S) 

UR 

ITE 

CS 

,17) 

02  Cl) 

UR 

ITE 

C  6 

,18) 

UR 

ITE 

C  6 

,19) 

MSECC 

UR 

ITE 

CS 

,67) 

UR 

ITE 

CS 

,68) 

TAMB, 

WRITE 

CS 

,21) 

UR 

ITE 

CS 

,22) 

KDUCT 

IF 

CDELPAD . 

LE.O.) 

UR 

ITE 

CS 

,23) 

KPAD 

15 

UR 

ITE 

CS 

,24) 

UR 

ITE 

CS 

,26 ) 

TINT, 

IF 

CT 

YPE.EQ 

.1)  UR 

*  PROPERTIES,  AND  BOUNDARY  CONDITIONS 
1 1 , A  AUG , AUGL ,DELR,DELPAD 
1) ,TTJET*PT JET, GAMMA 

, N2C 1 ),C02(1 ) ,C0C1>,H20C1),H2C1) ,FLAME 
1 ) ,T  Cl ,2) 

HAMB 

GO  TO  15 

KAINT,TEXT,KAEXT 
ITEC  6,27) 


SIMPLIFYING  AND  PRELIMINARY  CALCULATIONS 
DELR=DELR/12. 

DELPAD=DELPAD/12. 

DTUBE=OTUBE/12. 

DELZ=AUGL/CN-1 ) 

LUM= I  NT C FLAME/DEL Z) 
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<RAUG  =  SQRT< AAUG/PI  ) 

R  JET<  1  ) =SQR  T  <  A  JE  T  (1  )/PD 
AREAR=2.*PI *RAUG*OELZ 
AREAZ=PI*C  CRAUG+OELR) **2-RAUG**2> 

K=<OELR*DELPAO> /< DELR/KDUCT*DELP AD/KP AD) 

PTJET=PT  JET  *144. 

T<1 *2)=T<1 ,2)*460. 

TTJET=TTJET ♦460. 

TINT=TINT*460. 

TEXT=TEXT*4G0. 

TAMB=TAMB*460. 

RATIO  =  « 1 •♦ESC/1  00 .1/(1. *01  SC/1  00.) 

3l=K*AREAR/2./(DELR*OELPAD> 

»2=AREAR/2. 

P3=K0UCT*AREAZ/DELZ 
?4=K/f DELR*DEL?AO> 

?5-2. *P1 

?6=KDUCT/DELZ 

?7=SIGMA*AR£AR 

IF  IITYPE.EQ.l>.OR. {TYP£.EQ.3>>  30  TO  20 
I  FRONT  =  INT  (FRONT/DELZ  )  *1 
IBACK=INT<3ACK/0ELZ)*1 
IF  (  I  SACK  •  G  T  •  N  >  I8ACK  =  N 
C 

C  CALCULATE  BULK  VELOCITIES  AT  AUGMENTED  INLET 
20  RHOS=PRESS/RAIR/T(l *2) 

VSEC=MSEC<1 l/RHOS/t  AAUG- AJET< 1  )  ) 

P3UER=<GAMMA-1 .) /GAMMA 

7<1 ,1 )=TTJET*(PRESS/PTJET)**POWER 

RHOJ=PRESS/RAIR/T  tl  ,1 > 

VJET=MJET(l)/RHOJ/AJET(l  ) 

UJETt 1 >-VJE T 
UJETt  2 1 =UJET  <  1 ) 

C  CALCULATE  THE  CRAYA-CURTET  NUMBER 
A=AJET(1)/AAUG 
8=1. -A 
C=0.5-A 

CRAYA=tA*VJET*B*VSECJ/<A*VJET**2*C*VSEC**2-0.5*tA*VJET*R*VSEC>**2  > 
1**0.  5 

4RITEC6.34I  VJET  t  VSEC * CR  AY  A 
C 

C  CALCULATE  THE  EXPANSION  OF  THE  JET  USING  THE  EMPIRICAL  RELATIONSHIPS 
C  OF  BECKER*  ET  AL*  OTH  SYMP.tINT)  ON  C0MB.»l?o3 
C  THE  CONSTANTS  HAVE  BEEN  ADJUSTED  TO  FIT  TEST  CELL  DATA 
FACT0R=0. 

IF  CCRAYA.GT.0.7)  GO  TO  40 
XR=4.07*£XPt3.54*CRAYA> 

FACTOR=l./XR**l.&67 
A=SQRTC A> 

X=AUGL/RAUG 
DE  V-X 

DO  30  1  =  1 »  ITER 
DSAVE=DEV 
XSA VE=X 

FF=0.13l*(FACTOR*X**2.667*X)-l ,*A 
D':F0X  =  0.34S*FACTGR*X«*1.667*0.13l 
X=X-FF/OFFOX 


DEV=A3S(X-XSAVE) 

IF  (OEV  .LE.DSAVE)  30  TO  ?5 
WRITE(6,36) 

WRITK6.37) 

WRI TE(6,38) 

30  TO  500 
25  CONTINUE 

IF  (DEV. LE. TOD  GO  TO  35 
30  CONTINUE 
WRITE(6,36) 

WR  I  TE  (  6  »  39  )  TOL  » I TE  R 
WRITE (6*38) 

35  XSEP=X*RAUG 
30  TO  45 

40  XSEP=7.91*(RAUG-RJET(1>) 

45  CONTINUE 

C  DETER  MIME  IF  JET  EVER  REACHES  AUGMEMTER  WALL 
IF  (XSEP.GT.AUGL)  GO  TO  50 
\ISEP=INT(XSEP/OELZ)*2 
IF  (NSEP.GT.N)  NS  EP  =  N 
WRITE (6*41)  XSEP 
SO  TO  55 
50  M5EP-10000 
WRITE(6,42) 

55  CONTINUE 
C 

C  CALCULATE  THE  PROGRESSIVE  RADII  OF  THE  JET  EXHAUS T . . . . A3 AI N  USE  THE 
C  EMPIRICAL  RELATIONSHIPS  OF  BECKER... 

Z( 1)  =  0. 

DELHI  1)=0. 

DO  60  I =2  *  N 
Z( I >=(I-1)*DELZ 

RJETf I) =RJET(  l)+0 .1 31 *Z( I  )  *  (  1  .  «-F  AC T OR  *  ( Z  (  I  > /RA UG  )  *  *  1 . 667) 

IF  (CRAYA.GT.0.7)  R JE T ( I ) =R JET ( 1 ) *0 . 1 2S5*Z ( I ) 

IF  (RJETf I).GT.RAUG)  RJET(I)=RAUG 
C  CALCULATE  SECONDARY  FLOW  ENTRAINED  BY  THE  EXPANDING  JET 
DE LM< I )=RHOS*VSEC  *PI* (RJETf I ) • *2 -RJETf 1-1 )  **2) 

60  CONTINUE 
<K=N-1 

DO  490  1=1 »KK 

MJETC I*l)=MJETf I) ♦  DELMf I ) 

MSECf  I«-1)=MSEC(  I) -DELMf I +1) 

490  CONTINUE 

MJETfN  +  DsMJETfN) 

WRITE (5*9) 

11  =  1 

DO  65  1=1*5 
L I M  =1  0  *  I 

I3RINT=MINO(LIM,N) 

WRITE(6»43) 

WRITE (6*44)  (Z(J) » J  =  I I , I  PR  I NT ) 

WRITE(6»4S)  (RJETf J>* J=I  I*IPRINT) 

WRITE(S*47)  ( MJET (  J*1 ) * J  =  1 1  *  I  PR  I  NT  ) 

WRITE (6 *49)  (MSEC< J) * J=IT  ,  I  PR  I  NT ) 

IF  (IPRINT.GE.N)  GO  TO  70 
11= I3RINT*l 
6S  CONTINUE 
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70  CONTINUE 


C 

c  calculate  radiation  constants... .areas,  mean  beam  lengths,  and 

C  SHA  3E  FACTORS 
DO  75  I  =1  ,  N 

ARADC I D=2.*PI*RJETC I) *DELZ 
LC I >  =  1. D*R JET (  I  ) 

75  CONTINUE 

ARADCl>=ARADClD/2. 

AR ADC N) = AR AD  f  N ) /2 . 

CALL  SHAPECRJET  ,RAUG»DELZ»F»N  > 

CALL  SHAP£1CRAUG,QELZ  ,  F 1 *N  D 
C 

C  CALCULATE  THE  COMPOSITION  OF  THE  JET  EXHAUST  IN  MOLES/LB  OF  EXHAUST 
02C 1) =02( 1  )  /3200. 

N2ClD=N2ClD/2800. 

C0-2(1)=C02(1)/4400. 

H20C1 ) =  H2  0  C 1 D  /  1 80  0  • 

C0C1D=C3C1  )/2800. 

H2C 1) =H2C 1 )/200. 

DO  80  I  =2  »  N 

02  C  I)  =02<  1-1  D*0  .0072»DELMC  D/MJETC  I) 

N2C  ID  =N2C I-1D*0.0272»DELMC ID/MJET  C ID 
C  02 C I D=C02C  I-1D *  M  J  E  T  < 1-1 D/MJETC I > 

H20CI )=H20C 1-1 ) *MJETC 1-1 D/MJETC I D 
C0CID=C0CI-1D  *MJETC  1-1 D /MJETC I D 
H2C ID=H2C I-1D*MJETC 1-1 D/MJETC I D 
80  CONTINUE 
C 

C  USE  G A  US  S- SE I  DEL  ITERATION*  WITH  RELAXATION.  TO  SIMULTANEOUSLY  SOLV 
C  FOR  JET  AND  SECONDARY  FLOW  TEMPERAT JRE3  AS  tJELL  AS  TEMPERATURES  ON 
C  THE  INSIDE  AND  OUTSIDE  OF  THE  AUGMENTER  WALL 
C 

C  FIRST,  INITIALIZE  TEMPS  FOR  THE  ITERATION 

r3UESS=CMJETClD*TCl»lD+MSECCID*rCl»2D D/CMJETC1D*MSECC1DD 
TCI  ,3»  =  TSUESS 
TCl,4D=TAMB 
DO  85  1  =  2, N 

TCI ,1 >=TC1,ID-CTC 1 .ID -TGUES3D  »ZC I D/AUGL 
TCI,2)=TC1 ,2> 

TCI »  3  D  =T  SUE33 
T  C I ,4  D  =  TAM8 

85  CONTINUE 

C  AND  .MAKE  SOME  ADDITIONAL  PREPATORY  CALCULATIONS 
DO  86  I = I , N 
MJETC I D  =  5600. *MJET  C I D 
MSECC I D  =3600.*MSECC ID 
DELMC I)=3600.*Q£LMCI) 

UBCI  D=P1 
3 4 B  C l D  =P4 
358C I D  =P5 

IF  CCTYPE.E3.1D.0R.CTYPE.E3.3DD  GO  TO  86 
Ir  CC I .LT.  IFRONTD  .OR. C I  .GT  .I3ACKD D  GO  TO  3S 
P1BCI D=KDUCT*AREAR/2./CELR 
34B  C I D =KOUCT/DELR 
a5BCID=2.*PlHCID 

86  CONTINUE 
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MJET(N«-1)=MJET(N) 

MSEC(  N*  1 ) =M$EC (  N) 

IND1=MIN0(N,N3EP) 

IND3=MIN0(N-1 * NSEP-1) 

JK=NS£P-1 
DE V  MAX=  T( 1,1) 

C 

C  BEGIN  THE  ITERATION 
00  AO  0  1  =  1  *  ITER 
D3AVE=DEVMAX 
DEVMAX=0. 

C 

C  CALOJLATE  THE  CONVECTIVE  HEAT  TRANSFER  FILM  COEFFICIENTS 
C  USE  THE  RELATI ONSHIPS  OF  MCADAMS  FOR  TURBULENT  FLOW  THRU  A  DUC T  * 

C  HOWEVER*  USE  BOTH  VELOCITIES  AND  THE  EQUIVELENT  DIAMETER  OF  THE 
C  SECONDARY  FLOW  "DUCT "  TO  CALCULATE  REYNOLDS  NO. ...THE  CONSTANT  IS 
C  EMPIRICALLY  ADJUSTED  TO  FIT  TEST  CELL  TEMP  DATA 
DO  90  11=1* IND3 

- ACTOR  =  T< II *21 **1.5/(225. *T( 1 1 ,2) ) 

<AIR=9.26E-A*FACT0R 
MUAIR=2.694E-3*FACT0R 
RHOAIR=PRESS/RAIR/T (I  I *2  ) 

VC0NV1=VSEC*T( II, 2)/r (1 ,2) 

RHOJ=PRESS/RAIR/T (11,1) 

VC0NV2=MJET(II)/RH0J/PI/RJET( II)**2/36Q0. 

VCONV  =  A*VCONV2-*B*VCONV1 
DEQV  =  2.*(RAUG-!RJET(  II  ) ) 

RCONVC I  I )=EMPIR1*KAIR/DEQV**0 .2  * ( VC ON V *R H 0 A  I R / HUA I R > *  *  0 . 5 
HCONV( II)=HCONV(I I  )/II 
90  CONTINUE 

IF  (NSEP.GT.N)  GO  TO  100 

C  ONCE  THE  JET  REACHES  THE  WALL,  THE  PROBLEM  SIMPLY  BECOMES  THAT 
C  OF  TURBULENT  FLOW  THRU  A  CIRCULAR  DUCT 
DO  95  1 1  =  I ND1 , N 

-  AC  TOR  =  11,1 1 **1.5/(2  25. ♦Til  1,1)) 

KAIR=9.26E-4*FACT0R 
MUAIR=2.&94E-3*FACT0R 
RHOAIR  =  PRESS/RAIR/T(II,l  ) 

VCONV=( MSEC(l) ♦MJETIl ) )/AAUG/RHOA IR/3600 . 

HC0NV(II)=EMPIR2*KAIR  * ( V CONV* RHOA I R/ MUA I R ) *  *  0 . 3/R AUG* *0 .2 
95  CONTINUE 
100  CONTINUE 
C 

C  CALCULATE  BULK  VELOCITIES  OF  THE  JET 
DO  102  11=2, N 
RHOJ=PRESS/RAIR/T (11,1) 

AJET(II)=PI*RJET(II)**2 

J  JET  (11*1)  =  MJET  (I  I*l)/RHOJ/AJET(  ID/36  00. 

102  CONTINUE 
C 

SAVE=T(1,3) 

RADINC  1)  =  0. 

DO  105  I I = 1 , N 

EMSS=EMISS(T(II,1 ),L( II) ,11) 

ALFA  =  EMSS*SQRT(T(II,1)/T(1  ,3)  ) 

C  RADIATION  REACHING  THE  AUGMENTER  WALLS  HAS  TWO  SOURC ES . . . F ROM  THE 
C  EXHAUST  GAS  OR  FROM  OTHER  PARTS  OF  THE  WALL 
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C  NES.ECT  THE  FACT  THAT  SOME  RADIATION  BE  TWEEN  WALL  SEGMENTS  DOES  NOT 
C  PASS  THROUGH  THE  JET 

TAOINTl )=RAOIN( 1) +SISMA* ARAO( II ) *F( 1 1 tl > *  t EMSS*T< II,1)**4-ALFA 
1*T<  1 , 3)  **4  >«-P7*Fl  (I  I,  1)  *  <1  .-EMSS)  *<T  (  II,3>**4-T«1,3>**4> 

105  CONTINUE 

T  (1  ,3)-<PlB(l)*T<  l,4)-*-P2*HC0NVll)*T(l,2)+KAINT*TINT-»-33*T(2,3> 

1  «-RAOIN<  I)  )  /  (  °1 8  Cl  )  «-P2  *  HC  ONV  (  1  )  «-P3*KA  I  NT  > 

T  <  1  » .3  I  =  T<  1  *  3 )  *R  *<  T <  1  ♦  3  J  -  SA  VE  ) 

DEV=ABS(T(1,3)-SAVE) 

I-  (DEV-DEVMAX)  115,115, 110 
110  JEVMAX=OEV 

115  CONTINUE 

DO  125  J=  1 »  N 
3AVE=T< J,4> 

IF  <<TYPE.EQ.2>.0R.CTYPE.EQ.4>)  GO  TO  115 
r<J»A)  =  CHAMR*TAMB  +  PA*TCJ,3))/(HAMB«-:>A) 

33  TO  117 

116  IF  It J.GE. IFRONT) . A  NO . ( J . LE . I B AC K > )  T < J , A )  =  (  1 0 0 0 . *T ST M ( 1 ) +  P6 * 

1  T  <  J  »  3  )  )/(  1000  . *-P6) 

117  T<J,4)=T<J,4),-R*<T<J,4  >-SAVE> 

DEV=ABS(T (  J  ,  4  ) -SAVE  ) 

IF  (OEV-OEVMAX )  125,125,120 
120  DE V  MA  X=  DEV 
125  CONTINUE 

DO  140  J  =  2  » I  NO  3 
S  A  VE  =  T  (  J « 1 ) 

R  ADO'JTC  J>  =  3. 

EMSS=EMISStT<J»l)  ,L<J),J) 

DO  130  1 1  =  1 »  N 

ALFA=EMSS*SQRTCT< J,1)/TCII  ,3)  ) 

C  IT  IS  ASSUMED  THAT  THERE  IS  NO  RADIATION  BETWEEN  OI-'FERENT  LOCATIONS 
C  ALONG  THE  JET  EXHAUST  STREAM 

RADOUTC  J1=RADOUT<  J)  t-SIGMA*ARAD  <  J>*F(  J , 1 1 > M £ MSS *T C J  ,1  )  **4-ALFA 
1*T C II ,31**4) 

130  CONTINUE 

CPIN=SPHT(T( J-l ,1 ) »J-1  ) 

“POUT  =  SPHT (T  f  J  « 1) « J) 

Tf  J,  1  >= CMJETI  J*l>  *CP0UT* TDATUM+MJETC  J)*C;>IN*<TC  J-l,  1>-TDATUM>* 
10.24*DELM(J)*(T<J-1,2 ) -T D A TUM ) -R ADOU T ( J ) +2 . E -5  *  I M JE T ( J ) * 

2 JJETl J) **2-MJET(J+l >*UJE  Tt J  +  l ) * *2 ) ) /M JET ( J*1 ) /CPOUT 
r(J»l>=T(J,l)«-R*CT<Jtl  )  -  SAVE) 

DEV=A8S<T< J,1)-SAVE) 

IF  <DEV-DEVMAX)  140,140,135 
135  QEVMAX=DEV 
140  CONTINUE 

DO  150  J=2  ,  IND3 
3AVE=TC J.2) 

MJ,2>=(a.24»<<DrLM(J*l>*MSEC(J*l>>  * T D A T UM* MSEC < J > *  1 T C J-l » 2 ) 
l-TDATUM)  )*AREAR*HCONVCJ)  *T(  J,  3))  7(0.24*  (DELMt J*1 >  ♦MSEC  <  J«- 1  >  ) 
2*AREAR*HC0NV(J>  > 

T<J,2)  =  T<J,2)«-R*<TfJ,2)-SAVE) 

DEV=ABS(  T (J»2)-SAV£ ) 

I"  CDEV-DEVMAX)  150,150*145 
145  DEVMAX=OEV 
150  CONTINUE 

DO  16  5  J= 2, I  NO  3 
3AVE=TC J, 5) 


RADINf J>=0. 

DO  155  11=1, N 

:mss=em  iss<  r a  i ,1 > ,l< ii>  ,  ii> 

ALFA=EMSS*SQRT<TCII ,1 >/TtJ,3) > 

RADIN<J>=RADIN<J>*SI3MA*ARA0<II>*F<II,J)*(ZMSS*T(II,1>**4-ALFA 
1*T< J, 3>  **4>*P7*F1< 1 1, J>  *  <1 .-EMS3>*f  T< II ,3) **4-T« J, 3 1**41 

155  CONTINUE 

T f J*3>=<P53( J> *  T(J,4)*P3*<T<J*1,3>«-T(J-1,3>>  *ARE  A  R  *  HC  ON V {  J >  * T <  J  , 2  > 
l*RAOINt J> )/<P5B<J)*2.*P3+AREAR*HCONV< J) ) 

T <J,3>  =  Tf J,3»*R*<  TCJ,3>-SAVE> 

DEV=A8S<T< J,3)-3AVE> 

IF  (OEV-DEVMAX)  165,165,160 
160  DEVMAX=OZV 
165  CONTINUE 

C  THE  DOWNSTREAM  CALCULATIONS  DEPEND  UPON  WHETHER  THE  JET  EVER  REACHES 
C  AUGMENTER  WALL 

l-  <  NSEP-N  )  170  ,230,260 

C  THE  CASE  WHERE  THE  JET  REACHES  THE  WALL  UPSTREAM  FROM  AJGMENTER  EXIT 
170  DO  185  J=NSZP,KK 

RAOOUTC J»=SISMA*ARADt J)*EMISStT< J,l) ,L(J),J)*(T(J,1)**4 
1-3QRT (T  CJ,l )/T (J, 3) )  *  T  <J«3)**4  ) 

CP  I N=SPHT ( T<J-1  ,1 )  ,J-1  ) 

C3OUT=SPHT ( T( J, 1) , J) 

3AVE-TC J,1 » 

T< J,1 )  =  CMJET( J*l>  *C?OUT*TDATUM+AREAR*HC3NV( J) *T< J , 5  >*MJET  < J) 
l*CPIN*(TfJ-l,l)-TDATJM)*0»24*DtLM<J)*<T(J-l*2  ) -TD  AT  UM  ) -R  A  DO’J  T  <  J  ) 
2*2.E-5*<MJZT< J) *UJET( J) *  * 2-M JE T { J* 1 ) *U JE T  C J *  1 ) *  * 2 ) ) ✓(  M JET  C J* 1 )  * 

3  2  POUT  *AREAR*HCONV< J) ) 

T ( J  ,  1 >  =  T<J»1)«-R*<T<J,1)-3AVE) 

DEV=ABS(T  <J, ll-SAVE) 

Ir  OEV-DEVMAX)  195,185,180 
1B0  DE V MA  X  =  OE  V 
185  CONTINUE 

DO  200  J=  NSEP , KK 
RADINf J)=0  . 

DO  190  11=1, N 

ZMSS=EMISS(T  C  II  ,1  )  *  L (  I  I  )  ,  I  I  ) 

ALFA  =  EMSS*SQRT(T< II ,1 )/T(J,3>  ) 

R AO  INC J)=RADINC J) *SIGMA* ARADf II ) *F< 1 1  * J) *(EMSS*T( II  ,1 )**4 
1-ALFA*TCJ,3>**4)*P7*F1CI I,J»*  C 1 . -EMSS > * < T ( I I , 3 > **4-TC J,3>**4> 

190  CONTINUE 

SAVE=T<J,3) 

TC J,3) =<P5BC J) *TC J,4) *P3* (T< J*1 ,3)*T( J-l , 3) > ♦ AR E A R * HC ON V C J ) * T ( J  ,  1 1 
l *RAOI N<  J) > /<P53f J  >  *2. *P3*ARrA  R*HC0NV ( J) ) 

Tf J,3)=T< J,3)*R*( T< J, 3) -SAVE) 

DEV  =  ABS<T (  J  ,  3  ) -SAVE  ) 

IF  <DE V-DE  VMAX I  200  ,200,195 
195  DEVMAX=OEV 
200  CONTINUE 

3 AVE  =  T  C  N, 3 ) 

RADINfN>=0. 

DO  205  11=1, N 

EMSS=EMISS<T< 1 1 ,1 )  ,L<  1 1)  ,  1 1) 

A LF  A  =  EMSS*SQR  T ( TC 1 1 , 1 ) /T  < N  »  3) > 

RADIN(N)=RAQIN(N) *3  I G  M A  *  ARADf  II) *F( 1 1 , N ) * f EMSS * tf II  ,1 >**4-ALFA 
l*T  t  N»  3 ) *  *4 ) *P7 *F 1 f 1 1 , N ) *  f 1  •  -EMSS  >  *  f  T  t 1 1 » 3 >  *  *4-T  tN  , 3  > • *4  ) 

205  CONTINUE 
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TCN,?>=tPlB(N)*TtN,4)  ♦P3*T(IM-1  ,3)«-KAFXT*TEXT«-'32MC0NV(N>*T<N,l> 
1«-RA0INCN))/<P13<\)  ♦P3+KAEXT*P2*HC0NV(N)  > 
T(N.3>=T(N,3>«-R*<T(N,3)-SAVE> 

3EV=ABS<T{N.5)-SAVE> 

IF  (DEV-3EVMAX)  215.215*21!) 

210  3  E  V  MA  X  =  DE  V 
215  CONTINUE 

SAVE  =  MN.1> 

C30UT=SPHTCT<N,1) »N) 

:pin=sphkt<n-i »i ) .N-i) 

EMSS  =  EMISS(T(N.1>  ,L<N>»N) 

RADOUT<N)  =  SIGMA*ARAC)<N>*EMSS* (T{N»t)**4-SQRTCT(N»l>/T(N,3>> 
1*T<N,3>**4> 

T  (N.1  )  =  (M  JE  T  <  N-*- 1 )  *CPOUT  *  TOATUM-RAOOUT  ( N) +M JE  T  (N)*CPIN*  CTTN-1  .1  > 
1-TDATUM)«-P2*HC0NV<N)*T<N  , 3 > +2 . E -5  * < M JE T < N >  *J  JET < N > * *2 -M JET C N*1 > 

2  *  J  JET ( N*1 )**2»»/C MJET C N* 1 > *CP OU T *P2 *HCONV ( N > ) 

TfN.l)  =  TCN.n*R*<T(N»l)-SAVE) 

3EV=A8S<T(N,1)-SAVE) 

IF  (0EV-OEVMAX)  225.225.220 
220  3EVMAX=DEV 
225'  CONTINUE 

30  226  J=NSEP»N 
TC J,2)=TCJ,l) 

226  CONTINUE 
30  TO  305 

THE  CASE  WHERE  THE  JET  REACHES  THE  WALL  AT  THE  AUGMENTED  EXIT 
230  3 AVE=T ( N. 3) 

RADIN<N)=0. 

30  235  11=1, N 

■MSS=EMISSCT(II*1 >,LCII>,II) 

ALFA=EMSS*SQRTCT< II ,1 )/T (N.3) ) 

RAD  ININ) =RAOININ)+SIGMA* ARAOC  II) *FC II ,N) *t£MSS*Tl II  ,1 )*»4-ALFA 
l*TtN,3)**4)*P7*FlIII,N)*tl.-EMSS>*ITtII,3)**4-TlN,3)**4> 

235  CONTINUE 

TIN,3)=IP131N)  *TIN,4)«-P3*TIN-1 , 3 > +KA E XT  *  TEXT *P2 *HCONV  IN)*TIN,1 ) 

1 *R ADI Nl N) )  /IP1B  IN  ) +P  3  ♦KA  E  X  T  *P  2  *HC  ON  V  l  N  )  ) 

T  CN»3)  =  T(N,5)*R*<T(N, 31-SAVE) 

3EV=ABSITIN,3)-SAVE ) 

IF  IOEV-DEVMAX)  24 5,245.240 
240  3EVMAX=0EV 
245  CONTINUE 

3 A  VE - T (N.l) 

CD0UT=SDHT(TIN,1) ,N) 

C  3  I N=SP  HT I T I N- l .1 ) »N-1 ) 

:mss=emissct<n,i)  .linj.N) 

RAD0UTCN)=SIG^A*ARAD(N)*EMSS*CTCN,1)**4-SQRTCT<N,1)/TCN.3)> 

1*TIN,3)**4) 

TCN.l  )  =  IMJETIN«-1>  *CPOUT*  TDATUM*MJET  (N>  *CPIN*  (TC  N-l,  1 ) - T D AT UM ) ♦ a  2 
l*HCONVIN)*TIN,3)*0.24*DELMtN) *  I  TIN-1 , 2 ) -T D AT UM ) -R ADOU T f N) 
2«-2.E-5*IHJETIN)*UJETIN)*  *  2-MJE  T I  N  «■  1 )  »U JE T  IN* 1 ) * *2 ) ) /I M JET  I N* 1 ) 
3*CP0UT*P2*HC0NVIN) ) 

TIN, l)  =  TIN,l)*iR*ITIN,l)- SAVE) 

3EV=ABStT (N.l )-SAVE ) 

IF  IDEV-OEVMAX)  255,255,250 
250  QE V MAX=  OE V 
255  CONTINUE 

TCN,2)=T<N,1 ) 


30  TO  305 

c  th-:  cas e  where:  the  jet  never  reaches  the  augmenter  wall 

260  S A VE=T  t  N,  3  ) 

RADIN<N)=0. 

30  265  11  =  1. N 

EHSS  =  EMISSCT  CII *1 ),L<  II) »II) 

4LFA=EMSS*SQRT(T<  II  .1  1/MN.3)  ) 

R AOIN(N)=RADIN<N) ♦SIGMA* ARADt  II) *Ff 1 1 . N) *  I EMSS* T { 1 1  ,1 ) ** 4- ALFA 
l*T«N,3) **4)*P7*Fl ( IItN) *  C1.-EMSS)*CT  C II .3) **4-T<N,3)**4  ) 

265  CONTINUE 

TIN,  3)=<P1BCN)*  TC  N.  4)  *P3  *TCN-1 .3)  *KAEXT*TE  XT  *P2 * HCONV C N ) * T t N ,  2) 
1*RADIN<N) )/(P18(N)*P3*KAEXT*P2*HC0NV(N) ) 

TCN,3)=T(N,3)*R*I T<  N. 3 ) -  SAVE ) 

DEV=ABS<TCN,3)-SAVE> 

IF  tDEV-OEVMAX)  275,275.270 
270  OEVMAX=DEV 
275  CONTINUE 

SAVE=T(N*1) 

RAD OUT  C  N)  =  Q« 

EMSS=EMIS3(T<N,1)  »L(N  )  ,N  ) 

30  280  11=1, N 

ALFA  =  EMSS*SQRT(T(N,1)/T(  II  ,3)  ) 

RAD0UT<N)=RAD0UT<N)*SIGMA*ARAD(N)*F<N,II)*<EMSS*T(N»1)**4-ALFA 
1*T< II ,3>**4) 

280  CONTINUE 

CPIN=SPHT(TCN-1,1 ),N-1) 

C30UT  =  SPHT(T(N,1)  ,N) 

T«N,1 )  =  CMJETC N+l) *CPOUT*TDATUM*MJETCN) *CPIN*CT{ 1-1, 1)-T DATUM)* 
13.24*DELM(N) *<T(N-1,2 ) -T DA TUM ) -R ADOU T C N > *2 . E -5  * t M JE TC N > * 
2JJETCN) **2-MJET(N*l )*UJ£TCN*1 >  *  *  2 >  > / M JE T < N* 1 ) /CPOUT 
TCN»1)=T(N,1)*R*(T<N, 11-SAVE) 

OEV=ABS(T(N,l)-SAVE> 

IF  (DEV-DEVMAX)  290,290.285 
285  OE VMAX=DE V 
290  CONTINUE 

SAVE=T ( N, 2 ) 

TIN,?) - <0.24* MSEC  <N*1 > • T DATUM* MS EC < N ) * l T < N-l , 2 ) -T DA T J M ) *P 2 
1* HCONV <N)*T<N,3))/<0. 24 *MSEC<N*1)*P2*HC0NVCN)) 

T (N»2)=TCN,2)*R*CT  CN,2 ) - S AVE ) 

3EV=ABS<T (N, 2)- SAVE) 

IF  (QEV-OEVMAX)  300,300,295 
295  DEVMAX=DEV 
300  CONTINUE 
305  CONTINUE 
C 

C  CHECK  TO  SEE  IF  ITERATION  IS  CONVERGING 
IF  (DEVMAX.LE.DSAVE)  GO  TO  310 
WRITE <6,36 ) 

WRITEI6.49) 

WRITEC6.33) 

30  TO  500 

C  CHECK  TO  SEE  IF  ITERATION  HAS  CONVERGED  TO  WITHIN  TOLERANCE 
310  IF  (DEVMAX.LE.TOL)  GO  TO  405 
C 

400  CONTINUE 
C 

WRITEC6.36) 
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WRITE(6»51>  .'OL.ITER 
WRITEI6.38) 

405  CONTINUE 

WRITE<6,89)  I *  DEV  M A  X 

C  CONI/  ERT  TEMaS  FROM  RANKINE  TO  FAHRENHEIT  AND  aR  INT  RESUlTS 
DO  410  M=1»N 
DO  410  J=l«4 
TEMP(M*J)=T (M»J)-460. 

410  CONTINUE 
11  =  1 
I SEC=0 

DO  415  M=  l  »  5 
lIM=10*M 

IPRINT=MINOILIM*N) 

IaRNT=MINOILIM.NSEP-l ) 

WRITEC6.52) 

WRITEC6.44)  <ZCJJ » J  =  1 1 . I  PR  I NT ) 

WR I TE I  6  *53 )  (TEMPI J,1 ) *J=II .IPRINT) 

IF  (ISEC.GT.Q)  GO  TO  412 

WRITE  I  6 *54)  (TEMP(J*?).J=II.IPRNT> 

412'  WR I TE I  6  *56  )  (  TEMP  I  d  ,  3  >  »  J  =  1 1  ,  I  PR  I  NT  > 

WRITE  I  6  *57)  (TEM?(J*4)*J=II»IPRINT) 

IF  (IPRINT.GE.N)  GO  TO  420 
11=  IPR INT  +  1 

IF  (IPRNT.GE.JX  )  ISIC=1 
415  CONTINUE 

420  CONTINUE 

C  PRINT  JET  VELOCITT  AT  AUGMENTER  EXIT 
iWR  I  TE 1 6  *69)  UJETINI 

C  PRINT  THE  HEAT  TRANSFER  CHARACTERISTICS  OF  THE  JET  EXHAJST 
DO  421  M  =  1,N 

EMIS(M)=£MISS(T(M,1 ) »L(M  )*M) 

421  CONTINUE 
11  =  1 

DO  422  M=l»5 
L I M=1 0*M 

I?RINT=MINO(LIM,N) 

WRI TE 1 6  *72 ) 

WRITE  16*44)  (ZC J) »J  =  I I  * IPRINT > 

WRITE  I  6*73)  (HCONVC J) * J= 1 1  * IPRINT ) 

WRITE  16*74)  (EMISI J) » J= II  .IPRINT  ) 

IF  (IPRINT.GE.N)  GO  TO  423 
1 1  =  IPR I NT  +  1 

422  CONTINUE 

423  CONTINUE 
C 

C  CALCULATE  THE  STEAM  GENE R AT  I  ON. . • . A SSUME  FEED  HATER  ENTERS  THE  BOILER 
C  IN  A  SATURATED  CONDITION 

IF  ITYPE.EQ.l)  GO  TO  500 
IF  CTYPE.EQ.2)  WRITEI6.28) 

I~  ITYPE.E0.3)  WRITE(6»29) 

Ir  (TYPE.E0.4)  URITEC6*31) 

IF  ITYPE.EQ.3)  GO  TO  455 
C  STEAM  GENERATION  IN  WATER  ’WALL  BOILER 
DO  424  I =  1  * N 
3CONVII )=0. 

MSTM1 C I )=0 • 


* ADST (  I  )  =  0. 

424  CONTINUE 

A30IL(l)=A3QILtl>/< I 3 ACK - I F RONT ) 

MSTM1 <  I FRONT)  =  ABO  IL(1  )/2  .*P43  C  IFRONT  >  *(  T<  IFRONT  ,3  )-T(  IFRONT, 4)  ) 
1/HSTMC1 ) 

MSTM1(IBACK)=AB0IL( 1 )  /  2  •  »P43 1 IBACK)*(T( IBACK»3)-T  (I  BACK, 4)) 
l/HSTMt 1 ) 

STEAMlr  M3TMK  I FRONT > ♦ M ST  Ml <  I  BACK) 

J-RONT=IFRONT*l 

J3 ACK= I BACK-1 

30  425  I=JFRONT  ,  JBACK 

M3TM1CI)=AB0ILC1>*P43(I)*<TCI,3>-T(I,4>>/H3TM<1> 
STEAM1=STEAM1+MSTM1<I  ) 

425  CONTINUE 
rSTM(l>=TSTMtl)-460. 

A  30  IL ( 1 )=A8QIL( 1) *  ( I3ACK-IFRONT) 

3T0T=STEAM1*HSTM< 1 ) 

JR  I TE (6 , 32  )  TSTM(1)»PSTM(1),HSTM<1)*AB0IL<1> 

JRITEC6.77)  FRONT  ,  8  AC  K 
JRITE<6,53>  STEAM1 
JR  I TE  <  6  »  76  )  QTOT 
J  R I TE  (  6  ,  59 ) 

11  =  1 

30  430  1=1,5 
LIM=10*I 

IPR INT=MINO<LIM,N) 

JRIT~(6«44)(Z(J)«U=II »IPRINT> 

JRITEC6,61XMSTMl<JI,  J=  1 1,  IPR  I  NT) 

IF  C  IPRINT.GE.N)  GO  TO  435 
II=IPRINT*1 
430  CONTINUE 
435  CONTINUE 

C  CALCULATE  AND  PRINT  FRACTION  0-  STEAM  GENERATION  DUE  TO  CONVECTION 
C  AND  TO  RADIATION 

3C0NVI IFRONT >=HC0NV( IFRONT >*P2*(Tt IFRONT, 2 >-T( I  FRONT, 3) ) 
RADST(IFRONT)=RADINf IFR0NT3/2. 

3TOT  =  ABS(QCONVCIFRONT)  > ABS  C R  ADS T  <  I FR ONT  )  > 

3C0NVI I FRONT  I =QCONV II  FRONT) /QTOT 
RAOSTI  IFRONT  XRAOSTH  FRONT) /QTOT 

3C0NVI IBACK)=HCONV< 13 ACK  I  *P2*  <  T C IBACK,2>-T(  I  BACK,  3)  ) 

RADST! IBACK)=RADIN( I  3 ACK  )/2. 

3T0T=ABS(QC0NV( IRACK) ) *A BS < R A DST ( 13 AC K ) ) 

3C0NVI I3ACK)=QC0NV< IBACK  3/QT3T 
RADST (IBACK)=RADST< IBACK) /QTOT 
30  440  I=JFRONT, JBACK 

3C0NV ( I )=HCONV  < I) * ARE AR * C T ( I , 2 ) - T ( I , 3 ) ) 

RADST  <I)=RADIN<I) 

3T0T=ABS<  QCONV (  I )  )  ♦■ABS  (R  ADST  (  I  >  ) 

RADST (  I>  =  RADST<I)/QTOT 
QCONV < I )=QCONV(I)/QTOT 
RAO  INI  I)=RAOINI I) /QTOT 
440  CONTINUE 
11  =  1 

30  445  1=1,5 
L I M  =  1 0  *  I 

13R INT=MINO<LIM,N) 

WRITE(6,62> 
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i 


n  o 


.'RITF{£,,44>  (ZCJ) *J=:i,IPRINT> 

WRITE (6 *63)  OCONIV(Jl*J=II«r?RINT> 

WR I  TEC  6  *64)  <RADST<J),J=II,IPRINT> 

I-  C  I  PR  INT  .GE  •  N  )  GO  TO  *50 
1 1  =  I°R I  NT ♦ 1 
445  CONTINUE 
450  CONTINUE 

IF  (TYPE.EQ.2)  30  TO  460 

C  STEAM  GENERATION  IN  CONVECTION  BCILER.. ..ASSUME  FLOW  OVER  STAGGERED 
C  WATER  TUBES. ...PROPERTIES  WILL  BE  BASED  UPON  FILM  CONDITIONS 
4  55  CPAVG  =  (  0.24*MSEC(  N>  +CPOU  DMJETCN)  )  /(  MSEC  <  N )  «•  M  JE  T  (  N  )  ) 

T  AVG=TDATUM«-C0.24  *M3EC(N)*CT(  N  *  2  ) -T  D  AT  UM  >  ♦C3  OU  T  *M  JE  T  (  N  )  *  C  T  C  N  *  1  ) 
l-TDATUM)  >/(C3AVG*  (MSEC  (N  >«-KJET  CN> ) ) 

T  G AS=T  A VG 

TAVG=0.5*(TAV5  +  TSTM€2)  > 

C  USE  THE  EUCKEN  EQUATIONS  TO  ESTIMATE  CONDUCTIVITY  AND  VISCOSITY 
rACT0R=TAVG**1.5/C225.+TAVG> 

<AIR=9.26E-4*FACT0R 

MUAIR=2.694E-3*FACT0R 

RH0AIR=PRES3/RAIR/TAVG 

VCONV  =  CMSEC(N)-*-MJET(N ) ) / AAUG/RHO A  I R 

HTUBES=0.2R7*KAIR/DTUBE*  t  R  HOA  I  R  *  V  C  ONV  *  DT'JBE/ MU  A IR  )*  *0.6 
C  IT  HAS  BEEN  ASSUMED  THAT  THE  JET  AND  SECONDARY  FLOW  ARE  COMPLETELY 
C  MIXED  BY  THE  TIME  THEY  REACH  THE  CONVECTION  BOILER 
3STM2  =  HTUBES*ABOI LC2)  * C  T GAS-T S T M ( 2 ) ) 

STEAM2=QSTM2/HSTM(2) 
rSTM(2)=TSTM(2)-460. 

T3AS=TGAS-460. 

DTU8E=12.*DTUBE 

WRITE  16*33)  TSTMC2) *PSTM(2),HSTMC2>,A80IL(2> *DTU3E 
W  R I TE ( 6  *  66 )  STEAM2 
WRI TE (6  *71 1  TGAS*HTUBES,QSTM2 

460  CONTINUE 

DETERMINE  THE  ECONOMICS  OF  GAS  TURBINE  TEST  CELL  HEAT  RECOVERY 
IF  <  T YPE.EQ .4 )  WRITEC6.9) 

WRITE (6*79) 

461  IF  CTYPE.EQ.3)  GO  TO  475 

C  FIRST  FDR  THE  WATER-WALL  HEAT  EXCHANGER 

IF  CIIND2.EQ.1). AND. (TYPE. NE.4)>  WRITEC6*9) 

Ic  CC IN02.EQ.2I .AN0.CTYPE.EQ.4  )  )  URITEC6.9) 

IF  (  ( IND2.EQ.0).0R.(I ND2  .  GE . 7  )  >  WRITEC6»31>  PRICE1 

IF  (  ( IND2.GT .0).AND.C IN0  2.LT.7) .AND. (TYPE.EQ.4)  )  URITEC6.91) 

C  CALCULATE  THE  VALUE  OF  THE  STEAM  GENERATED 
VALUE 1=0.365E-3*STEAM1*HSTM( 1 ) *T I  ME* COST 
C  CALCULATE  THE  PRESENT  VALUE  OF  THE  SAVINGS  ACCUMULATED  OVER  THE 
C  ECONOMIC  LIFE  OF  THE  HEAT  RECOVERY  HEAT  EXCHANGER 
3VC 1>=VALUE1*RATI 0 
IF  (PVC1I.LT. 0.1  PVC1)=0. 

SIR(1)=PV(1»/PRICE1 
30  465  1=2*25 

:»V(I)=PV(I-1)*VALUE1*  RAT  I  0**1 
IF  (PV(I).LT. 0.)  PVCII=0. 

3IRCI»=PVC IJ/PRICE1 
465  CONTINUE 
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C(IND2.EQ.3).0R.(IND2.GE.7>>  WRITECS.82)  0 1 SC * C OST « E SC 

WRITE(S*R3> 

DO  470  1=1*5 
J  =  5  *  I 

WRITECS*R4>  J,TIME*PV(J) *SIRC J> 

470  CONTINUE 

475  IF  (TYPE. £0,2)  GO  TO  495 
THEM  FOR  THE  CONVECTION  HEAT  EXCHANGER 

1=  (C IND?.E3.3>.AN0.C TYPE.NE.4) )  WRITE(S,9> 

Ir  C(IN02.EQ.4>.AN0.(TYPE.E3.4>>  WRITEC6,9> 

IF  (CINO?.£-3.0).OR.CIND2.GE.7))  WRITEC6.38)  PRICE2 

IF  <( IND2.GT.0) .AND.(IN02.LT.7).AN0.(TYPE.E3.4))  WRITE(6*92> 

AS  490VE*  CALCULATE  THE  VALUE  OF  THE  STEAM  GENERATED 
VALUE1  =  0.3  65E-3*^TEAM2*HSTM(2  >*T I  ME* COST 
3  V  ( 1)=VALU£1 *RATI 0 
IF  CPVC 1) .LT.O. )  PV(1)=0. 

SrRCl)=PV(l)/PRICE2 
DO  490  1=2*25 

PVC I) =P VC  I -1 1 ♦  VALUE  1*  RAT  I  0**1 
IF  C?V( I).LT.O.  )  PVCI)=0. 

SIRCt  )  =  PVC  D/PRICE2 
480  CONTINUE 

IF  CC IND2.EQ.0) .OR. Cl ND2.GE.7) >  WRITECS*87)  01  SC* COST *ESC 
WRITE (6.83) 

DO  485  1=1*5 
J  =  5  *  I 

WRITE(S*94)  JfTIME.PVCJ) *SIRCJ) 

485  CONTINUE 
495  TIME=2.*TIME 

IF  (IN02.E0.3)  TIME=3. 

Ic  (IND2.EQ.4)  TIM£=5. 

IF  (IN02.EQ.5)  TIME=8. 

I  NO 2  =  I N02 ♦ 1 

IF  CIND2.LT. 7)  GO  TO  461 
500  STOP 

1  r  OR  MA  T  (  15  > 

2  rORMAT  C  3F10.0*2F1 0.3) 

3  r  OR  MA  T  C  2F 1  0 . 1 1 

4  r  OR  MA  T  (  4F 1 0  •  l  ) 

6  FORMATC6F10.2I 

7  r  OR  MA  I  (  2C 1  0  •  3  ) 

8  r0RMATC4F10.1,c10.3) 

9  -0RMATC1H1) 

11  "ORMATC 1H1 *5X*"AUGMENTER  GEOMETRY") 

12  "ORMATC  IHO»10X*"INITIAL  AREA  OF  JET  =  ",F5.2*"  S3F T " « /I IX  * " CRO S3  S 
1ECTI0NAL  AREA  OF  AUGMENTER  =  "*F5.1»"  SOF T " * /I 1 X , "LENG TH  OF  AUGMEN 
2TER  =  "  ,F5  •  1 ♦ "  FT",/1 1 X  * " THIC KNESS  OF  AUGMENTER  WALLS  =  ",F5.3*"  I 
3N"*/11X, "THICKNESS  OF  ACOUSTIC  PILLOWS  =  ",F5.2*"  IN"//) 

13  r  OR  MA  T  C 1H  * 3X * "PHYS IC AL  CHARACTERISTICS  OF  JET  EXHAUST") 

14  rORMAT(lH0*10X*"MASS  FLOW  RATE  OF  JET  =  ",F4.0."  LB/SEC "♦ /I 1 X *" TOT 
l AL  TEMPERATURE  =  ",F5.0»"  DEGF" * /11X *"TOT AL  PRESSURE  =  "*F4.1*"  »S 
21 A"«/ll X, "RATIO  OF  SPECIFIC  HEATS  =  "*F4.2) 

15  r  OR  MA  T  C 1H  *10X, "COMPOSITION  OF  EXHAUST  (PERCENT  WEIGHT)") 

17  rORMATC  1H  »15X»"QXYGEN . "»P5.2*/16X*"NITR0GEN . 

1.. ."*F5.2»/16X* "CARSON  0  I  OX  IDE. . . . . . " » F5 . 2* /1 6X *• CARR ON  MONOXIDE.. 

a..."»F5.2*/l6X,"WATER  VAPOR . " »F5 . 2  * /I 5 X  ♦" HY DR OGEN . 

3.. ."*F5.2*/11X*"LENGTH  OF  LUMINOUS  JET  =  ",F5.l*"  FT"//) 


13 

r  DR  MA  T ( 

1H 

,  5  X  ,  " 

PH 

19 

r  OR  MA  T ( 

1H0 

,10X, 

"S 

1 

C  TEMPE 

RATURE  = 

•i 

21 

r  ORMA  T ( 

1H 

,5X," 

HE 

22 

r OR  MA  T ( 

1  HO 

,10X  , 

"T 

1 

3TU/HR 

-c  T 

-DEGF 

") 

23 

r  OR  MA  T ( 

1H 

,10X, 

"T 

Y3ICAL  CHARACTER  IS  T I C3  OF  3 
ECONDAR Y  FLOW  RATE  =  ",r5.0, 

»  F  4.  0  »  "  0E5F “//  ) 

AT  TRANSFER  PROPERTIES") 

HERMAL  CONDUCTIVITY  OF  AUGMENTER 


ECONOARY  FLOW") 

"  LB/SEC " , /I  IX,  "STAT  I 


IERMAL  CONDUCTIVITY 

1"  8TU/HR-FT-0EGF" ) 

24  r  OR  M4 T ( 1H  *  1  OX , "EFFECT  IV E 


WALLS  =  "»F4.1»" 
OF  ACOUSTIC  PILLOWS  =  ",F5.3, 


HUSH  HOUSE  CHARACTERISTICS") 


13X, "CONDUCTANCE  AT  AUGMENTER  INLET  = 


»/l 
»  "  T 
AT 


3  AUGMENTER  EXIT  = 

27  rORMAT(lH  ,5X,"TH£ 

28  FORHATC 1H1,5X»"THE 
1  THE  AUGMENTER"//) 

29  ;  (JR  MA  T(  1H1»5X,"THE 
HE  GAS  PATH  AT  THE 

31  r  OR  MA  T( 1H1*5X*" THE 

1  ROUND  THE  AUGMENTER  AND*,/6X,"A  STAGGERED 


"♦E9.3."  BTU/HR-DEGF",/16> 

»  r  >,  •  u  *  "  OE3F",/13X,  "CONDUCTANCi 
3TU/HR-DEGF"//) 

FACILITY  HAS  NO  BOILERS"/) 

FACILITY  HAS  A  WATER-WALL  BOILER  WRAPPED  AROUND 


WATER-TUBE  BOILER  IN  T 


FACILITY  HAS  A  STAGGERED 
AUGMENTER  EXIT") 

FACILITY  HAS  BOTH  A  WATER-WALL  BOILER  WRAPP 
«./cv.H4  ^Tir.c'Brn  WATER-TUBE  BOILER  IN 


2  GAS  PATH  AT  THE  AUGMENTER  EXIT"//) 

32  r  OR  MA  T ( 1H  ,/6X, "WATER-WALL  BOILER  CH AR ACT ER I ST  I CS " » /I  6 X  * "STE AM 
HERATURE  =  ",F4.0,"  DE  GF  "  * /1 5  X  *  "  S  TE  A  M  PRESSURE  =  ",F5.1,"  °SIA 
2SX,"HEAT  OF  VAPORIZATION  =  ", F6.lt"  B TU/LB" » /1 6 X , "T OT A L  HEAT  T 


ED  A 

the 

TEM 
"  »/l 

BANS 


33  rORMAT( 1H0»//6X, "CONVECTION  WATER-TUBE  BOILER  CHARACTERISTICS" 
1 X,  "STEAM  TEMPERATURE  =  ",F4.0t"  DE3F » * /I 6 X  ,  " S TE AM  PRESSURE  =  " 
2 1 » "  PSIA" ,/16X,"HEAT  OF  VAPORIZATION  =  ", F6.lt"  BTU /LB " t / 16X t " 


3L  HEAT  TRANSFER  AREA  = 

43  =  ",F4.2t"  IN"/) 


»F5.0  t "  SOFT", /16X, “OUTER  DIAMETER  OF 


34  -ORMATC1HO 
1ITY  =  ",F5 
2/SEC" » //6X 

36  r  OR  MA  T ( 1H0 
!•**““*** 

37  rCRMAT(lH 
1ALYSIS  IS 

38  r  OR  MA  T ( 1H 
1*“*»»**** 

39  rORMAT C 1H 

I  TO  WITHIN 

41  -ORMATtlHO 
1DISTANCE  OF 

42  ; 3R  MA  T ( 1H0 

43  r  OR  MA  T  < 1H0 
1 T I  ON  (FT) 


INLET", /HXt"UET  V 
VELOCITY  =  ",F4.0, 


5X , "BULK  VELOCITIES  AT  AUGMENTER 
0,"  FT/SEC",/11X," SECONDARY  FLOW 
"CRAYA-CURTET  NUMBER  =  ",F4.2/) 

///GX, "******************************************** 


t/16 

,F5. 

TOTA 

TUBE 

ELOC 

"  FT 


*  **  * 

***.»»**•»••***»*. .•••ft**.************,'/) 

FOR  JET  SEPARATION  POINT  IS  DIVERGING,  AN 


5X, "ITERAT ION 
ERMINATED"/) 
5X, "*****»***- 


FOR  JET 
AFTER  " 


5X, "ITERATION 
"  ,F  3 .1 , "  FEET 
5  X  ,  " JE  T  EXHAUST  REACHES 
",F4.0,"  FEET  FROM  THE 

N 


SEP ARAT I  ON 

■ 


,13,"  HERAT 
THE  AUGME 
INLET"/) 


44  -ORMATdH 
46  r  DRMAT( 1H 


/6X  ,  "MASS  FLOW  RATES  (LB/SEC ) 
AND  JET  RADIUS  (FT)"/) 

5X,"DIST  FM  INLET", IX, 10F10.2) 
5X, "RADIUS  OF  JE T "  ,  1 X , 1  OF  1 0 .2 ) 


TEMPERATURES 


VAR 


POINT  DID  NOT  CONV 
IONS"/) 

NT  ER  WALL  AT  AN  AX 

AUGMENTER  WALLS"/ 
I  AT  ION  WITH  AXIAL 


) 

ERSE 

I  AL 
) 

LOCA 


47 

-OR 

MA 

T  (  1H 

»5X , " 

JE 

T  FL 

OW" 

,  GX  ,1 

48 

r  OR 

MA 

T (  IH 

,5X," 

SE 

CONO 

ARY 

FLOW 

49 

r  OR 

MAT( 1H 

,  5X  ,  " 

IT 

ERAT 

ION 

FOR 

1 

S  TER 

MINAT 

ED"/) 

51 

FOR 

MAT ( 1H 

,5X," 

IT 

ERAT 

ION 

FOR 

1 

IN 

"* 

F  3 . 1  » 

"DEGF 

AFTER 

", 

13," 

52 

r  DR 

MA 

T  (  1H0 

,  /  SX  , 

"T 

EMPE 

RATURE  ( 

l 

-T) 

• 

the 

WALL 

00 

INCLUDES  AC 

53 

r  OR 

MA 

T(  1H 

|5X," 

JE 

T  TE 

MP" 

» GX  » 1 

TEMPERATURES 

ITERATIONS"/) 


S 

IS 

DIVE 

R5 

ING, 

ANA 

LYS 

IS 

I 

S 

DID  NOT 

C 

ON  VE 

RGE 

TO 

WITH 

/) 

T  I 

3N 

WITH 

A 

XIAL 

LOC 

ATI 

ON 

( 

ous, 

IF  A 

NY 

."/) 
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54  -ORMATdH  ,5X,"SEC  FLOW  TEM°" *  1 X , 1 0" 1 0 . 1 > 

56  r  OR  MA  T  < 1H  ,5X,"AUG  WALL  ID  "  .  3  X  ,  1  OF  1  0 . 1  > 

57  r  OR  MA  T (  1H  ,  5  X  ,  "  A  U  G  WALL  OD ” ♦ 3 X » 1  OF 1 0 . 1  ) 

58  r  OR  MA  T  d  H  ,/6X, "TOTAL  STEAM  GENERATION  IN  WATER-WALL  BOILER  -  ",E1 
10.4."  LB/HR "/ ) 

59  rORMAT(lHO«5X," WATER-WALL  STEAM  GENERATION  CL3/HR)  VARIATION  WITH 
1AXIAL  LOCATION  (FT  PROM  AUGMENTER  INLET)"/) 

61  r ORMA T < 1H  .5X» "STEAM  GEN " « 5X . 1 0" 1 0 . 1  ) 

62  -  OR  MA  T ( 1H  ,/6X . "FRACT ION  OF  WATER-WALL  STEAM  GENERATION  BY  CONVECT 
1 1  ON  AND  BY  RADIATION"/) 

63  -ORMATdH  . 5X, "CONVECTION", 4X  .13F10. 3) 

64  r  OR  MA  T (1H  ,  5X  ,  "R  A  0 1  AT  I  ON  "  ,  5X  ,  1  OF  1  0 . 3  ) 

66  r OR MAT ( 1H  ,  /6X, "TOTAL  STEAM  GENERATION  IN  CONVECTION  WATER-TUBE  B 
10ILER  =  ",£10.4,"  LB/HR") 

67  rORMAT( 1H  ,5X, "AMBIENT  CONDITIONS") 

63  rORMAT( 1H0 , 10X, "TEMPERATURE  =  ",F4.0,"  DEGr /I IX  ."CONVECT IVE  FILM 

I  COEFFICIENT  =  ",F5.1,"  B TU/H R-SQF T-OEGF "// ) 

69  rORMATdHO  ,/6X.  "BULK  VELOCITY  OF  JET  AT  AUGMENTER  EXIT  =  ",F5.0," 
lrT/SEC"/> 

71  -ORMAT( 1H0.5X, "AVERAGE  GAS  TEMPERATURE  AT  BOILER  INLET  =  ",F5.0," 
1DESF" ,/6X, "BOILER  OVERALL  HEAT  TRANSFER  COEFFICIENT  =  ",F5.1,"  BTU 
2/HR-SQFT-OEGF"  ,/6X, "TOTAL  HEAT  TRANSFERRED  =  ".E10.4,"  BTU/rfR") 

72  -ORMATC1HO./6X, "EXHAUST  GAS  EMISSIVITY  AND  CONVECTIVE  FILM  COEFFIC 

I I  ENT  (BTU/HR-30FT-DE3F )  BY  LOCATION"/) 

73  r  DR  MA  T ( 1H  ,5X,"FILM  C OEFF" ,4X , 1  )F1 0 . 2 > 

74  rORMAT(lH  , 5X , "EMI  SSI  V  I T Y" ,4X , 1 QF 1 0 . 4  ) 

76  -ORMATdH  ,5X, "TOTAL  HEAT  TRANSFERRED  =  ",E10.4,"  BTU/HR"/) 

77  -ORMATdH  ,5X , "  THE  WATER-WALLS  ARE  LOCATED  ALONG  THE  SECTION  BEGIN 
INI NG  " , F5 • 1 , "  FEET  FROM  THE  AUGMENTER  INLET " ,/SX, "AND  ENDING  AT  A 
2DISTANCE  OF  ",F5.1,"  FEET  FROM  THE  AUGMENTER  INLET"/) 

73  FOR  MAT (F10.1 , 1 1 0 , 2F 15  .2) 

79  -ORMAT(1HO,//6X  ,"THF  ECONOMICS  OF  GAS  TURBINE  TEST  CELL  ENERGY  REC 
10VERY") 

81  -ORMAT(1HO,//6X,"TH£  PURCHASE  PRICE  OF  THE  WATER-WALL  HEAT  EXCHANG 
1ER  IS  *",F10.2,"»  PAYABLE  WHEN  INSTALLATION  IS  COMPLETE"/) 

82  rORMATClHO,5X,"THE  ACCUMULATED  PRESENT  VALUE  OF  SAVINGS  RESULTING 

1  ROM  THE  ADDITION  OF  THE  WATER-WALL  BO ILE R " , /6X , " AS SU M I NG  A  DISCOU 
2NT  RATE  OF  «,F4.1,"  PERCENT,  A  CURRENT  STEAM  PRTCE  OF  $",F5.2,«/MB 
3 T J ,  AND  A", /SX, "STEAM  ESCALATION  RATE  OF  "F4.1,"  PERCENT") 

83  rORMATClHO,6X, "ECONOMIC", 6X , " OPE R AT  I  ON " ♦ 1 1 X , "ACCOM  PV  OF",12X,"SIR 
1"»/6X,"LIFE  (YRS) ",5X, "(HRS/DAY)", 10X, "STEAM  GEN  ($>") 

84  r ORMAT( 1H  , 8X , I  2 , 1 2 X , F 4. 2 , 10X  ,F I  4 . 0 , 1  OX , F6 . 2 ) 

87  rORMAT(lHO ,5X,"THE  ACCUMULATED  PRESENT  VALJE  OF  SAVINGS  RESULTING 
1-ROM  THE  ADDITION  OF  THE  CONVECTION  B 0 1 LER " , /6 X , " AS SUM  I NG  A  DISCOU 
2NT  RATE  OF  ",F4.1,"  PERCENT,  A  CURRENT  STEAM  °RICE  0-  $",F5.2,"/M8 
3TU,  AND  A", /6X, "STEAM  ESCALATION  RATE  OF  ",F4.1,"  PERCENT") 

38  ' 0RMATC1H0,//6X,"THE  PURCHASE  PRICE  OF  THE  CONVECTION  HEAT  EXCHANG 
1ER  IS  i",F10.2»",  PAYABLE  WHEN  INSTALLATION  IS  COMPLETE"/) 

89  rORMAT(lHO,2X,"  SOLUTION  CONVERGED  AFTER  ",I3,"  ITERATIONS  TO  WITH 
1IN  ",F4.2,"  DEGREES"/) 

91  -ORMATC1HO  ,/6X, "WATER-WALL  HEAT  EXCHANGER") 

92  FORMATClHO»/6Xt"CONVECTION  (WATER  TUBE)  HEAT  EXCHANGER") 

END 
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rUNCTION  SPHTCT.I  » 

C  PROS  R A M  TO  calculate:  THE  MEAN  SPECIFIC  HEAT  BETWEEN  T  AND  SO  3EGF 
C  OF  A  MIXTURE  OF  GASES. ... THS  GASES  CONSIDERED  ARE  OXYGEN*  NITROGEN 
C  CARSON  MONOXIDE.  HYDROGEN*  WATER  VAPOR*  AND  CARBON  DIOXIDE 
C  SPECIFIC  HEAT  CALCULATED  IN  8TU/LR-DEGF 
REAL  MU *MOLES«N2*NITC 50) ,MONOX(50) 

COMMON  OX Y, NIT, 01  OX, WATER*  MONOX ,HYD» LUM 
DIMENSION  OX Y ( 5 0 ) * D IOX ( 5 0 ) * WA TER C 50 )  .HY0C50) 

DIMENSION  CP  <  6 ) 

DATA  TDATUM/520 •  / 

SPHT=0. 

DELTA=T-TDATUM 
32  =  0XYC  I) 

N2  =  NITC  I> 

CO=MONOXCI  ) 

H2  =  HYD<  I) 

HZO=WATERCI) 

C  02  =D I OXC I ) 

C  FIRST  CALCULATE  THE  MEAN  SPECIFIC  HEAT  OF  INDIVIDUAL  COMPONENTS 
C  USING  THE  RELATIONSHIPS  OF  SUEIGERT  AND  BEARDSLEY,  REF.  GEORGIA 
C  INST.  OF  TECH.  BULLETIN  2  <1938) 

CP< 1>=< 11 .515 *0ELTA-344. *(SDRTI  T > -SQR T f TO A TUM ) > +1 5 3 0. 
1*<AL0G(T)-AL0GC  TDATUM) ) >/DELT A* 02 
Eot2)=(9.47*D£LTA-3.47E3*(AL0G(T)“AL0GCT0ATUM) >-l.lSE6 
I*( I. /T- I. /TDATUM) D/OELT  A*N2 
CPC  3) =  f 9.46 *0ELTA-3. 2 9E3 * C ALOG < T > - ALOG C TDATUM) )-l .07“ 6 
1*< 1  ./T- 1. /TDATUM) )/DELTA*CO 
C°<4)=C5.76*0£LTA*2.89E-4*CT*  *  2- TO  A  TUM*  *  2  )  *-4  0  .  *(SQRTCT) 

1 -SORT C TDATUM))) /DEL TA*H2 

CP  C5)  =  <  19.86*0ELT  A- 11 94.  * C SQR T C T ) -S QR T f  TDATUM)  )«-7  50  0. 

1*( ALOG(T)-ALOG<  TO  ATOM) ))/OELT A*H20 
CPC6)=< 16.2*D£LTA-6.53E3MAL0GCT>-AL3G(TDATUM> ) -1 .41E6 
1M1./T-1./T0ATUM) ) /DELTA *C 02 
C  THE  TOTAL  NUMBER  OF  MOLES  IN  THE  MIXTURE 
M0LES=02*N2*C0>H2 ♦H20+C02 
C  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE 

MW=C02*32.*N2*28.*C0*28.  «-H2*2  .  *820*1  8.  *C02*  44.  ) /MOLES 
C  ALLOWING  THE  MEAN  SPECIFIC  HEAT  TO  3E  DETERMINED 
DO  5  J=l,6 

SPHT=SPHT*CP< J) /MOLES /MW 
5  CONTINUE 
RETURN 
END 
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cUMCTION  EMISS(T,L,I> 

C  PR03  R A  M  TO  CALCULATE  THE  EMISSIVITY  OF  A  MIXTURE  OF  COMBUSTION  GASES  • 
C  CARBON  DIOXIDE,  WATER  VAPOR,  AND  CARBON  MONOXIDE  ALL  CONTRIBUTE  ANO 
C  THEIR  INDIVIDUAL  EMISSIVITIES  ARE  OBTAINED  FROM  CURVE  FITS  OF  THE 
C  TABLES  OF  MCADAMS,  "HEAT  TRANSMISSION",  MCGRAW-HILL, (1054) 

REAL  L , MOLES *N2 ,NIT<50>  , MONOX <50 ) 

COMMON  OXY,NIT,OIOX,WATER,MONOX,HYD*LUM 
DIMENSION  OX Y ( 50 I ,DIOX<50),UATER<50)  ,HYD<50> 

DATA  OELTA/O./ 

C 

C  IF  SEGMENT  "I"  OF  THE  JET  EXHAUST  IS  LUMINOUS ,  SET  THE  EMISSIVITY 
C  EQUAL  TO  ONE 

IF  (I.GT.LUM)  GO  TO  5 
EMISS=1. 

30  TO  50 
C 

5  02=OXY<I> 

N2=NIT< II 
C  0  =  MONO  X  < I  ) 

M2=HYD<I) 

H20=WATER<I I 
C02=0 IOX( 1 1 

C  CALCULATE  THE  PARTIAL  PRESSURE  OF  INDIVIDUAL  GASES  IN  ATMOSPHERES 
C  ASSJMING  A  TOTAL  PRESSURE  OF  ONE  ATMOSPHERE 
M0LES  =  C02*H20  +  C0*02-*N2-*H2 
3C02=C02/M3LES 
-*H20-H2  0/ MOLES 
?CO=C  0/ MOLES 

C  THE  INDEPENDENT  VARIABLE  IS  PARTIAL  PRESSURE  TIMES  MEAN  BEAM  LENGTH 
*C  02L-PCQ2 *L 
3H20L=PH20*L 
3C0L=PC0*L 
C 

C  CALCULATE  EMISSIVITY  OF  CARBON  DIOXIDE 
IF  <PC02L.LE.O. >  GO  TO  10 
3OUER=-0.721*0.215*ALOG10<PC02L> 

EMAX  =  -0.0A«-10.**POWER 
TMAX=26  J0.*800.*ALOG10<EMAX> 

3OUER  =  <  CT-TMAX 1/2  80 0. ) **2 
30WER=EXP<-P0UERI-1. 

ECO2=EMAX*10.**P0WER 
C  CORRECT  FOR  OVERLAP  WITH  WATER  VAPOR 
EC02=EC 02 -DELTA 
30  TO  15 
10  EC02=0. 

C  CALCULATE  EMISSIVITY  OF  WATER  VAPOR 
15  IF  (PH20L.LE.O. I  GO  TO  20 
EMAX=2.*ALOG10<PH20L> 

EMAX=-2.097*0.R21 *EMAX-Q  .  1 1 15 *E MAX* < EM AX -1 . I -0 . 005S7*EMAX 
1*<E MAX-1.) *< EMAX-2,  I 
3QWER=(EMAX-0.59)*(T*5000.)/700Q.*0.59 
EH20=10. “POWER 

C  CORRECT  FOR  PARTIAL  PRESSURE  OTHER  THAN  ZERO 
EH20=EH20*(l.-*<  0.62-0.1*  AL0G<PH20L>  )*PH20) 

30  TO  25 
20  EH20=0. 

C  CALCULATE  EMISSIVITY  OF  CARBON  MONOXIDE 
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25  IF  ( 3  COL*  LE  ■  0 • )  30  TO  30 

:>OWER=-0. 8477*0.1 809* ALOGIO(PCOL) 

E MAX=-0. 0 403  +  1 0. *  *POW£R 
TMAX= 1600. *3280.  *(0.122  -EMAX  ) 

30WER=(  (T-TMAX) /1800. >  **2 
30WER=-1.*EXP< -POWER) 

EC0=EMAX*10.**P0UER 

C  CORRECT  EMISSIVITY  OF  CO  FOR  OVERLAP  WITH  C02 
E C  0=0 •  7  *EC 0 
30  TO  35 
30  -:co=o. 

C 

C  CALCULATE  THE  TOTAL  EMI3SIVITY  OF  THE  COMBUSTION  GASES 
35  EMISS=EC02*EH20*EC0 

I F  (EMISS.LE.l. )  GO  TO  50 
WRITE (6*1) 

WRl TE (6»2) 

WRITE(S«3) 

50  RETURN 

1  =ORMAT( 1H0  »////SX ,"****..****•**•*••**•* ........... 


1................ ........ .......................................... 

2  r  OR  MA  T  f 1H  »//6X,"EMISSIVITY  GREATER  THAN  ONE . THE  GAS  COMPOSITI 

ION  IS  PR0BA9LY  INCORRECT") 

3  F0RMATC1H  « //6X »" **•* ************************ ***.**•*.....*.*.•** * 


1. ............................................................ .u// > 

ENO 
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SUBROUTINE  SHAPE< RI ,R 0,L »F  ,N> 

C  PROGRAM  TO  CALCULATE  THE  SHAPE  FACTORS  FOR  RAO  I  ANT  HEAT  TRANSFER 
C  BETWEEN  CONCENTRIC  CYLINDERS  OF  VARIOUS  LENGTHS  AND  ASPECTS 
C  THE  EXPRESSIONS  USED  HERE  WERE  DEVELOPED  FROM  EXPRESSIONS  FOR  SHA^E 
C  FACTORS  DERIVED  BY  LEUENBERGER  AND  3ERSON,  "C 3MP I L AT  ION  OF  RADIATION 
C  SHA 3E  FACTORS  FOR  CYLINDRICAL  ASSEMBLIES",  ASME  PAPER  5S-A-144  (1956) 
REAL  L»L1 

DIMENSION  F(50,50 > ,RI (50 > ,F1(2> 

DATA  PI/3. 1416/ 

C 

C  FIRST  CALCULATE  SHAPE  FACTORS  OF  DIRECTLY  OPPOSED  SEGMENTS 
00  5  1  =  1, N 

1=  (R I ( I) .LT.RO>  GO  TO  2 
-(1,1 >=1.00 
SO  TO  5 
2  ?=L 

IF  (( I.EQ.1>.0R.( I.EO.N) >  Z  =0  •  5  *  L 
A=Z**2-RO**2+RI(I>**2 
3=Z**2*R0**2-RI(I) *  *  2 
:  =  A/3 

D=Z**2*R0**2*RI  ( I ) *  *2 

r( I ,1 >=l.-l./PI*( ACOS(C> -0.5/RI(I)/Z*(SQRT(D**2-<2. *RI(I)*R0>**2> 

1  *AC03(C*RI ( I)/RO>  *A*ASIN(R I (I >/R0> -0.5*PI*3>  > 

5  CONTINUE 

C  AND  THEN  ALL  OTHERS 
00  40  1=1, N 
03  35  J=l «N 
IF  (I.E9.J)  GO  TO  35 
I-  (RI(I>.LT.RO>  GO  TO  10 
r ( I » J)  =  0. 

SO  TO  35 
10  CONTINUE 

IF  (J.GT.I)  GO  TO  20 
C 

C  FIRST,  SEGMENTS  RADIATING  "UPSTREAM" 

JJ=J*1 

11=1-1 

Ll  =  ( I-J  +  l >  *L 
IF  (I.EO.N)  LI =L1  -0 .5 *L 
00  15  K  =1 , 2 

A=L1**2-R0**2*RI(I>**2 

3=L1**2*R0**2-RI(I>**2 

C  =  A/B 

3  =  L1**2*.R0**2*RI(I>  **2 

rl(K)=RI(I)*Ll/PI/(R0**2-Rl(I)**2) * ( A C 0 3 ( C ) -0 . 5/R I (  I  > /L 1  * ( S ORT 
l(0**2-(2.*RI(I>*RO)**2>*AC0S(C*RI(I>/R0)*A*ASIN(RI(I>/R0> 

2-0 .5*PI *B>  > 

L1=L1-L 
15  CONTINUE 
Z=L 

Ir  (I.EO.N)  Z=0.5*L 

-(I,J)=0.5*(1.-F(  1,1) -(RO**2-RI(I)**2)/RI(I)/Z*(rl(l)-Fl(2)> ) 

C  NOTE  THAT  THESE  ARE  COMBINED  SHAPE  FACTORS  AND  NEED  TO  BE  SEPARATED 
SO  TO  35 
C 

20  CONTINUE 

C  AND  THEN  SEGMENTS  RADIATING  "DOWNSTREAM" 
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AD- A  1 39  457 


THERMAL  ENERGY  RECOVERY  IN  GAS  TURBINE  ENGINE  TEST 
CELLS(U)  NAVAL  CIVIL  ENGINEERING  LAB  PORT  HUENEME  CA 
C  A  KODRES  NOV  83  NCEl-TN-1679 


UNCLASSIFIED 


F/G  10/1 


NL 


n  o 


1 1  =  !♦ I 
JJ=J-1 

L1=<J-I*1 >*L 

IF  (I «EQ.l )  L1=L1-0.5*L 

33  25  K=  1  »  2 

3:L1**2^0**2-SI(IJ»*2 

c=A/a 

]:L1**2*R0»*2*RI(I)*»2 

rl<K)=RI<I>*Ll/Pr/fP0**2-RI< I)**2I*CACOSCC)-0.5/RI<I)/L1‘CSQRT 
1CD**2-C2.*RI<I>*RO»**2)*ACOS<C*RIC I ) /RO > ♦ A *AS I NC R I ( I)/RO) 
2-0.5*PI*B) > 

L1=CJ-I >*L 
25  CONTINUE 
Z=L 

I-  CI.EQ.l)  Z=3.5*L 

rU,J>  =  0.5*<l.-F<I,I)-<RO**2-RItI>**2>/RI<I)/Z*(Fl<l>-Fl<2>)> 
IF  CJ.LE.II)  GO  TO  35 
30  30  K=  1 1  »  J J 
r<ItJ)=F<I,J)-F(It<) 

30  CONTINUE 
35  CONTINUE 
40  CONTINUE 

SEPARATE  THE  "UPSTREAM"  SHAPE  FACTORS 
DO  55  1=3, N 
J J= I -2 
LL=I-1 

DO  50  M=1,JJ 
J= J J-M+1 
<  K  =  J* 1 

DO  45  K=KK,LL 
r<I,J)=FCI,JI-FCI,K) 

45  CONTINUE 
50  CONTINUE 
55  CONTINUE 
RETURN 
END 
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3JBR0UTINE  SHAPE1  <R  ,Z  ,F,  N  ) 

C  PR03RAM  TO  CALCULATE  THE  SHAPE  r ACTOR  rOR  RADIANT  HEAT  TRANSFER 
C  BETJEEN  FINITE  RING  AREAS  ON  INTERIOR  OF  RIGHT  CIRCULAR  CYLINDER 
C  REF...HOTTEL  AND  SAROFIM,  "RADIATIVE  TRANSFER",  MCGRAW-HILL,  (1967) 
DIMENSION  F<«S0,50) 

DATA  PI/3. 1416/ 

C 

DO  10  1=1, N 

DO  10  J  =1 »  N 

IF  (I.EQ.J)  GO  TO  5 

A  =  I ABSC I-JI  *Z 

B=< IABS<I-J>-1. >*Z 

:  =  ( IABSCI-J)*1.  >*Z 

-  ( I ,J>=PI/2.*<8**2-B*SQRT(4.*R**2*B**2>-2.*CA**2-A*SQRT(4.*R**2 
1*A**2>I*C**2-C*S0RT(4.*R**2*C**2)> 

30  TO  8 

5  -ri,J)=PI*C2.*R*Z+Z**2“Z*SQRTC4.*R**2+Z**?>) 

8  AREA=2.*PI *R*Z 
r  (  I ,J)=F< I, JD/ARFA 
10  CONTINUE 
RETURN 
END 
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TYPICAL  RESULTS 


This  particular  summary  is  for  the  J79  at  A/B  being  tested  in  the 
hush  house  at  the  Marine  Corps  Air  Station,  El  Toro,  Calif.  Waterwalls 
and  a  water  tube  boiler  are  examined  together  in  order  to  illustrate  the 
capabilities  of  the  model. 


augmented  geometry 

INITIAL  AREA  OF  JET  =  4.20  SOFT 

CROSS  SECTIONAL  AREA  OF  AUGMENTER  =  65.0  SOFT 

LENGTH  OF  AUGMENTER  =  67.0  FT 

THICKNESS  OF  AUGMENTER  WALLS  =  .111  IN 

THICKNESS  OF  ACOUSTIC  PILLOWS  =  6.00  IN 


PHYSICAL  CHARACTERISTICS  OF  JET  EXHAUST 

hass  flow  rate  of  jet  =  lao.  lb/sec 
TOTAL  TEMPERATURE  =  3060.  OEGF 
TOTAL  PRESSURE  :  36.8  PSIA 
RATIO  OF  SPECIFIC  HEATS  =  1.40 
COMPOSITION  OF  EXHAUST  IPERCENT  WEIGHTI 

OXYGEN .  0.00 

NITROGEN . 75.00 

CARBON  DIOXIDE . 16.00 

CARBON  MONOXIDE .  1.00 

WATER  VAPOR .  8.00 

HYDROGEN .  0.00 

LENGTH  OF  LUMINOUS  JET  =  0.0  FT 


PHYSICAL  CHARACTERISTICS  OF  SECONDARY  FLOW 

SEC3NDARY  FLOW  RATE  =  1130.  LB/SEC 
STATIC  TEMPERATURE  =  75.  OEGF 


AMBIENT  CONDITIONS 

TEMPERATURE  =  75.  DEGF 

CONVECTIVE  FILM  COEFFICIENT  =  5.0  BTU /HR-SOFT -DEGF 


MEAT  TRANSFER  PROPERTIES 

thermal  conductivity  of  augmenter  walls  =  30.0  btu/hr-ft-degf 

THERMAL  CONDUCTIVITY  OF  ACOUSTIC  PILLOWS  =  .050  BTU/HR-FT-DEGF 

EFFECTIVE  HUSH  HOUSE  CHARACTERISTICS 

TEMPERATURE  AT  AUGMENTER  INLET  =  70.  DEGF 

CONDUCTANCE  AT  AUGMENTER  INLET  =  .500E*03  BTU/HR-DEGF 

TEMPERATURE  AT  AUGMENTER  EXIT  =  150.  DEGF 
CONDUCTANCE  AT  AUGMENTER  EXIT  =  .500E»03  3TU/HR-DEGF 


BULK  VELOCITIES  AT  A JGHENTER  INLET 
JET  VELOCITY  =  2928.  FT/SEC 
SECONDARY  FLOW  VELOCITY  =  259.  FT/SEC 

CRAYA-CURTET  NUMBER  =  .62 


JET  FXHAJST  REACHES  THE  AUGMENTER  WALL  AT  AN  AXIAL  DISTANCE  OF  24.  FEET  FROM  THE  INLET 
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H»SS  FLOJ  RATES  TLB/SCC)  VARIATION  WITH  AXIAL  LOCATION  TFT)  AND  JET  URDUS  IFT> 


i 

i 

i 


t 


DISf  FM  INLET 

Q.00 

3. 72 

7.44 

11.17 

14.99 

18.61 

22.33 

26.06 

29.78 

33.50 

RAOIUS  0:  JET 

1*16 

1.64 

2.14 

2.64 

3.14 

3.66 

4.19 

4.  48 

4.48 

4.48 

JET  FLOW 

140.00 

262. SI 

375.07 

518.06 

693.92 

908.57 

1159.60 

1310.00 

1310.00 

1310.00 

SECONDARY  -LOW 

1130.00 

1047.39 

934.93 

791 .14 

614.08 

401  .43 

150.40 

•  00 

•  00 

•  00 

MASS  FLOW 

RATES 

I  LB/ S  E  C  ) 

VARIATION 

WITH  axial 

LOCATION 

f  FT  >  AND  JET 

RADIUS 

IFT1 

0 1 ST  FM  INLET 

37.22 

40.94 

44.67 

48.39 

52.11 

55.83 

59.56 

63.28 

67.00 

RAOIUS  Or 

JET 

4.48 

4.48 

4.48 

4.48 

4.48 

4.48 

4.48 

4.48 

4.40 

JET  FLOW 

1310.00 

1310.00 

1310  .00 

1310.00 

1310.00 

1310.00 

1310.00 

1310.00 

1310.00 

SECONOART 

FLOW 

.00 

•  00 

•  00 

.00 

.00 

•  00 

.00 

.00 

.00 

jlution  CONVERGED 

AFTER  32  ITERATIONS  TO  WirMfR  .84  0CG9CES 

Tf HPER ATjRf  C  DCGF 1 

)  V4RIA7ION  WITH 

AXIAL  LOCATION 

(F  T  1 . 

the  WALL  00 

INC LUDES 

ACOUSTIC 

PILLOWS,  IF 

ANT. 

OIST  FM  INLET 

0.00 

3.72 

7.44 

11.17 

14.99 

IB. 61 

22.33 

26.06 

29.78 

33.50 

JET  TEMP 

2?4  9.9 

1830.4 

141 7.1 

1096.2 

861.3 

690.0 

563.4 

509.6 

509.  T 

509.7 

sec  flow  temp 

75.0 

75.3 

75.6 

75.0 

75.9 

76.1 

76.5 

AUG  WALL  10 

117.4 

235.2 

264, r 

229.7 

1  90.6 

169.6 

173.9 

493.6 

498.9 

497.3 

AUG  WALL  00 

7  3.0 

75.0 

75.0 

75.0 

75.0 

7  5.0 

75.0 

75.0 

75.0 

75.0 

TEMPERATJRE  (OCGF  f 

variation 

WITH  AXIAL 

LOCATION  «FT). 

The  WALL  00 

INCLU  DES 

ACOUSTIC 

PILLOWS,  IF 

ANT. 

OIST  FM  INLET 

37.22 

40.94 

44.67 

48.39 

52.11 

55.83 

59.56 

63.28 

67.00 

JET  TEMP 

509.4 

509.1 

500.8 

508.5 

508.3 

508.1 

509.0 

50  8.2 

508.0 

AUG  WALL  10 

353.2 

353.2 

353.2 

353.2 

353.2 

353.2 

353.2 

353.2 

352.6 

AUG  WALL  00 

352.0 

352.0 

352.0 

J52.0 

352.0 

352.0 

352.0 

352.0 

352.0 

bulk  velocity  of  je 

T  AT  AUGMENTER  EXIT 

=  508. 

FT  fS€C 

EXHAUST  ilS 

EMISSIVITT  AND 

CONVECTIVE  FILM  COEFFICIENT 

IBTU/HR-SOFT-DEGF)  BT 

LOCATIOM 

OIST  FM  INLET 

0.00 

3.72  7.44 

11.17 

14.89  19.61 

22.33 

26.06 

29.78 

33.50 

FILM  CQErF 

5  5.15 

17.97  10.71 

7.70 

6.19  5.45 

5.58 

23.61 

23.61 

23.61 

E  MISS  IV  If  T 

.1711 

.2048  .1997 

•  1940 

.1666  .1510 

.1376 

•  1241 

•  1164 

•  1164 

EXHAUST  GAS  EMISSIVITT  AND 

CONVECTIVE  FILM  COEFFICIENT 

f  B  TU/HR-SQF  T -OEGFl  BT 

LOCATION 

OIST  FM  INLET 

37.22 

40.94  44.67 

48.39 

52.11  55.83 

59.56 

63.28 

67.00 

FILM  COErF 

23.61 

23.61  23.61 

23.61 

23*60  23*60 

23.60 

23.61 

23.61 

EMISSIVITT 

•  1164 

•1164  .1164 

•  1164 

.1164  .1164 

•  1164 

•  1164 

•  1164 
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THE  FACILITY  HAS  BOTH  A  VATER-UALL  BOILER  URAPPEO  AROUNO  THE  AUGMENTER  AMO 
A  STAGGERED  WATER-TUBE  BOILER  IN  THE  GAS  PATH  AT  THE  AUGMENTER  EXIT 


UATER-UALL  BOILER  CHARACTERISTICS 

STEAM  TEM^RATURE  =  552*  OCGF 

STEAM  PRESSURE  =  138.2  PSIA 

HEAT  OF  VALORIZATION  =  069.1  9TU/LB 

TOTAL  HEAT  TRANSFER  AREA  =  950.  SOFT 

THE  UATEt-UALLS  ARE  LOCATED  ALONG  THE  SECTION  BEGINNING  40.0  CEET  FROM  THE  AUGMENTER  INLET 
AND  ENDING  AT  A  OISTANCE  OF  67.0  FEET  FROM  THE  AUGMENTER  INLET 

TOTAL  STEAM  GENERATION  IN  UATER-UALL  BOILER  =  .400SE»04  LB/HR 

TOTAL  HEAT  TRANSFERRED  =  .3481E+07  BTU/HR 

UATER-UALL  STEAM  GENERATION  ( iB/HR )  VARIATION  UITH  AXIAL  LOCATION  (FT  F RDM  AUGMENTER  INLET) 


0 1ST  FM  INLET  0.00 
STEAM  GEN  0.0 
DIST  FM  INLET  37.22 
STEAM  GEN  2S9.7 


3.72  y .44  11.17  14.89  18.61 

0.0  0.0  0.0  0.0  0.0 

40.94  44.67  4R.39  52.11  55.83 

523.3  516.8  514.0  512-6  512-2 


FRACTION  OF  UATER-UALL  ^T CAM  GENERATION  BY  CONVECTION  AND  BY  RADIATION 


DIST  FM  INLET  0.00 
CONVECTION  0.000 
RAOIATION  0.000 


3.72  7.44  11.17  14.89  18.61 
0.000  0.000  0.000  0.000  0.000 
0.000  0.000  0.000  0.000  0.000 


FRACTION  OF  UA TER-UA. L  STEAM  GENERATION  BY  CONVECTION  AND  BY  RADIATION 


DIST  FM  INLET 

CONVECTION 

RADIATION 


40.94  44.67  48.39  52.11  55.83 

.961  .97?  .975  .977  .978 

.039  .028  .025  .023  .022 


CONVECTION  WATER-TUBE  BOILER  CHARACTERISTICS 
STEAM  TEMPERATURE  =  352.  DEGF 
STEAM  PRESSURE  =  138.2  PSIA 
HEAT  OF  VA»ORIZATION  =  869.1  BTU/LB 

TOTAL  HEAT  TRANSFER  AREA  s  1520.  SOFT 
OUTER  DIAMETER  OF  TUBES  =  3.00  IN 


TOTAL  STEAM  GENERATION  IN  CONVECTION  UATER-TURE  BOILER  =  .1320005  LB/H* 

AVERAGE  GAS  TEMPERATURE  AT  BOILER  INLET  =  509.  DEGF 

BOILER  OVERALL  HEAT  TRANSFER  COEFFICIENT  =  48.1  BTU/HR -SOFT- OEGF 

TOTAL  HEAT  TRANSFER*  i  0  =  «1147E»08  BTU/HR 


THE  ECONOMICS  OF  GAS  TURBINE  TEST  CELL  ENERGY  RECOVERY 


THE  PURCHASE  PRICE  Q-  THE  WATER-WALL  HEAT  EXCHANGER  IS  S1000000.00#  PAYABLE  WHEN  INSTALLATION  IS  CONPLCT 


The  ACCUMULATED  PRESENT  VALUE  OF  SAVINGS  RESULTING  FROM  THE  ADDITION  OF  THE  WATER-WALL  BOILER 
ASSUMING  A  DISCOUNT  RATE  OF  10.0  PERCENT*  A  CURRENT  STEAM  PRICE  OF  S  9. 71/MBTU •  AND  A 
STEAM  ESCALATION  RATE  OF  8.0  PERCENT 


ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  CYRS) 

( MRS/ 3A Y  > 

STEAM  SEN  <S> 

5 

•  25 

19652* 

.01 

10 

.25 

27981 • 

.03 

15 

•  25 

90159. 

•  09 

20 

•  25 

51271  . 

•  05 

25 

•  25 

61 9 OR  • 

•  06 

THE  PURCHASE  PRICE  0r  THE  CONVECTION  HEAT  EXCHANGER  IS  S2000000.00#  PAYABLE  WHEN  INSTALLATION  IS  COHPLET 


THE  ACCUMULATED  PRESENT  VALUE  OF  SAVINGS  RESULTING  FROM  THE  ADDITION  OF  THE  CONVECTION  BOILER 
ASSUMING  A  DISCOUNT  RATE  OF  10.0  PERCENT#  A  CURRENT  STEAM  PRICE  OF  t  9.75/MBTU#  AND  A 
STEAM  ESCALATION  RATE  OF  8.0  PERCENT 


ECONOMIC 

OPE  RAT  I  ON 

ACCUM  PV  OF 

SIR 

LIFE  I YR5  > 

CHRS/3AY) 

STEAM  GEN  IS) 

5 

•  25 

98203. 

•  02 

10 

.25 

92180. 

.05 

15 

.25 

132301. 

•  07 

20 

.25 

168906. 

•  08 

25 

.25 

202302. 

•  10 

UATER-WALL 

MEAT  EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  CYRSI 

C  HR  S/3A  Y  > 

STEAM  GEN  <*1 

5 

.50 

29263. 

•  03 

10 

.50 

55961. 

•  06 

15 

.50 

80319. 

•  08 

20 

.50 

102591. 

•  to 

25 

.50 

122816. 

.12 

CONVECTI3N 

C WATER  TU1E)  HEAT 

EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  CYRS 1 

CMRS/OAT) 

STEAM  GEN  CS> 

5 

•  50 

96905. 

•  05 

10 

.50 

189359. 

.09 

15 

•  SO 

269603. 

•  13 

20 

•  50 

337812. 

.17 

75 

•  59 

909609 . 

.’0 

D-28 
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WATER-WALL  HEAT  EXCHANGER 


economic 

OPERATION 

ACCOM  PV  OF 

SIR 

LIFE  I YRS ) 

(HRS/OAY) 

STEAM  GEN  (S) 

5 

1.00 

5952?. 

.06 

10 

1.00 

111923. 

.11 

15 

1.00 

160639. 

.16 

20 

1.00 

205083. 

.21 

25 

1.00 

245631. 

.25 

CONVECTION 

(WATER  TU9E  )  HEAT 

EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  (YRS) 

(HRS/OAY) 

STEAM  GEN  (S> 

5 

1.00 

192810 . 

.10 

10 

1.00 

368718. 

.19 

15 

1.00 

529206. 

.26 

20 

1.00 

675625. 

.34 

25 

1.00 

B0920R. 

.40 

WATER-WALL 

HEAT  EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  (YRS) 

(HRS/OAY) 

STEAM  GEN  (S) 

5 

2.00 

117053. 

.12 

10 

2.00 

223845. 

.2? 

15 

2.00 

321276. 

.32 

20 

2.00 

410165. 

.41 

25 

2.00 

491262. 

.49 

CONVECTION 

(WATER  TU3E  )  HEAT 

EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  (YRS) 

(HRS/OAY) 

STEAM  GEN  (S) 

5 

2.00 

385621. 

.19 

to 

2.00 

737437. 

.37 

15 

2.00 

1058412. 

.53 

20 

2.00 

1351249. 

.68 

25 

2.00 

1618416. 

.81 

WATER-WALL 

HEAT  EXCHANGER 

ECONOMIC 

OPERATION 

ACCUM  PV  OF 

SIR 

LIFE  (YRS) 

(HRS/OAY) 

STEAM  GEN  ($> 

5 

3.00 

175580. 

.19 

10 

3.00 

335768. 

.34 

15 

3.00 

491914. 

.48 

20 

3.00 

615248. 

.62 

25 

3.00 

736893. 

.74 
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convection  hater  tjtei  heat  exchanger 


ECONOMIC  OPERATION 


LIFE  trRS) 

1  HR  S/0  A  T  ) 

3 

3.00 

10 

3.00 

15 

3.00 

20 

3.00 

25 

5.01 

ACCUM  PV  OF 
STEAM  GEN  <*» 
578431  . 
1106155. 
1587  6 1 A • 
2026874. 
2427625. 


NATER-WAlL  heat  exchanger 


ECONOMIC 

OPERATION 

ACCUM  PV  OF 

LIFE  TYRS) 

THRS/OATl 

STEAM  GEN  C»> 

5 

5.00 

292633 . 

10 

5.00 

559613. 

15 

5.00 

803159. 

20 

5.00 

1025413. 

25 

5.01 

1228155. 

CONVECTION  HATER  TU3EI  HEAT  EXCHANGER 


ECONOMIC 

OPERA! I  ON 

ACCUM  PV  OF 

LIFE  TYRSI 

IHRS/1AY) 

STEAM  GEN  (»> 

5 

5.00 

964052. 

10 

5.00 

1843592. 

1  5 

5.01 

2646030. 

20 

5.00 

33  78 1 2  3 . 

25 

5.00 

4046039. 

UATER-WA.L 

HAT  EXCHANGER 

ECONOMIC 

OPERATION 

LIFE  IYR3I 

tmrs/jayi 

3 

8.00 

10 

8.00 

15 

8.00 

211 

a. no 

2  i 

a.  on 

ACCUM  PV  OF 
STEAM  GEN  (l) 
468215. 
895381 . 
12R5105. 
1640660. 
1965  04  9.' 


CONV  (  I  I  IN 


(  U  A  1 1  14  IUI  I  HT  A  T  EXCIfANOfR 


CC0M9MI; 

OPE  RATION 

ACCUM  PV  OF 

LIFE  CYRS! 

THRS/OATl 

STEAM  GEN  C* 

1 

a. on 

1542485. 

10 

a. oo 

2949748. 

1*) 

8.00 

4233649. 

20 

8.00 

5404957. 

PS 

8.00 

6473663. 

DISTRIBUTION  LIST 


ARMY  Fal  Engr,  Letterkenny  Army  Depot,  Chambersburg,  PA 
AF  SM/ALC7DEEEN  (J.  Pestillo)  McClellan  AFB,  CA 

AF  ENERGY  LIAISON  OFF-SERI  OFESC/OL-N  (Capt  B  Tolbert)  Golden  CO 

AF  HO  AF/LEY  Washington,  DC;  LEEH  (J  Stanton)  Washington.  DC;  PREES  Washington  DC  ;  Traffic 
Mgmnt  Cargo  Br  (LTC  S.  Rohrbough)  Washington, DC 
AFB  (AFIT/LDE),  Wright  Patterson  OH;  (RDVA)  AFESC/R&D  Tyndall.  FL;  82ABG/DEMC.  Williams  AZ; 
ABG/DEE  (F.  Nethers),  Goodfellow  AFB  TX;  AERO  Propulsion  Lab  (Energy  Conv)  Wright  Patterson  OH; 
AFESOTST.  Tyndall  FL;  AUL/LSE  63-465,  Maxwell  AL;  CESCH,  Wright-Patterson;  Elec  Sys  Div.  Code 
OCRS  (LTC  M.  Traister)  Hanscom,  MA;  HO  MAC/DEEE.  Scott.  II;  HO  Tactical  Air  Cmd/DEMM 
(Schmidt)  Langley.  VA;  MACOS-XOXC  (Col  Lee)  Scott.  IL;  SAMSO/MNND.  Norton  AFB  CA;  Samso/Dec 
(Sauer)  Vandenburg,  CA;  Stinfo  Library.  Offutt  NE;  WPNS  Safety  Div,  Norton.  CA;  Wright-Patterson. 
Encrgv  Conversion.  Dayton.  OH 
AFESC  DEB.  Tvndall.  FL;  HO.  RDVA  &  RDVCW 
AFWL  NTS  A.  Kirtland  AFB  NM 

AMMO  HAZARDS  REVIEW  BD  Chrmn.  NAVSAFEC'EN  (Code  43)  Norfolk.  VA 
ARCTICSUBLAB  Code  54,  San  Diego.  CA 

ARMY  (SA1L-FM)  OSA  Washington,  DC;  36th  Engineer  Group  AFVK-C  (Capt  P  Topp),  Ft  Benning.  GA; 
ARRADCOM.  Dover,  NJ:  ARRCOM.  Rock  Island.  IL:  BMDSC-RE  (H.  McClellan)  Huntsville  AL: 
Contracts  -  Facs  Engr  Directorate.  Fort  Ord.  CA;  DAEN-CWE-M.  Washington  DC:  DAEN-MCE-D. 
Washington.  DC:  DAEN-MPE-D  Washington  DC;  DAEN-MPO-U.  Washington.  DC':  DAEN-MPR.  Chief 
of  Engrs  Sol  Therm/Sol  Htg  &  Cool  Washington;  DAEN-MPU.  Washington  DC;  Def  Ammo  Cen  &  Scol 
(SARAC'-DEV)  R.  Hill.  Savanna.  IL:  ERADCOM  Tech  Supp  Dir.  (DELSD-L)  Ft  Monmouth.  NJ;  Engr 
District  (Memphis)  Library.  Memphis  TN;  FESA-TS  Fort  Belvoir  VA:  HO  DA.  Wash..  DC;  HO 
DAMA-CSC-T  Washington.  DC:  Hq  USAREUR.  Kaefertal.  Ger.;  Natick  R&D  Command  (Kwoh  Hu) 

Natick  MA;  Nuclear  &  Chemical  Agency,  Springfield.  VA;  OTRMSTR  Scol  (Code  ATSM-CD)  Ft.  Lee, 

VA;  R&D  Command  DRDNA-UST)  (J.  Siegel)  Natick.  MA:  STEYP-SE.  Yuma  Proving  Ground.  Yuma. 

AZ;  Savanna  Army  Depot  (Ernst).  IL:  Tech.  Ref.  Div.,  Fort  Huachuca.  AZ 
ARMY  -  CERL  Energy  Systems.  Champaign.  IL;  Library,  Champaign  IL 

ARMY  AMMUNITION  PLANT  SARHW-FEM  Hawthorne.  NV;  SARHW-FET  Hawthorne.  NV;  Sarhw  -  FEM 
Hawthorne.  NY 

ARMY  COE  Philadelphia  Dist.  (LIBRARY)  Philadelphia.  PA 

ARMY  CORPS  OF  ENGINEERS  Fae  Engr  Supp  Agency.  Ft.  Belvoir.  VA.  MRD-Eng  Div..  Omaha  NE; 

ORLCD-I  (R  Dockery)  Louisville.  KY;  Seattle  Dist.  Library.  Seattle  WA 
ARMY  CRREL  A.  Kovacs,  Hanover  NH:  G.  Phetteplace  Hanover.  NH:  Library.  Hanover  NH 
ARMY  CRREL  R  A  Eaton 

ARMY  DARCOM  Code  DRCMM-CS  Alexandria  VA 

ARMY  ENG  DIV  ED-CS  (S.Bolin)  Huntsville.  AL;  HNDED-CS.  Huntsville  AL;  HNDED-FD.  Huntsville.  AL: 
Huntsville.  AL 

ARMY  ENG  WATERWAYS  EXP  STA  (S  Kiger)  Vicksburg.  MS:  Library.  Vicksburg  MS 
ARMY  ENGR  DIST  Library.  Portland  OR 

ARMY  ENVIRON.  HYGIENE  AGCY  Bio  Acous  Div  Aberdeen  Proving  Ground  MD.  Dir  Env  Qual 

Aberdeen  Proving  Ground  MD;  HSE-EW  Water  Qual  Eng  Div  Aberdeen  Prov  Grnd  MD;  HSE-RP-HG/Pest 
Coord.  Arberdeen  Proving  Ground,  MD;  Librarian.  Aberdeen  Proving  Ground  MD 
ARMY  LOGISTICS  COMMAND  Code  ALC/ATCL-MS  (Morrissctt)  Fort  Lee.  VA 
ARMY  MAT  SYS  ANALYSIS  ACT  Code  DRXSY-CM  (M  Ogorzalck)  Aberdeen  Proving  Grnd  MD 
ARMY  MATERIALS  &  MECHANICS  RESEARCH  CENTER  Dr.  Lcnoe.  Watertown  MA 
ARMY  MISSILE  R&D  CMD  SCI  Info  Cen  (DOC)  Redstone  Arsenal.  AL 
ARMY  MOBIL  EOUIP  R&D  COM  DRDME-GS  Fuel  Tech  Br.  Ft  Belvoir.  VA 
ARMY  MPTRGN  CEN  Phys  Sec  D  A. L  E  T.  FT  McClellan  Al 
ARMY  MTMC  Trans  Engr  Agency  MTT-CE.  Newport  News.  VA 
ARMY  TRANSPORTATION  SCHOOL  Code  ATSPO  CD-TE  Fort  Eustis.  VA 
ARMY  TRNG  &  DOCTRINE  CMD  Code  ATCD-SL  Fort  Monroe.  VA 
ARMY-DEPOT  SYS  COMMAND  DRSDS-AI  Chambersburg.  PA 

ARMY-MERADCOM  DRDME-GE  (Heavy)  Ft  Belvoir  VA:  DRDME-MR  (J.  Sargent)  Ft  Belvoir  VA; 

DRDME-ZPS.  FT  BELVOIR  VA 
ASO  PWD  (ENS  M  W  Davis).  Phildadelphia.  PA 

ASST  SECRETARY  OF  THE  NAVY  Spec.  Assist  Submarines.  Washington  DC 
BUMED  Security  Offr.  Washington  DC 

BUREAU  OF  RECLAMATION  (J  Graham).  Denver.  CO;  Code  1512  (C  Selandcr)  Denver  CO 

CHNAVRES  Code  002C  New  Orleans.  LA 

CINCLANT  CIV  ENGR  SUPP  PLANS  OFFR  NORFOLK.  VA 

CINCLANTFLT  Code  J54I  Asst.  Nuc.  Warfare  Plans  Off  Norfolk.  VA;  NUC  Wpns  Sec  Offr  Norfolk.  VA 
CINCPAC  Fac  Engmg  Div  (J44)  Makalapa.  HI 


CINCUSNAVEUR  NUC  Wpns  Sec  Offr  London,  England 
CNAVRES  Code  13  (Dir.  Facilities)  New  Orleans.  LA 

CNM  Code  03462  Washington  DC:  Code  03462.  Washington  DC;  Code  03ZB  Washington  DC;  Code  043 
Washington  DC;  Code  043  Washington  DC;  Code  MAT-04.  Washington,  DC;  Code  M AT-OKE, 

Washington.  DC;  MAT  08T245  (Spalding).  Washington.  DC:  MAT-0718.  Washington,  DC;  NMAT  -  044. 
Washington  DC 

CNO  Code  N OP-464,  Washington  DC;  Code  OP  323.  Washington  DC:  Code  OP  405.  Washington  DC;  Code  OP 
405.  Washington,  DC:  Code  OP  414.  Washington  DC;  Code  OP  97  Washington  DC;  Code  OP  97 
Washington.  DC;  Code  OP  987  Washington  DC;  Code  OP-413  Wash.  DC;  Code  OP323  Washington  DC; 

Code  OP403.  Washington  DC;  Code  OPNAV  09B24  (H);  Code  OPNAV  22.  Wash  DC;  Code  OPNAV  23. 
Wash  DC:  OP-098.  Washington.  DC;  OP-23  (Capt  J.H  Howland)  Washinton.  DC;  OP-403DI,  Washington. 
DC;  OP-41  IF.  Wash  DC:  OP987J.  Washington.  DC 
C'OMASWWINGPAC  Security  Offr.  San  Diego.  CA;  Weapons  Offr.  San  Diego,  CA 
COMCBPAC  Code  CB50.  Makalapa.  Oahu.  HI:  Makalapa.  Oahu.  HI:  Operations  Off.  Makalapa  HI 
COMDEVGRUONE  CMDR  San  Diego.  CA 
COMFAIRWESTPAC  Security  Offr.  Misawa  Japan 

COMFLEACT.  OKINAWA  PWD  -  Engr  Div.  Saseho.  Japan:  PWO.  Kadcna.  Okinawa;  PWO.  Sasebo,  Japan: 
Security  Offr.  Japan 

COMNAVAIRLANT  Code  335  Norfolk.  VA;  NUC  Wpns  Sec  Offr  Norfolk.  VA 

COMNAVAIRPAC  FOR  Wpns  Offr  (Code  93)  NAS  North  Island:  NUC  Wpns  Sec  Offr  NAS.  North  Island 
COMNAVBEACHPHIBREFTRAC.ru  ONE  San  Diego  CA 
COMNAVMARIANAS  Code  N4.  Guam 
COMNAVSURFLANT  Norfolk.  VA 

COMNAVSURFPAC  NUC  Wpns  Sec  Offr  NAB.  Coronado 
COMOCEANSYSLANT  PW-FAC  MC.MNT  Off  Norfolk.  VA 
COMOCEANSYSPAC  SCE.  Pearl  Harbor  HI 
COMRNCF  Nicholson.  Tampa.  FL;  Nicholson.  Tampa.  FL 
COMSUBDEVGRUONE  Operations  Offr.  San  Diego.  CA 
COMSUBLANT  NUC  Wpns  Sec  Officer  Norfolk.  VA 
COMUSNAVCENT  N42.  Pearl  Harbor.  HI 
NAVSURFLANT  NUC  Wpns  Sec  Officer  Norfolk.  VA 
NAVSURFPAC  Code  N-4.  Coronado 

COMOPTEVFOR  C  MDR.  Norfolk.  VA:  Code  705.  San  Diego.  CA 
DEFENSE  CIVIL  PREPAREDNESS  AGENCY  Washington.  DC 
DEFENSE  INTELLIGENCE  AGENCY  DB-4CI  Washington  DC 

DEFFUELSLIPPCEN  DFSC-OWE  (Term  Engrng)  Alexandria.  VA:  DFSC-OWE.  Alexandria  VA 

DOE  E  Hightower.  Washington.  DC:  G.F.  Hedrick.  Washington.  DC 

DIRSSPO  OP  Engrng  Sect  (Lopata).  Washington.  DC;  Security  Officer.  Washington.  DC 

DLA  DLA-TP.  Alexandria.  VA:  Industrial  Security.  Alexandria.  VA 

DLSIE  Army  Logistics  Mgt  Center.  Fort  Lee.  VA 

DNA  Field  Command.  Kirtland  AFB  NM;  STTL.  Washington  DC 

DOD  Explos  Safety  Board.  Wash  DC;  Explosives  Safety  Board  (Library).  Washington  DC:  Staff  Spec.  Chcm. 
Tech.  Washington  DC 

DOE  (G.  Boyer).  Washington.  DC;  Div  Ocean  Energy  Svs  Cons/Solar  Energy  Wash  DC:  INEL  Tech.  Lib. 

(Reports  Section).  Idaho  Falls.  ID;  S&S  Div  (G.  Hedrick)  Washington  DC 
DOT  Office  of  Hazardous  Matl  Ops.  Wash  DC 
DTIC  Defense  Technical  Info  Ctr<' Alexandria.  VA 

DTNSRDC  Anna  Lab  (Code  119)  Annapolis  MD;  Anna  Lab  (Code  1568)  Annapolis  MD:  Anna  Lab.  Code 
2724  (D  Bloomquist)  Annapolis.  MD:  Anna  Lab.  Code  2842  (Mr.  Sherman)  Annapolis.  MD:  Anna  Lab. 

Code  4121  (R  A  Rivers)  Annapolis.  MD 
DTNSRDC  Code  1706.  Bethesda  MD;  Code  172  (M.  Krenzke).  Bethesda  MD 
DTNSRDC  Code  284  (A.  Rufolo),  Annapolis  MD 

DTNSRDC  Code  4111  (R.  Gierich).  Bethesda  MD;  Code  42.  Bethesda  MD 

DTNSRDC  Code  522  (Library).  Annapolis  MD;  Pattison.  Code  1706.  Bethesda.  MD;  Rispin.  Code  1706, 

Bethesda.  MD 

ELECTRONIC  WARFARE  LAB.  Intel  Mat  Dev  &  Spt  off.  Ft  Geo  G.  Meade.  MD  20755;  Intel  Mat  Dev/Supp 
Off  Ft  Geo  G  Meade  MD 

ENVIRONMENTAL  PROTECTION  AGENCY  ANR-458  (Dr  R.S.  Dyer)  Washington.  DC:  Off.  of  Air 
Programs  (E.  Stork).  Ann  Arbor  MI;  Reg.  Ill  Library,  Philadelphia  PA;  Reg.  VIII.  8M-ASL.  Denver  CO 
FAA  ACS-200.  Washington.  DC;  ACS-250  Washington  DC;  ACS-300  Washington.  DC 
FAA  Civil  Aeromedical  Institute.  Oklahoma  City.  OK 
FED  RAILROAD  AMIN  Administrator.  Wash  DC 
FLTCOMBATTRACENLANT  PWO.  Virginia  Bch  VA 
FMFLANT  CEC  Offr.  Norfolk  VA 
FMFPAC  CG(FEO)  Camp  Smith.  HI 

GSA  Assist  Comm  Dcs  &  Cnst  (FAIA)  D  R  Dibner  Washington.  DC  ;  FPS-DIR  Supp  Services.  Washington. 


/ 


DC;  FPS-Dir  Supp  Serv  (Q  Y  Lawson)  Washington  DC;  Off  of  Des  &  Const-PCDP  (D  Eakin)  Washington, 
DC;  Off  of  Fed  Prot  Serv  Mgmt  Washington,  DC;  PUB  BMg  Serv  (W  D  Rust  Jr)  Washington  DC;  Pub  Bldg 
Serv.  Washington  DC 
HCU  ONE  CO,  Bishops  Point,  HI 
HC  &  RS  Tech  Pres.  Service.  Meden.  Washington,  DC 
HO  UNC/USFK  Korea  Security  Officer 

JOINT  CRUISE  MISSILE  PROJ  OFF  JCM22  Washington,  DC;  PM  3  Technical  Library,  Washington,  DC 
KWAJALEIN  MISRAN  BMDSC-RKL-C 

LIBRARY  OF  CONGRESS  Washington,  DC  (Sciences  &  Tech  Div) 

MARCORPS  Code  LME-1  (Riggs)  Washington,  DC 

MARINE  BARRACKS  CO  Brunswick,  ME;  CO  Colts  Neck.  NJ;  CO  Yorktown.  VA,  CO.  Adak,  AK,  CO. 
Alameda,  CA;  CO.  Barbers  Point.  HI;  CO,  Bermuda;  CO.  Cecil  Fid.  FL;  CO.  Charleston.  SC;  CO. 

Concord.  CA;  CO.  Jacksonville.  FL;  CO.  Keyport.  WA;  CO.  Lualualei.  Waianae.  HI;  CO.  Moffett  Fid.  CA; 
CO.  Seal  Beach.  CA;  NAS.  Norfolk.  VA 

MARINE  CORPS  BASE  1st  For  Serv  Supp  Gru  (CSS-5)  Camp  Pendleton  CA;  Code  406.  Camp  Lejeune,  NC; 

M  &  R  Division.  Camp  Lejeune  NC;  Maint  Off  Camp  Pendleton,  CA;  PWD  -  Maint.  Control  Div.  Camp 
Butler.  Kawasaki.  Japan;  PWO  Camp  Lejeune  NC;  PWO,  Camp  Pendleton  CA;  PWO.  Camp  S.  D.  Butler. 
Kawasaki  Japan 

MARINE  CORPS  HQS  Code  LFF-2.  Washington  DC;  Code  LM-2.  Washington.  DC;  Code  MMEA 
Washington,  DC;  Code  MPH-50.  Washington.  DC;  Code  OTTI  Washington  DC 
MARITIME  ADMIN  (E.  Uttridge),  Washington.  DC 

MCAS  Facil.  Engr.  Div.  Cherry  Point  NC;  CO.  Kaneohe  Bay  HI;  Code  S4.  Quantico  VA;  Facs  5  )ept  - 
Operations  Div.  Cherry  Point;  PWD  -  Utilities  Div,  Iwakuni.  Japan:  PWD.  Dir.  Maint.  Contro 
Iwakuni  Japan:  PWO.  Iwakuni,  Japan;  PWO.  Yuma  AZ;  SCE.  Futema  Japan 
MCDEC  M&L  Div  (LTC  R.  Kerr)  Quantico.  VA;  M&L  Div  Quantico  VA;  NSAP  REP,  Ouantic  A 
MCLB  B520,  Barstow  CA:  Maintenance  Officer.  Barstow.  CA:  PWO.  Barstow  CA 
MCRD  SCE.  San  Diego  CA 

MILITARY  SEALIFT  COMMAND  Code-M-OIFB.  Washington.  DC;  Washington  DC 
NAF  PWD  -  Engr  Div,  Atsugi.  Japan;  PWO,  Atsugi  Japan 
NALCOEASTPAC  Security  Officer,  Alameda.  CA;  Security  Offr 
NALF  OINC.  San  Diego.  CA 

NARF  Code  100.  Cherry  Point.  NC:  Code  612.  Jax.  FL:  Code  640.  Pensacola  FL;  Equipment  Engineering 
Division  (Code  61000).  Pensacola.  FL;  SCE  Norfolk.  VA 
NAS  AUW  Officers.  Brunswick.  ME:  CO.  Guantanamo  Bay  Cuba:  Code  114.  Alameda  CA;  Code  183  (Fac. 

Plan  BR  MGR):  Code  187,  Jacksonville  FL:  Code  18700.  Brunswick  ME;  Code  18U  (ENS  P.J.  Hickey). 
Corpus  Christi  TX;  Code  502  (SSF  Officer).  Norfolk.  VA;  Code  6234  (G.  Trask).  Point  Mugu  CA;  Code  70. 
Atlanta.  Marietta  GA;  Code  8E,  Patuxent  Riv.,  MD;  Dir  of  Engrng.  PWD,  Corpus  Christi.  TX;  Dir  of 
Security.  Cecil  Fid.  FL:  Dir.  Maint.  Control  Div..  Key  West  FL:  Dir.  Util.  Div..  Bermuda:  Grover.  PWD. 
Patuxent  River.  MD;  Lakchurst,  NJ;  Lead.  Chief.  Petty  Offr.  PW/Self  Help  Div.  Beeville  TX:  OIC.  CBU 
417,  Oak  Harbor  WA;  PW  (J.  Maguire).  Corpus  Christi  TX:  PWD  -  Engr  Div  Dir.  Millington.  TN;  PWD  - 
Engr  Div.  Gtmo,  Cuba:  PWD  -  Engr  Div,  Oak  Harbor.  WA:  PWD  -  Maint.  Control  Dir.  Millington.  TN; 
PWD  Maint.  Cont.  Dir..  Fallon  NV;  PWD  Maint.  Div..  New  Orleans.  Belle  Chasse  LA:  PWD.  Code  I821H 
(Pfankuch)  Miramar.  SD  CA.  PWD.  Maintenance  Control  Dir..  Bermuda:  PWD.  Willow  Grove  PA:  PWO 
Belle  Chasse.  LA;  PWO  Chase  Field  Beeville.  TX;  PWO  Key  West  FL;  PWO  Lakchurst.  NJ:  PWO  Sigonella 
Sicily:  PWO  Whiting  Fid.  Milton  FL;  PWO,  Dallas  TX;  PWO.  Glenview  IL;  PWO.  Kingsville  TX;  PWO. 
Millington  TN:  PWO.  Miramar.  San  Diego  CA;  PWO.,  Moffett  Field  CA;  SCE  Norfolk.  VA;  SCE.  Barbers 
Point  HI;  SCE,  Cubi  Point.  R.P:  Security  Officer.  Norfolk.  VA;  Security  Offr.  Alameda  CA;  Security  Offr. 
Barbers  Point.  HI;  Security  Offr.  Bermuda;  Security  Offr,  Guam;  Security  Offr.  Key  West  FL:  Weapons 
Officer.  Norfolk.  VA:  Weapons  Offr.  Agana,  Guam.  MI:  Weapons  Offr.  Alameda.  CA;  Weapons  Offr. 
Barbers  Point.  Oahu.  HI:  Weapons  Offr.  Moffett  Fid.  CA;  Weapons  Offr.  North  Island:  Wpns  Offr.  Cecil 
Fid.  FL:  Wpns  Offr.  Key  West.  FL 
NASA  LBJ  Space  Cen  (Alexander)  Houston.  TX 
NATL  ACADEMY  OF  SCIENCES  Polar  Res  Com  (DeGoes) 

NATL  BUREAU  OF  STANDARDS  B-348  BR  (Dr.  Campbell).  Washington  DC;  Kovacs.  Washington.  D  C; 

Law  Enforce  Stds  Lab  Washington.  DC;  R  Chung  Washington.  DC 
NATL  RESEARCH  COUNCIL  Naval  Studies  Bd,  Washington.  DC;  Naval  Studies  Board.  Washington  DC 
NATL  SEC  AGENCY  Ft  Geo  G  Meade  MD.  H  V  Powell  Jr  FT  Geo  G.  Meade  MD 
NAVACT  PWO.  London  UK 
NAVACTDET  PWO.  Holy  Lock  UK 
NAVAEROSPREGMEDCEN  SCE.  Pensacola  FL 

NAVAIRDEVCEN  Chmielewski.  Warminster,  PA;  Code  813.  Warminster  PA:  PWD.  Engr  Div  Mgr, 

Warminster,  PA;  Security  Offr.  Warminster  PA 
NAVAIRPROPTESTCEN  CO.  Trenton.  NJ 

NAVAIRSYSCOM  Code  417I2A  (Michols)  Washington.  DC:  PMA  258  Washington,  DC 
NAVAIRTESTCEN  PATUXENT  RIVER  PWD  (F  McGrath).  Patuxent  Riv  .MD 
NAVAVIONICFAC  PW  Div  Indianapolis.  IN;  PWD  Deputy  Dir.  D/701.  Indianapolis.  IN 


NAVAVNWPNSFAC'  Sec  Offr  St.  Mawgan,  England:  Wpns  Offr.  St.  Mawgan.  England 

NAVCHAPGRU  CO  Williamsburg  VA;  Engineering  Officer.  Code  60  Williamsburg.  VA;  Operations  Officer. 

Code  30  Williamsburg,  VA;  Supply  Officer.  Code  70  Williamsburg.  VA 
NAVCOASTSYSCEN  CO.  Panama  City  FL;  Code  47.3  Panama  City.  FI.;  Code  713  (J  Ouirk)  Panama  City.  FE: 
Code  713  (J.  Mittleman)  Panama  City.  FL,  Code  710.  Panama  City.  FL.  Code  772  (C  B  Koesy)  Panama  City 
FL.  Library  Panama  City.  FL.  PWO  Panama  City.  FL 
NAVCOASTSYSTCTR  Security  Offr.  Panama  City  FL. 

NAVCOMMAREAMSTRSTA  Code  W-bO.  Elec  Engr.  Wahiawa.  Ill;  Maim  Control  Div.  Wahiawa.  HI;  PWO. 

Norfolk  VA;  SCE  Unit  I  Naples  Italy;  SCE.  Wahiawa  HI 
NAVCOMMAREAMSTRSTA  Security  Offr.  Norfolk  VA 

NAVCOMMSTA  CO.  San  Miguel.  R.P.:  Code  401  Nea  Makri.  Greece;  PWD  •  Maim  Control  Div.  Diego 
Garcia  Is.:  PWO.  Exmouth.  Australia;  SCE  Kcflavik.  Iceland;  SCE  Rota  Spain;  SCE.  Balboa.  CZ;  Security 
Offr.  Stockton  C'A 

NAVCONSTRACEN  Code  00U13.  Port  Hueneme  CA;  Curriculum  Instr.  Stds  Offr.  Gulfport  MS 
NAVEDTRAPRODEVCEN  Technical  Library.  Pensacola.  FL 

NAVEDUTRAC'EN  Engr  Dept  (Code  42)  Newport.  Rl:  PWO.  NAVBASE.  Newport.  R1 
NAVELEXSYSCOM  Code  PME  124-61.  Washington.  DC:  Commander.  Washington.  DC:  PME  124-612.  Wash 
DC;  PME-I2I  Washington.  DC 

NAVENENVSA  Code  II  Port  Hueneme,  CA;  Code  11 1 A  (Winters)  Port  Hueneme  CA 
NAVENVIRHLTHCEN  CO.  NAVSTA  Norfolk.  VA 
NAVEODTECHCEN  Code  603.  Indian  Head  MD 

NAVFAC  PWO.  Brawdy  Wales  UK;  PWO.  Centerville  Bch,  Ferndalc  CA;  PWO.  Point  Sur.  Big  Sur  CA 
NAVFACENGCOM  Alexandria.  VA;  Code  03  Alexandria.  VA;  Code  03T  (Essoglou)  Alexandria.  VA:  Code 
(U4B)  Alexandria.  VA.  Code  0431  (P  W  Brewer)  Alexandria.  Va;  Code  0431.  Alexandria.  VA:  Code 
0432D.  Alexandrai.  VA;  Code  0433  (D  Potter)  Alexandria.  VA:  Code  0433C.  Alexandria.  VA:  Code  0434B 
Alexandria.  Va:  Code  046;  Code  0461D  (V  M  Spaulding)  Alexandria.  VA:  Code  04A1  Alexandria.  VA: 

Code  04B  (M  Yachnis)  Alexandria.  VA:  Code  04B3  Alexandria.  VA:  Code  031.  Alexandria.  VA;  Code 
031A  Alexandria.  VA:  Code  06.  Alexandria  VA.  Code  063T.  Alexandria.  VA;  Code  09V134.  Tech  Lib. 
Alexandria.  VA:  Code  100  Alexandria.  VA:  Code  1002B  (J.  Leimanisl  Alexandria.  VA:  Code  1113. 
Alexandria.  VA;  Code  1 1  IB  Alexandria.  VA:  Code  1 1  IN  Alexandria.  VA.  Code  461D.  Alexandria.  VA; 
code  (1ST  Alexandria.  VA 

NAVFACENGCOM  -  CUES  DIV.  Code  03.  Wash.  DC:  Code  10  11  Washington  DC;  Code  101  Wash.  DC: 

Code  403  Washington  DC;  Code  403  Wash.  DC;  Code  407  (D  Schcesele)  Washington.  DC;  Code  FPO-IC 
Washington  DC;  Code  FPO-IE.  Wash  DC;  Contracts.  ROICC.  Annapolis  MD;  FPO-I  Washington.  DC: 
FPO-IEA3  Washington  DC;  FPO-IP  IP3  Washington.  DC:  Library.  Washington.  D  C 
NAVFACENGCOM  -  LANT  DIV  Code  04  Norfolk  VA  Norfolk  VA.  Code  05.  Norfolk.  VA;  Code  II. 

Norfolk.  VA:  Code  111.  Norfolk.  VA.  Code  403.  Norfolk.  VA:  Code  403  Ci.  '  Engr  BR  Norfolk  VA;  Code 
40S.  Norfolk.  VA;  Eur  BR  Deputy  Dir.  Naples  Italy;  Library,  Norfolk.  VA:  FDT&ELO  I02A.  Norfolk. 

VA 

NAVFACENGCOM  -  NORTH  DIV  (Boretsky)  Philadelphia.  PA:  CO;  Code  04  Philadelphia.  PA:  Code  04AL. 
Philadelphia  PA;  Code  03.  Phila.  PA;  Code  09P  Philadelphia  PA:  Code  111  Philadelphia.  PA:  Code  114  (A 
Rhoads)  Philadelphia.  PA:  Code  114  (A.  Rhoads):  Code  403  Philadelphia.  PA:  ROICC.  Contracts.  Crane 
IN 

NAVFACENGCOM  -  PAC  DIV.  (Kyi)  Code  101.  Pearl  Harbor.  HI:  CODE  OOP  PEARL  HARBOR  HI;  Code 
04  Pearl  Harbor  HI;  Code  03.  Pearl  Harbor.  HI:  Code  II  Pearl  Harbor  HI:  Code  2011  Pearl  Harbor.  HI; 
Code  402.  RDT&E.  Pearl  Harbor  HI;  Code  403  Civil  Engr  BR  Pearl  Harbor  HI:  Code  406  Pearl  Harbor.  HI 
Commander.  Pearl  Harbor.  HI:  Library,  Pearl  Harbor.  HI 
NAVFACENGCOM  -  SOUTH  DIV.  Code  04.  Charleston.  SC;  Code  03.  Charleston.  SC;  Code  00  (Watts) 
Charleston.  SC;  Code  11.  Charleston.  SC;  Code  403.  Gaddy.  Charleston.  SC;  Code  403  Charleston.  SC; 

Code  411  Soil  Mech  &  Paving  BR  Charleston.  SC;  Code  90.  RDT&ELO.  Charleston  SC:  Library. 

Charleston.  SC 

NAVFACENGCOM  -  WEST  DIV.  102;  AROICC.  Contracts.  Twentynine  Palms  CA;  Code  04.  San  Bruno.  CA; 
Code  04B  San  Bruno.  CA;  Code  03.  San  Bruno.  CA:  Code  101.6  San  Bruno.  CA;  Code  114C.  San  Diego 
CA;  Code  402  San  Bruno.  CA:  Library.  San  Bruno.  CA;  O9P/20  >.„t  Bruno.  CA;  RDT&ELO  San  Bruno. 
CA;  Security  Offr.  Seattle  WA 

NAVFACENGCOM  CONTRACTS  AROICC.  NAVSTA  Brooklyn.  NY:  AROICC.  Quantico.  VA.  Colts  Neck. 
NJ:  Contracts.  AROICC.  Lemoore  CA;  Dir.  Eng.  Div..  Exmoulh.  Australia:  Eng  Div  dir.  Southwest  Pac. 
Manila,  PL  Engr  Div.  (F  Hein).  Madrid.  Spain;  NAS.  Jacksonville.  FL:  OICC.  Southwest  Pac.  Manila.  PI: 
OICC-ROICC.  NAS  Oceana.  Virginia  Beach.  VA;  OICC/ROICC.  Balboa  Panama  Canal;  ROICC  AF 
Guam;  ROICC  Code  493  Portsmouth  VA.  ROICC  Key  West  FL;  ROICC  MCAS  El  Toro;  ROICC  Rota 
Spain:  ROICC'.  Diego  Garcia  Island;  ROICC.  Kcflavik.  Iceland;  ROICC.  NAS.  Corpus  Christi.  TX; 

ROICC.  Pacific.  San  Bruno  C'A;  ROICC.  Point  Mugu.  CA;  ROICC.  Silverdale.  WA;  ROICC.  Yap. 
ROICC'-OIC’C-SPA.  Norfolk.  VA 
NAVFORCARIB  Commander  (N42).  Puerto  Rico 
NAVHOSP  PWD  -  Engr  Div.  Beaufort.  SC 

NAVINVSERVHQ  NIS-243  (J  OHara)  Washington.  DC;  Security  Officer.  Alexandria.  VA 


NAVMAG  PWD  •  Engr  Div.  Guam;  SCE.  Guam;  SC  I'.  Subic  Bay.  R.P  ;  Security  Offr.  Guam 
NAVMAG  Security  Offr,  Lualualci.  Ill  Waianae.  Ill 
NAVMAG  Wpns  Offr.  Guam.  Mariana  island 

NAVOC'EANO  Code  .143:  (J.  DoPalma).  Bay  St.  Louis  MS;  library  Bay  St  Louis.  MS 
NAVOC'EANSYSC'EN  C  ode  (N  (Talkington).  San  Diego.  CA:  C  ode  4473  If. o side  Library.  San  Diego.  CA; 

Code  4473B  (Tech  Lib)  San  Diego.  C'A;  Code  53(14  (J  Stachiw).  San  Diego,  CA.  Code  5314  (II  Wheeler). 
San  Diego  CA;  Code  5331  (R Jones)  San  Diego  Ca;  Code  533  (Hurley).  San  Diego.  CA;  Code  5333 
(Baehman)  San  Diego.  CA;  Code  6700.  San  Diego.  CA;  Code  SI  I  San  Diego.  CA;  llauaii  Lab  (R  Yumoril 
Kailua.  HI;  Hi  Lab  Tech  Lib  Kailua  HI 
NAVORDFAC  CO  (Code  66).  Sasebo.  Japan 
NAVORDMISTESTFAC  PWD  -  Engr  Dir.  White  Sands.  NM 

NAVORDSTA  PWD  -  Dir.  Engr  Div.  Indian  Head.  Ml);  PWC).  I  omsville  KY;  Seeunty  Oflr.  Indian  Head  MD 
Security  Offr,  Indian  Head  MD 

NAVPETOFF  C'ode  30.  Alexandria  VA;  Code  XDI07.  Alexandria.  VA 
NAVPETRES  Director.  Washington  DC 

NAVPGSCOL  (Dr.  G.  Haderlie).  Monterey.  CA;  C  Morers  Monterey  CA;  Code  nlWI.  (O  Wilson)  Monterey 
CA;  F.  Thornton.  Monterey  CA 

NAVPHIBASE  CO.  AC'B  3  Norfolk.  VA;  COMNAVBEAC  IKiRC  I  WO  Norlolk  VA;  C  ode  S3 1 .  Norfolk  VA; 
Dir  Amphib.  Warfare  Brd  Staff.  Norfolk.  VA;  Harbor  Clearance  l  nit  Two.  Little  Creek.  VA;  SCT 
Coronado.  SD.CA.  Security  Offr.  Norfolk  VA 
NAVRADRECFAC  PWO.  Kami  Seva  Japan 

NAVRECiMEDCEN  Chief  of  Police.  Camp  Pendleton  CA;  Code  3d.  I  ns  Health  Sers .  (Al  Bryson)  San  Diego. 
CA 

NAYHOSP  CO.  Millington.  IN 

NAVRECiMEDCEN  PWD  -  Engr  Div.  Camp  l.ejeune.  NC;  PWO  Portsmouth.  VA;  PWO.  Camp  l.ejcune.  NC' 
NAVRECiMEDCEN  PWC).  Okinawa.  Japan 

NAVRECiMEDCEN  SCE:  SCE  San  Diego.  CA.  SCE.  Camp  Pendleton  CA;  SCT.  Guam.  SCI  .  Newport.  Rl; 
SCE.  Oakland  C'A 

NAVRECiMEDCEN  SC'E.  Yokosuka.  Japan 
NAVRECiMEDCEN  Security  Offr.  Oakland  CA 
NAVREGMEDC'I  INK'  A  Watanahe.  Pearl  Harbor.  HI 

NAVSCOLCECOFF  C35  Port  Ilucnemc.  C'A;  CO.  C  ode  C44A  Port  llueneme.  CA 
NAVSCSOL  PWO.  Athens  CiA 
NAVSEASUPPCENI.ANT  Code  S(K)  Portsmouth  VA 

NAVSEASYSCOM  (R  Seal.  Washington.  DC;  Code  04J33  Washington.  IK';  Code  05E1.  Wash.  DC;  Code 
051.313  (Ncwhouse)  Wash  DC;  Code  6313.  Wash  DC;  Code  033  (Mr  J  Peters)  Washington.  DC;  Code 
OOC-D.  Washington.  DC;  Code  PMS  305  A  3.  Washington.  DC.  Code  PMS  305  A3.  Washington.  DC; 

Code  PMS  3%. 3311  (Rekas),  Wash..  DC;  Code  SEA  OOC  Washington.  DC;  E  Wuigk.  Washington.  DC; 
MO  NO.  Washington.  DC;  PMS- 305  Al.  Washington.  DC;  PMS.W5-A3.  Washington.  DC:  SEA  (ME  (I. 

Kess)  Washington.  DC;  SEA  05L313.  Washington.  DC;  SEA-0411.  Wash  DC.  SEA-05R3  (J  Freund) 
Washington.  DC;  SEA-5433.  Washington.  DC:  SEA-0043  Washington.  DC:  SEA-00611.  Washington.  DC; 
SEA05E1.  Washington.  D  C. 

NAVSEC  Code  6156D.  Washington.  DC;  Code  6I57D.  Washington.  DC 

NAVSECGRUACT  Facil.  Off..  Galeta  (s.  Panama  (  anal:  PWO.  Adak  AK.  PWO.  Ed/ell  Scotland.  PWO. 

Puerto  Rico;  PWO.  Torri  Sta.  Okinawa;  Security  Offr.  Winter  Harbor  ME 
NAVSECSTA  PWD  -  Engr  Div.  Wash  .  DC 
NAVSHIPREPFAC  Library.  Guam;  SCE  Subic  Bay 

NAVSHIPYD  (Code  453)  Shop  03  Portsmouth.  VA;  Bremerton.  WA  (Carr  Inlet  Acoustic  Range):  Code  134. 
Pearl  Harbor.  HI;  Code  303.4.  Long  Beach  CA;  Code  303.5  (Library)  Puget  Sound.  Bremerton  WA;  Code 
3X0.  Mare  Is..  Vallejo.  CA;  Code  3X0. 38  (Goodwin).  Vallejo.  CA:  Code  380.  Portsmouth.  VA;  Code  383.3. 
Pearl  Harbor,  111:  Code  400.  Puget  Sound.  Code  410.  Mare  Is.,  Vallejo  CA;  Code  440  Portsmouth  NH; 

Code  440,  Norfolk;  Code  440.  Puget  Sound.  Bremerton  WA;  Code  453  (Util.  Supr),  Vallejo  CA. 

Commander.  Philadelphia.  PA;  I  D  Vivian:  Library.  Portsmouth  NH:  PW  Dept.  Long  Beach.  CA.  PWD 
(Code  430)  Dir  Portsmouth.  VA;  PWD  (Code  450-HD)  Portsmouth.  VA;  PWD  (Code  453-HD)  SHPO  03. 
Portsmouth.  VA;  PWD  (Code  457-HD)  Shop  07.  Portsmouth.  VA;  PWD  (Code  460)  Portsmouth.  VA: 

PWO.  Bremerton.  WA;  PWO,  Mare  Is.;  PWO.  Puget  Sound;  SCE.  Pearl  Harbor  HI;  Security  Offr. 
Bremerton  WA;  Security  Offr.  Charleston  SC;  Security  Offr.  Vallejo  CA;  Tech  Library.  Vallejo.  CA 
NAVSTA  Adak.  AK.  CO  Roosevelt  Roads  P  R.  Puerto  Rico;  CO.  Brooklyn  NY;  Code  4.  12  Marine  Corps 
Dist.  Treasure  Is..  San  Francisco  CA;  Dir  Engr  Div.  PWD.  Mavporl  FL;  Engr.  Dir..  Rota  Spain:  Long 
Beach.  CA;  Maint  Cont.  Div  .  Guantanamo  Bay  Cuba;  PWD  (LTJG.P.M.  Motolenich).  Puerto  Rico;  PWD 
-  Engr  Dept.  Adak.  AK;  PWD  -  Engr  Div,  Midway  Is.;  PWO.  Guantanamo  Bay  Cuba;  PWO.  Kcflavik 
Iceland.  PWO.  Mayport  FL;  SCE.  Guam;  SCE.  Pearl  Harbor  HI:  SCE.  San  Diego  CA;  SCE.  Subic  Bay. 

R  P  ;  Sec  Offr.  Roosevelt  Roads.  PR;  Security  Offr.  Guam;  Security  Offr.  San  Francisco.  CA;  Utilities  Engr 
Off  Rota  Spain 

NAVSUBASE  Code  23  (Slowev)  Bremerton.  WA;  SCE,  Pearl  Harbor  HI 
NAVSUPPACT  CO.  Naples.  Italy.  Engr  Div  (F.  Mollica).  Naples  Italy;  PWO  Naples  Italy 


NAVSUPPFAC  PWD  -  Main!  Control  Div.  Thurmont.  MD 
NAVSUPPO  PWO,  La  Maddalena.  Italy 

NAVSURFWPNCEN  G-52  (Duncan)  Dahlgrcn,  VA;  PWO,  White  Oak.  Silver  Spring.  MD;  Silver  Spring  MD 
NAVTECHTRACEN  SCE.  Pensacola  FL 

NAVTELCOMMCOM  Code  53,  Washington.  DC;  Security  Offr,  Washington  DC 
NAVUSEAWARENGSTA  Security  Spec  (Code  01A)  Key  >rt.  WA 
NAVWARCOL  President.  Newport.  R1 

NAVWPNCF.N  Code  24  (Dir  Safe  &  Sec)  China  Lake.  CA;  Code  2636  China  Lake;  Cmdr,  China  Lake.  CA; 
Code  26605  China  Lake  CA;  Code  3276  China  Lake.  CA;  Code  623  China  Lake  C'A;  PWO  (Code  266)  China 
Lake.  CA;  ROICC  (Code  702),  China  Lake  CA 

NAVWPNEVALFAC  CO  Kirtland  AFB.  NM;  Sec  Offr.  Kirtland  AFB.  NM;  Technical  Library.  Albuquerque 
NM 

NAVWPNSTA  (Clebak)  Colts  Neck.  NJ;  Code  04.  Concord.  CA;  Code  042.  Colts  Neck  NJ;  Code  042.  Concord 
CA;  Code  042A.  Seal  Beach.  CA;  Maint.  Control  Dir..  Yorktown  VA;  Nav  Wpns  Hndlng  Cen  (Code  8052) 
Colts  Neck  NJ;  Nuc  Safe  Offr.  Yorktown,  VA 
NAVWPNSTA  PW  Office  Yorktown.  VA 

NAVWPNSTA  PWD  -  Maint  Control  Div.  Charleston.  SC;  PWD  -  Maint.  Control  Div.,  Concord.  CA;  PWD  - 
Supr  Gen  Engr,  Seal  Beach.  CA;  PWO.  Charleston.  SC;  PWO.  Seal  Beach  CA;  Safety  Br.  Concord.  CA; 

Sec  Offr.  Charleston.  SC;  Security  Offr.  Colts  Neck  NJ;  Security  Offr.  Concord  CA;  Security  Offr.. 

Yorktown.  VA;  Special  Wpns  Offr.  Concord.  CA;  Wpns  Fac  Offr.  Colts  Neck,  NJ;  Wpns  Offr.  Charleston. 

SC:  Wpns  Offr.  Colts  Neck.  NJ;  Wpns  Offr.  Yorktown.  VA 
NAVWPNSUPPCEN  Code  04  Crane  IN 
NCBU  405  OIC.  San  Diego.  CA 
NCTC  Const.  Elec  School.  Port  Hueneme.  CA 

NCBC  Code  10  Davisville.  Rl:  Code  15.  Port  Hueneme  CA;  Code  155.  Port  Hueneme  CA;  Code  1552  (Brazele) 
Port  Hueneme.  CA;  Code  155B  (Nishimura)  Port  Hueneme.  CA;  Code  156.  Port  Hueneme.  CA;  Code  156F 
(Volpe)  Port  Hueneme.  CA:  Code  1571.  Port  Hueneme.  CA;  Code  25111  Port  Hueneme.  CA;  Code  400. 
Gulfport  MS:  Code  430  (PW  Engrng)  Gulfport.  MS;  Code  470.2.  Gulfport.  MS:  NEESA  Code  252  (P 
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